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Abstract 

Strongly correlated Fermi systems are fundamental systems in physics that are best studied experimentally, which until 
very recently have lacked theoretical explanations. This review discusses the construction of a theory and the analysis of 
phenomena occurring in strongly correlated Fermi systems such as heavy-fermion (HF) metals and two-dimensional (2D) 
Fermi systems. It is shown that the basic properties and the scaling behavior of HF metals can be described within the 
framework of a fermion condensation quantum phase transition (FCQPT) and extended quasiparticle paradigm that allow 
us to explain the non-Fermi liquid behavior observed in strongly correlated Fermi systems. In contrast to the Landau 
paradigm stating that the quasiparticle effective mass is a constant, the effective mass of new quasiparticles strongly 
depends on temperature, magnetic field, pressure, and other parameters. Having analyzed collected facts on strongly 
correlated Fermi systems with quite different microscopic nature, we find these to exhibit the same non-Fermi liquid 
behavior at FCQPT. We show both analytically and using arguments based entirely on the experimental grounds that the 
data collected on very different strongly correlated Fermi systems have a universal scaling behavior, and materials with 
strongly correlated fermions can unexpectedly be uniform in their diversity. Our analysis of strongly correlated systems 
such as HF metals and 2D Fermi systems is in the context of salient experimental results. Our calculations of the non- 
Fermi liquid behavior, the scales and thermodynamic, relaxation and transport properties are in good agreement with 
experimental facts. 
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1. Introduction 



Strongly correlated Fermi systems, such as heavy fermion (HF) metals, high-Tc superconductors, and two-dimensional 
(2D) Fermi liquids, are among the most intriguing and best experimentally studied fundamental systems in physics. 
However until very recently lacked theoretical explanations. The properties of these materials differ dramatically from 
those of ordinary Fermi systems UlSBBBSSISHIiS U . 121- For instance, in the case of metals with heavy fermions, 
the strong correlation of electrons leads to a renormalization of the effective mass of quasiparticles, which may exceed 
the ordinary, "bare", mass by several orders of magnitude or even become infinitely large. The effective mass strongly 
depends on the temperature, pressure, or applied magnetic field. Such metals exhibit NFL behavior and unusual power 
laws of the temperature dependence of the thermodynamic properties at low temperatures . Ideas based on quantum and 
thermal fluctuations taking place at a quantum critical point (QCP) have been put forward and the fascina ting behavior of 
these systems known as the non-Fermi liquid (NFL) behavior was attributed to the fluctuations HKslIiS, 3 15, 171. 
Suggested to describe one property, the ideas failed to do the same with the others and there was a real crisis and a new 
quantum phase transition responsible for the observed behavior was required 10. S& 13 11. 13, 181. 

The Landau theory of the Fermi liquid has a long history and remarkable results in describing a multitude of prop- 
erties of the electron liquid in ordinary metals and Fermi liquids of the ^He type | lS.2^,2iJ- The theory is based on 
the assumption that elementary excitations determine the physics at low temperatures. These excitations behave as quasi- 
particles, have a certain effective mass, and, judging by their basic properties, belong to the class of quasiparticles of a 
weakly interacting Fermi gas. Hence, the effective mass M* is independent of the temperature, pressure, and magnetic 
field strength and is a parameter of the theory. 

The Landau Fermi liquid (LFL) theory fails to explain the results of experimental observations related to the depen- 
dence of M* on the temperature T, magnetic field B, pressure, etc.; this has led to the conclusion that quasiparticles do 
not survive in stro ngly correlated Fermi systems and that the heavy electron does not retain its identity as a quasiparticle 
excitation Q H S E [IE E [H [H . 



1.1. Quantum phase transitions and the non-Fermi liquid behavior of correlated Fermi systems 

The unusual properties and NFL behavior observed in high-Tc superconductors, HF metals and 2D Fermi systems 
are assumed to be determined by various magnetic quantum phase transitions |[lll2.[3iS.|5lS0j[SS U. 12. 12. 14]. 
Since a quantum phase transition occurs at the temperature T - 0, the control parameters are the composition, electron 
(hole) number density x, pressure, magnetic field strength B, etc. A quantum phase transition occurs at a quantum critical 
point, which separates the ordered phase that emerges as a result of quantum phase transition from the disordered phase. 
It is usually assumed that magnetic (e.g., ferromagnetic and antiferromagnetic) quantum phase transitions are responsible 
for the NFL behavior. The critical point of such a phase transition can be shifted to absolute zero by varying the above 
parameters. 

Universal behavior can be expected only if the system under consideration is very close to a quantum critical point, 
e.g., when the correlation length is much longer than the microscopic length scale, and critical quantum and thermal fluc- 
tuations determine the anomalous contribution to the thermodynamic functions of the metal. Quantum phase transitions 
of this type are so widespread 



Tol U . 12 , 12 ] that we call them ordinary quantum phase transitions 1 22 1 ■ In this 



case, the physics of the phenomenon is determined by thermal and quantum fluctuations of the critical state, while quasi- 
particle excitations are destroyed by these fluctuations. Conventional arguments that quasiparticles in strongly correlated 
Fermi liquids "get heavy and die" at a quantum critical point commonly employ the well-known formula based on the 
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assumptions that the z-factor (the quasiparticle weight in the single-particle state) vanishes at the points of second-order 
phase transitions [18]. However, it has been shown that this scenario is problematic ll23li24ll . 

The fluctuations in the order parameter developing an infinite correlation and the absence of quasiparticle excitations 
is considered the main reason for the NFL behavior of heavy-fermion metals, 2D fermion systems and high-Tc super- 
conductors 12511 . This approach faces certain difficulties, however. Critical behavior in experiments with 
metals containing heavy fermions is observed at high temperatures comparable to the eff'ective Fermi temperature Tk- For 
instance, the thermal expansion coefficient a{T), which is a linear function of temperature for normal LFL, a(T) cc T, 
demonstrates the Vt temperature dependence in measurements involving CeNi2Ge2 as the temperature varies by two 
orders of magnitude (as it decreases from 6 K to at least 50 mK) 1 14]. Such behavior can hardly be explained within 
the framework of the critical point fluctuation theory. Obviously, such a situation is possible only as T — > 0, when the 
critical fluctuations make the leading contribution to the entropy and when the correlation length is much longer than the 
microscopic length scale. At a certain temperature T^, this macroscopically large correlation length must be destroyed by 
ordinary thermal fluctuations and the corresponding universal behavior must disappear. 

Another difficulty is in explaining the restoration of the LFL behavior under the application of magnetic field B, as 
observed in HF metals and in high-Tc superconductors fH [isl |26] . For the LFL state as T — > 0, the electric resistivity 
p(T) - pq +AT^, the heat capacity C{T) - joT, and the magnetic susceptibility - const. It turns out that the coefficient 
A(B), the Sommerfeld coefficient yoiB) cc M*, and the magnetic susceptibility ;t'(B) depend on the magnetic field strength 
B such that A(B) oc yliB) and A(B) oc x^(B), which implies that the Kadowaki-Woods relation K - A(B)/yl(B) |27] is 
B-independent and is preserved ll5ll . Such universal behavior, quite natural when quasiparticles with the effective mass 
M* playing the main role, can hardly be explained within the framework of the approach that presupposes the absence 
of quasiparticles, which is characteristic of ordinary quantum phase transitions in the vicinity of QCP. Indeed, there is no 
reason to expect that yo, x and A are affected by the fluctuations in a correlated fashion. 

For instance, the Kadowaki-Woods relation does not agree with the spin density wave scenario [15] and with the 
results of research in quantum criticality based on the renormalization-group approach [28]. Moreover, measurements 
of charge and heat transfer have shown that the Wiedemann-Franz law holds in some high-Tc superconductors [2^ 29ll 
and HF metals llsol Isil [32L Issll . All this suggests that quasiparticles do exist in such metals, and this conclusion is also 



corroborated by photoemission spectroscopy results 



The inability to explain the behavior of heavy-fermion metals while staying within the framework of theories based 
on ordinary quantum phase transitions implies that another important concept introduced by Landau, the order parameter, 
also ceases to operate (e.g., see Refsjainiliilil). Thus, we are left without the most fundamental principles of 
many -body quantum physics lfl9l 20l 21IL and many interesting phenomena associated with the NFL behavior of strongly 



correlated Fermi systems remain unexplained. 

NFL behavior manifests itself in the power-law behavior of the physical quantities of strongly correlated Fermi sys- 
tems located close to their QCPs, with exponents different from those of a Fermi liquid 1 36l 37]. It is common belief that 
the main output of theory is the explanation of these exponents which are at least depended on the magnetic character of 
QCP and dimensionality of the system. On the other hand, the NFL behavior cannot be captured by these exponents as 
seen from Fig. [T] Indeed, the specific heat CjT exhibits a behavior that is to be described as a function of both temperature 
r and magnetic B field rather than by a single exponent. One can see that at low temperatures CjT demonstrates the LFL 
behavior which is changed by the transition regime at which CjT reaches its maximum and finally CjT decays into NFL 
behavior as a function of T at fixed B. It is clearly seen from Fig. [T]that, in particularly in the transition regime, these 
exponents may have little physical significance. 

In order to show that the behavior of CjT displayed in Fig. [T]is of generic character, we remember that in the vicinity 
of QCP it is helpful to use "internal" scales to measure the effective mass M* cc CjT and temperature T |38, 3^. As 
seen from Fig. [1] a maximum structure in CjT cc M*^ at temperature Tm appears under the application of magnetic field 
B and Tm shifts to higher T as B is increased. The value of the Sommerfeld coefficient C/T = yo is saturated towards 
lower temperatures decreasing at elevated magnetic field. To obtain the normalized effective mass M*^, we use MJ^ and 
Tm as "internal" scales: The maximum structure in C/T was used to normalize C/T, and T was normalized by Tm- 
In Fig. |2]the obtained - M* /M*^ as a function of normalized temperature T^ - T /Tm is shown by geometrical 
figures. Note that we have excluded the experimental data taken in magnetic field B - 0.06 T. In that case, as will be 
shown in Subsections 19. 31 and |9.4. 51 Tm — » and the corresponding Tm and MJ|^ are unavailable. It is seen that the LFL 
state and NFL one are separated by the transition regime at which M*^ reaches its maximum value. Figure |2] reveals the 
scaling behavior of the normalized experimental curves - the curves at different magnetic fields B merge into a single 
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Figure 1 : Electronic specific lieat of YbRh2Si2, C/T, versus temperature T as a function of magnetic field B 13611 shown in the legend. 
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Figure 2: The normalized effective mass versus normalized temperature T^- M'^ is extracted from the measurements of the specific heat C/T on 
YbRh2Si2 in magnetic fields B 11361 listed in the legend. Constant effective mass M'^ inherent in normal Landau Fermi liquids is depicted by the solid 
line. 



one in terms of the normalized variable y = T/Tm- As seen from Fig. |2] the normalized effective mass M*^(y) extracted 
from the measurements is not a constant, as would be for a LFL, and shows the scaling behavior over three decades in 
normalized temperature y. It is seen from Figs. [T]and|2]that the NFL behavior and the associated scaling extend at least 
to temperatures up to few Kelvins. Scenario where fluctuations in the order parameter of an infinite (or sufficiently large) 
correlation length and an infinite correlation time (or sufficiently large) develop the NFL behavior can hardly match up 
such high temperatures. 

Thus, we conclude that a challenging problem for theories considering the critical behavior of the HF metals is to 
explain the scaling behavior of M^(y). While the theories calculating only the exponents that characterize M*^(y) at y » 1 

5 



deal with a part of the observed facts related to the problem and overlook, for example, consideration of the transition 
regime. Another part of the problem is the remarkably large temperature ranges over which the NFL behavior is observed. 

As we will see below, the large temperature ranges are precursors of new quasiparticles, and it is the scaling behavior 
of the normalized effective mass that allows us to explain the thermodynamic, transport and relaxation properties of HF 
metals at the transition and NFL regimes. 

Taking into account the simple behavior shown in Fig. |2l we ask the question: what theoretical concepts can replace 
the Fermi-liquid paradigm with the notion of the effective mass in cases where Fermi-liquid theory breaks down? To date 
such a concept is not available |3]. Therefore, in our review we focus on a concept of fermion condensation quantum 
phase transition (FCQPT) preserving quasiparticles and intimately related to the unlimited growth of M*. We shall show 
that it is capable revealing the scaling behavior of the effective mass and delivering an adequate theoretical explanation 
of a vast majority of experimental results in different HF metals. In contrast to the Landau paradigm based on the 
assumption that M* is a constant as shown by the solid line in Fig. |2l in FCQPT approach the effective mass M* of 
new quasiparticles strongly depends on T, x, B etc. Therefore, in accord with numerous experimental facts the extended 
quasiparticles paradigm is to be introduced. The main point here is that the well-defined quasiparticles determine as before 
the thermodynamic, relaxation and transport properties of strongly correlated Fermi-systems in large temperature ranges 
(see Section|9]and Subsection l9.4b . while M* becomes a function of T, x, B etc. The FCQPT approach had been already 
successfully applied to describe the thermodynamic properties of such different strongly correlated systems as ^He on the 
one hand and complicated heavy-fermion (HF) compounds on the other 23, 40(1 . 



1.2. Limits and goals of the review 

The purpose of this review is to show that diverse strongly correlated Fermi systems such three dimensional (3D) 
and 2D compounds as HF metals and 2D strongly correlated Fermi liquids exhibit a scaling behavior, which can be 
described within a single approach based on FCQPT theory 41, 42]. We discuss the construction of the theory and 



show that it delivers theoretical explanations of the vast majority of experimental results in strongly correlated systems 
such as HF metals and 2D systems. Our analysis is in the context of salient experimental results. Our calculations of the 
non-Fermi Uquid behavior, the scales and thermodynamic, relaxation and transport properties are in good agreement with 
experimental facts. We shall also focus on the scaling behavior of the thermodynamic, transport and relaxation properties 
that can be revealed from experimental facts and theoretical analysis. As a result, we do not discuss the specific features 
of strongly correlated systems in full; instead, we focus on the universal behavior of such systems. For instance, we ignore 
the physics of Fermi systems such as neutron stars, atomic clusters and nuclei, quark plasma, and ultra-cold gases in traps, 
in which we believe fermion condensate (FC) induced by FCQPT can exist [43, 44, 45, 46, 47, 48] . Ultra-cold gases in 
traps are interesting because their easy tuning allows selecting the values of the parameters required for observations of 
quantum critical point and FC. We do not discuss also microscopic mechanisms of quantum criticality related to FCQPT. 
Such mechanisms can be developed within FC theory. For example, the mechanism of quantum criticality as observed in 
f-electron materials can take place in systems when the centers of merged single-particle levels "get stuck" at the Fermi 
surface. One observes that this could provide a simple mechanism for pinning narrow bands in solids to the Fermi surface 
On the other hand, we consider high-Tc superconductors within a coarse-grained model based on the FCQPT theory 
in order to illuminate their generic relationships with HF metals. 

Experimental studies of the properties of quantum phase transitions and their critical points are very important for un- 
derstanding the physical nature of high-Tc superconductivity and HF metals. The experimental data that refer to different 
strongly correlated Fermi systems complement each other. In the case of high-Tf superconductors, only few experiments 
dealing with their QCPs have been conducted, because the respective QCPs are in the superconductivity range at low 
temperatures and the physical properties of the respective quantum phase transition are altered by the superconductivity. 
As a result, high magnetic fields are needed to destroy the superconducting state. But such experiments can be conducted 
for HF metals. Experimental research has provided data on the behavior of HF metals, shedding light on the nature of 



critical points and phase transitions (e.g., see Refs |15, 26, 29, 31, 32, 13413511 ). Hence, a key issue is the simultaneous 
study of high-Tf superconductors and the NFL behavior of HF metals. 

Since we are concentrated on properties that are non-sensitive to the detailed structure of the system we avoid difficul- 
ties associated with the anisotropy generated by the crystal lattice of solids, its special features, defects, etc.. We study the 
universal behavior of high-Tc superconductors, HF metals, and 2D Fermi systems at low temperatures using the model of 
a homogeneous HF liquid ll38i[39ll . The model is quite meaningful because we consider the scaling behavior exhibited 
by these materials at low temperatures, a behavior related to the scaling of quantities such as the effective mass, the heat 
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capacity, the thermal expansion, etc. The scahng properties of the normalized effective mass that characterizes them, are 
determined by momentum transfers that are small compared to momenta of the order of the reciprocal lattice length. The 
high momentum contributions can therefore be ignored by substituting the lattice for the jelly model. While the values of 
the scales like the maximum M*^ of the effective mass and Tm at which M*^ takes place are determined by a wide range 
of momenta and thus these scales are controlled by the specific properties of the system. 

We analyze the universal properties of strongly correlated Fermi systems using the FCQPT theory 10, 41, 42, 49ll . 
because the behavior of heavy-fermion metals akeady suggests that their unusual properties can be associated with the 
quantum phase transition related to the unlimited increase in the effective mass at the critical point. Moreover, we shall 
see that the scaling behavior displayed in Fig. |2]can be quite naturally captured within the framework of the quasiparticle 
extended paradigm supported by FCQPT which gives explanations of the NFL behavior observed in strongly correlated 
Fermi systems. 



2. Landau theory of Fermi liquids 

One of the most complex problems of modern condensed matter physics is the problem of the structure and properties 
of Fermi systems with large inter particle coupling constants. Theory of Fermi liquids, later called "normal", was first 
proposed by Landau as a means for solving such problems by introducing the concept of quasiparticles and amplitudes that 
characterize the effective quasiparticle interaction fl9','2^. The Landau theory can be regarded as an effective low-energy 
theory with the high-energy degrees of freedom eliminated by introducing amplitudes that determine the quasiparticle 
interaction instead of the strong inter particle interaction. The stability of the ground state of the Landau Fermi liquid is 
determined by the Pomeranchuk stability conditions: stability is violated when at least one Landau amplitude becomes 
negative and reaches its critical value [20, 50]. We note that the new phase in which stability is restored can also be 
described, in principle, by the LFL theory. 

We begin by recalling the main ideas of the LFL theory 12111 . The theory is based on the quasiparticle paradigm. 



which states that quasiparticles are elementary weakly excited states of Fermi liquids and are therefore specific excitations 
that determine the low-temperature thermodynamic and transport properties of Fermi liquids. In the case of the electron 
liquid, the quasiparticles are characterized by the electron quantum numbers and the effective mass M* . The ground state 
energy of the system is a functional of the quasiparticle occupation numbers (or the quasiparticle distribution function) 
n(p, T), and the same is true of the free energy F(n(p, T)), the entropy S (n(p, T)), and other thermodynamic functions. 
We can find the distribution function from the minimum condition for the free energy F - E - TS (here and in what 
follows ks-h- 1) 

6{F-iiN) 1 - «(P, T) 

——— ^sip,T)-^iT)-T\n — — — =0. (1) 
5n(p, T) n(p, T) 

Here /i is the chemical potential fixing the number density 



and 

6E(n{p, T)) 

e:(P,T) = (3) 
Snip, T) 

is the quasiparticle energy. This energy is a functional of «(p, T), in the same way as the energy E is: £(p, T, n). The 
entropy 5(«(p, T)) related to quasiparticles is given by the well-known expression ifloiEoll 



= -2 J [«(p. 



5 («(p, r)) = -2 [«(p, D ln(«(p, D) + (1 - «(p, T)) 



c/p 



X ln(l-n(p,r))]7:;3;3, (4) 



which follows from combinatorial reasoning. Equation ([T]) is usually written in the standard form of the Fermi-Dirac 
distribution, 

-1 



«(p, J) = < 1 -I- exp 



(e(p, T)-^i) 
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(5) 



At r ^ 0, dill and Q have the standard solution n{p, T — > 0) — > 6{pf - p) if the derivative deip ^ pf)/dp is finite 
and positive. Here pf is the Fermi momentum and 0(pf - p) is the step function. The single particle energy can be 
approximated as E{p ^ pf)- fi ^ Pf(p - Pf)/M*^, and inversely proportional to the derivative is the eff'ective mass of 
the Landau quasiparticle, 

1 1 ds{p,T = Q) ^ 

777 7 'p=pf ■ W 

Ml p dp 

In turn, the effective mass is related to the bare electron mass m by the well-known Landau equation | 2^, 21 ] 



^fr,(r,(PF,Pl) 



dp, (5^- 

where F^r.^r, (Pf, Pi) is the Landau amplitude, which depends on the momenta pp and p and the spins cr. For simplicity, 
we ignore the spin dependence of the effective mass, because M*^ is almost completely spin-independent in the case of a 
homogeneous liquid and weak magnetic fields. The Landau amplitude F is given by 

F,T,,Ti (p, Pi , n) = ^ , , , -— . (8) 

The stability of the ground state of LFL is determined by the Pomeranchuk stabili ty cond itions: stability is violated when 
at least one Landau amplitude becomes negative and reaches its critical value MM 



F"^' = -(2L + 1). (9) 

Here F"^ and F'l are the dimensionless spin-symmetric and spin-antisymmetric Landau amplitudes, L is the angular mo- 
mentum related to the corresponding Legendre polynomials P^, 

Fipcr, Pi<Ti)^-Yu t^i^' ^1 + ■ (10) 

L=0 

Here is the angle between momenta p and pi and the density of states = M^pf l(2jp-). It follows from Eq. (|7|l that 

m; f\ 

— = 1 + ^. (11) 

m i 

In accordance with the Pomeranchuk stability conditions it is seen from Eq. ( fTTT i that > -3, otherwise the effective 
mass becomes negative leading to unstable state when it is energetically favorable to excite quasiparticles near the Fermi 
surface. In what follows, we shall omit the spin indices cr for simplicity. 

To deal with the transport properties of Fermi systems, one needs a transport equation describing slowly varying 
disturbances of the quasiparticle distribution function np(r, f) which depends on position r and time t. As long as the 
transferred energy cj and momentum q of the quanta of external field are much smaller than the energy and momentum 
of the quasiparticles, qpp jij M*j) <s: 1 and co/T <s; 1, the quasiparticle distribution function n(q, oj) satisfies the transport 
equation ill IIS im 

drip 

— + VpSpVrnp - VrSpVpHp = /[Wp]. (12) 

The left-hand side of Eq. ( fT2l i describes the dissipationless dynamic of quasiparticles in phase space. The quasiparticle 
energy ep(r, f) now depends on its position and time, and the collision integral /[wp] measures the rate of change of the 
distribution function due to collisions. The transport equation (fT2l l allows one to derive all the transport properties and 
collective excitations of a Fermi system. 

It is common belief that the equations of this subsection are phenomenological and inapplicable to describe Fermi 
systems characterized by the effective mass M* strongly dependent on temperature, external magnetic fields B, pressure 
P etc. On the other hand, facts collected on HE metals demonstrate the specific behavior when the effective mass strongly 
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depends on temperature T, doping (or the number density) x and applied magnetic fields B, while the effective mass M* 
itself can reach very high values or even diverge, see e.g. 13l|41]. As we have seen in Section[T]such a behavior is so unusual 
that the traditional Landau quasiparticles paradigm fails to describe it. Therefore, in accord with numerous experimental 
facts the extended quasiparticles paradigm is to be introduced with the well-defined quasiparticles determining as before 
the thermodynamic and transport properties of strongly correlated Fermi-systems, M* becomes a function of T, x, B, 
while the dependence of the effective mass on T, x, B gives rise to the NFL behavior 

As we shall see in the following Section[3] Eq. O can be derived microscopically and it becomes compatible with the 
extended paradigm. 



3. Equation for the effective mass and the scaling behavior 

To derive the equation determining the effective mass, we consider the model of a homogeneous HF liquid and employ 
the density functional theory for superconductors (SCDFT) Ii54] which allows us to consider £ as a functional of the 
occupations numbers n(p) 1381 pSi pg pTI] . As a result, the ground state energy of the normal state E becomes the 
functional of the occupation numbers and the function of the number density x, E - £'(n(p), x), while Eq. (O gives 
the single-particle spectrum. Upon differentiating both sides of Eq. (O with respect to p and after some algebra and 



integration by parts, we obtain ll23l[38ll55il5{ 



5e(p) p r ,dn(pi) dp 



m J 



. F(p,pi,«) ■ (13) 

op m J dp I {Iny 



To calculate the derivative 5e(p)/5p, we employ the functional representation 

dp 



dpdpi 

(2^ 



F(p,pi,n)|„^„n(p)n(pi)^^ + ... (14) 



It is seen directly from Eq. ( fTlT l that the effective mass is given by the well-known Landau equation 



----- r 

M* m J 



PfPi j,, .dn(pi) dpi 

-^F(pF,pun)— —3. (15) 

Pp dpi (2ny 



For simplicity, we ignore the spin dependencies. To calculate M* as a function of T, we construct the free energy 
F = E - TS , where the entropy S is given by Eq. Minimizing F with respect to n(p), we arrive at the Fermi-Dirac 
distribution, Eq. (|5]l. Due to the above derivation, we conclude that Eqs. ( fTlT l and ( fTSl l are exact ones and allow us to 
calculate the behavior of both ds{p)/dp and M* which now is a function of temperature T, external magnetic field B, 
number density x and pressure P rather than a constant. As we will see it is this feature of M* that forms both the scaling 
and the NFL behavior observed in measurements on HF metals. 

In LFL theory it is assumed that M*^ is positive, finite and constant. As a result, the temperature-dependent corrections 
to M^, the quasiparticle energy e(p) and other quantities begin with the term proportional to in 3D systems and with 
the term proportional to T in 2D one llssll . The effective mass is given by Eq. (|7]i, and the specific heat C is lfl9ll 



and the spin susceptibility 



^ In^NT dS 

c.^=r„r = r-, (16) 

^-^^^iTTF^y ^''^ 

where jUb is the Bohr magneton and yo cc M*^. In the case of LFL, upon using the transport Eq. (fT2l i one finds for the 
electrical resistivity at low T yl | 

p(T) ^ po + AT''^ (18) 
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where po is the residual resistivity, the exponent - 2 and A is the coefficient determining the charge transport. The 
coefficient is proportional to the quasiparticle-quasiparticle scattering cross-section. Equation ( fTST i symbolizes and defines 
ffie LFL behavior observed in normal metals. 

Equation ( fTsT i at T = 0, combined with the fact that n(p, T - 0) becomes 6(pf - p), yields the well-known result 
'|60l|6lll 

M* 1 



1 -FV3 



where - Nq/^ Nq = mpf/{2n^) is the density of states of a free Fermi gas and f^{pF,PF) is the p-wave component 
of the Landau interaction amplitude. Because x = /5^/37r^ in the Landau Fermi-liquid theory, the Landau interaction 
amplitude can be written as F^{pf,pf) - F^{x). Provided that at a certain critical point jc^Cj the denominator (1 -F'(x)/3) 
tends to zero, i.e., (1 - F^{x)l3) oc (x- xpc) + o.{x - xpc)^ + ... — > 0, we find ffiat ||62, 63 1 



M*(x) 



■ a\ + 



fl2 



(19) 



where a\ and 02 are constants and r - (x- Xfc)/xfc is the "distance" from QCP Xfc at which M*{x Xfc) — > 00. We 
note that the divergence of the effective mass given by Eq. ( fT9] l does preserve the Pomeranchuk stability conditions for 
positive, see Eq. (|9]). Equations (fTTT i and ST% seem to be different but it is not the case since F^ oc m, while F^ oc M* 
and Eq. (fTTT i represents an implicit formula for the effective mass. 




x(10""cm"') 



Figure 3: The ratio M*/M in a silicon MOSFET as a function of the electron number density x. The black squares mark the experimental data on 
the Shubnikov-de Haas oscillations. The data obtained by applying a parallel magnetic field are marked by black circles L65» ,66t,71.1 . The solid line 
represents the function (86). 



The behavior of M*(x) described by formula (fT9l l is in good agreement with the results of experiments 11641 l65l l66l] 
and calculations |67, 68, 69]. In the case of electron systems, Eq. (fT9] l holds for x > Xfc, while for 2D ^He we have 
X < Xfc so that always r > ll42i. 70,1 (see also Section|S]l. Such behavior of the effective mass is observed in HE metals, 
which have a fairly flat and narrow conductivity band corresponding to a large effective mass, with a strong correlation 
and the effective Fermi temperature Tk ~ p^/M*{x) of the order of several dozen degrees kelvin or even lower (e.g., see 
Ref. lE]). 

The effective mass as a function of the electron density x in a silicon MOSFET (Metal Oxide Semiconductor Field 
Effect Transistor), approximated by Eq. ( fT9] l. is shown in Fig. |3] The parameters oi, 02 and Xfc are taken as fitting. We 
see that Eq. (fT9] l provides a good description of the experimental results. 
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Figure 4: The ratio M' /M in 2D ^He as a function of the density x of the liquid, obtained from heat capacity and magnetization measurements. The 
experimental data are marked by black squares l64ll73l . and the solid line represents the function given by Eq. jl9t . where ai = 1.09, = 1.68 nm"^, 
and Xfc = 5.11 nm"^. 



The divergence of the effective mass M*{x) discovered in measurements involving 2D ^He 164 17211 is illustrated in 
Fig. |4] Figures [3] and ID show that the description provided by Eq. ( fT9] ) does not depend on elementary Fermi particles 
constituting the system and is in good agreement with the experimental data. 

It is instructive to briefly explore the scaling behavior of M* in order to illustrate the ability of the quasiparticle 
extended paradigm to capture the scaling behavior, while more detailed consideration is reserved for Section |9l Let us 
write the quasiparticle distribution function as ni(p) = n(p, T) - n(p), with n(p) being the step function, and Eq. (fTSl l then 
becomes 



M* J 



M*(T) M* J opi (zny 

At QCP X — > Xfc, the effective mass M*{x) diverges and Eq. (l20l i becomes homogeneous determining M* as a function of 
temperature while the system exhibits the NFL behavior. If the system is located before QCP, M* is finite, at low temper- 
atures the integral on the right hand side of Eq. ( l20l l represents a small correction to 1 /M* and the system demonstrates 
the LFL behavior seen in Figs. [T]and|2] The LFL behavior assumes that the effective mass is independent of temperature, 
M*{T) ^ const, as shown by the horizontal line in Fig. |2] Obviously, the LFL behavior takes place only if the second term 
on the right hand side of Eq. (|20] | is small in comparison with the first one. Then, as temperature rises the system enters 
the transition regime: M* grows, reaching its maximum M*^ RtT - Tm, with subsequent diminishing. As seen from Fig. 
|2] near temperatures T > Tm the last "traces" of LFL regime disappear, the second term starts to dominate, and again Eq. 
( l20b becomes homogeneous, and the NFL behavior is restored, manifesting itself in decreasing M* as a function of T. 



4. Fermion condensation quantum phase transition 

As shown in Section [3] the Pomeranchuk stability conditions do not encompass all possible types of instabilities 
and that at least one related to the divergence of the effective mass given by Eq. (1% was overlooked |41]. This type 
of instability corresponds to a situation where the effective mass, the most important characteristic of quasiparticles, 
can become infinitely large. As a result, the quasiparticle kinetic energy is infinitely small near the Fermi surface and 
the quasiparticle distribution function «(p) minimizing £'(n(p)) is determined by the potential energy. This leads to the 
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formation of a new class of strongly correlated Fermi liquids with FC ll4lll42il48i[73ll . separated from the normal Fermi 
liquid by FCQPT |22','74l|75|]. 

It follows from ( fT9l ) that atT - and as r the effective mass diverges, M*{r) — > oo. Beyond the critical point Xfc, 
the distance r becomes negative and, correspondingly, so does the effective mass. To avoid an unstable and physically 
meaningless state with a negative effective mass, the system must undergo a quantum phase transition at the critical point 
X - Xfc, which, as we will see shortly, is FCQPT [174 l75l 122 1. Because the kinetic energy of quasiparticles that are near 
the Fermi surface is proportional to the inverse effective mass, the potential energy of the quasiparticles near the Fermi 
surface determines the ground-state energy as x — > Xfc- Hence, a phase transition reduces the energy of the system and 
transforms the quasiparticle distribution function. Beyond QCP x = Xfc, the quasiparticle distribution is determined by 
the ordinary equation for a minimum of the energy functional lij]]: 

= e(p) - u. Pi < p < Pf. (21) 

Equation (ISTT i yields the quasiparticle distribution function «o(p) that minimizes the ground-state energy E. This 
function found from Eq. ( |2T] ) differs from the step function in the interval from pi to pj, where < no(p) < 1, and 
coincides with the step function outside this interval. In fact, Eq. (1211 1 coincides with Eq. (O provided that the Fermi 
surface ai p - pf transforms into the Fermi volume at pi < p < pf suggesting that the single-particle spectrum is 
absolutely "flat" within this interval. A possible solution no(p) of Eq. (|2TI) and the corresponding single-particle spectrum 
e(p) are depicted in Fig. |5] Quasiparticles with momenta within the interval {pf - pi) have the same single-particle 
energies equal to the chemical potential ji and form FC, while the distribution no(p) describes the new state of the Fermi 
liquid with FC l,4L,4Z ,73,1 . In contrast to the Landau, marginal, or Luttinger Fermi liquids \^J6cJT\, which exhibit the 




Figure 5: The single-particle spectrum e(p) and the quasiparticle distribution function nQ{p). Because no(p) is a solution of Eq. i2U . we have 
no(p < Pi) = 1, < noipi < p < Pf) < 1, and no(p > pf) = 0, while E(pi < p < pf) = p. The Fermi momentum pf satisfies the condition 
Pi <Pf < Pf- 

same topological structure of the Green's function, in systems with FC, where the Fermi surface spreads into a strip, the 
Green's function belongs to a different topological class. The topological class of the Fermi liquid is characterized by the 
invariant 1,46- .47, .73,1 

N = tr(£ -^Giioj, p)5/G"' (/w, p), (22) 



where "tr" denotes the trace over the spin indices of the Green's function and the integral is taken along an arbitrary 
contour C encircling the singularity of the Green's function. The invariant in ( l22l i takes integer values even when the 
singularity is not of the pole type, cannot vary continuously, and is conserved in a transition from the Landau Fermi liquid 



to marginal liquids and under small perturbations of the Green's function. As shown by Volovik 11461 1471I73I1 . the situation 
is quite different for systems with FC, where the invariant becomes a half-integer and the system with FC transforms 
into an entirely new class of Fermi liquids with its own topological structure. 

4.1. The order parameter ofFCQPT 

We start with visualizing the main properties of FCQPT. To this end, again consider SCDFT. SCDFT states that 
the thermodynamic potential <£ is a universal functional of the number density n(r) and the anomalous density (or the 
order parameter) K(r, ri), providing a variational principle to determine the densities. At the superconducting transition 
temperature Tc a superconducting state undergoes the second order phase transition. Our goal now is to construct a 
quantum phase transition which evolves from the superconducting one. 



Let us assume that the coupling constant /Iq of the BCS-Uke pairing interaction 117811 vanishes, with Aq — » making 
vanish the superconducting gap at any finite temperature. In that case, Tc ^ Q and the superconducting state takes place 
at r = while at finite temperatures there is a normal state. This means that at J = the anomalous density 

^(r,ri) = <^T(r)^'i(ri)> (23) 

is finite, while the superconducting gap 

A(r) = ^o I K(r,ri)dri (24) 



is infinitely small i6l l63h . In Eq. ( l23T l. the field operator ^^(^(r) annihilates an electron of spin cr, cr =f > i at the position 
r. For the sake of simplicity, we consider the model of homogeneous HF liquid |6]. Then at T = 0, the thermodynamic 
potential reduces to the ground state energy E which turns out to be a functional of the occupation number «(p) since 
in that case the order parameter k(p) - v(p)M(p) - V«(p)(l - n(p)). Indeed, 

«(p) = y2(p); ^(p) = y(p)«(p), (25) 



where u(p) and v(p) are normalized parameters such that v^(p) + M^(p) - 1 and K(p) - V«(p)(l - n(p)), see e.g. ||2C 

Upon minimizing E with respect to «(p), we obtain Eq. (|2TI) . As soon as Eq. (1211 1 has nontrivial solution no(p) then 
instead of the Fermi step, we have < n()(p) < 1 in certain range of momenta pi < p < pj with K{p) - V«o(p)(l - no(p)) 
being finite in this range, while the single particle spectrum £(p) is flat. Thus, the step-like Fermi filling inevitably 
undergoes restructuring and forms FC when Eq. (1211 1 possesses for the first time the nontrivial solution at x = Xc which is 
QCP ofFCQPT. In that case, the range vanishes, /?,■ — > Pf ^ Pf, and the effective mass M* diverges at QCP 16, 23, 41 1 



M%x ^ xc) ~ pF dp 1"-"-^-'--' ^ 

At any small but finite temperature the anomalous density k (or the order parameter) decays and this state undergoes the 
first order phase transition and converts into a normal state characterized by the thermodynamic potential Oq. Indeed, 
at r — > 0, the entropy S = -d<S>o/dT of the normal state is given by Eg. dUl. It is seen from Eq. ^ that the normal 
state is characterized by the temperature-independent entropy 5o ||6, 23, [tH. Since the entropy of the superconducting 



ground state is zero, we conclude that the entropy is discontinuous at the phase transition point, with its discontinuity 
6S = Sq. Thus, the system undergoes the first order phase transition. The heat q of transition from the asymmetrical to 
the symmetrical phase is ^ = TcSo - since T^ = 0. Because of the stability condition at the point of the first order phase 
transition, we have <l)o(n(p)) - <l>(/('(p)). Obviously the condition is satisfied since ^ = 0. 

4.2. Quantum protectorate related to FCQPT 

With FCQPT (as well as with other phase transitions), we have to deal with strong particle interaction, and there is no 
way in which a theoretical investigation based on first principles can provide an absolutely reliable solution. Hence, the 
only way to verify that FC exists is to study this state by exactly solvable models and to examine the experimental facts 
that could be interpreted as direct confirmation of the existence of FC. Exactly solvable models unambiguously suggest 
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that Fermi systems with FC exist (e.g., see Refs 1180118 IL I821I83I1 ). Taking the results of topological investigations into 
account, we can state that the new class of Fermi liquids with FC is nonempty, actually exists, and represents an extended 



family of new states of Fermi systems l>46l I47ll73ll . 

We note that the solutions no(p) of Eq. (l2TT l are new solutions of the well-known equations of the Landau Fermi- 
liquid theory. Indeed, atT - Q, the standard solution given by a step function, «(p, T — > 0) — > 6(pf - p), is not the only 
possible one. Anomalous solutions e(p) = yu of Eq. ([TJ can exist if the logarithmic expression on its right-hand side is 
finite. This is possible if < no(p) < 1 within a certain interval (p,- < p < pf). Then, this logarithmic expression remains 
finite within this interval as T — > 0, the product T ln[(l - «o(p))/no(p)]|r^o — > 0, and we again arrive at Eq. ( 1211 1. 

Thus, as r — > 0, the quasiparticle distribution function no(p), which is a solution of Eq. (1211 1. does not tend to the step 
function 0{pf - p) and, correspondingly, in accordance with Eq. (|4|l, the entropy 5(7") of this state tends to a finite value 
5o as r — > 0: 

S{T^Q)^Sq. (27) 

As the density x — » xpc (or as the interaction force increases), the system reaches QCP at which FC is formed. This 
means that and that the deviation (5n(p) from the step function is small. Expanding the function £'(n(p)) 

in Taylor series in 5n(p) and keeping only the leading terms, we can use Eq. (|2TI) to obtain the following relation that is 
valid within the interval pi < p < pf 



ii = e(p) = eo(p) + I F(p, pi)6n(pi)-^. (28) 

(2nY 



Both quantities, the Landau amplitude F(p, pi) and the single-particle energy eo(p), are calculated at n(p) = 0(pf - p). 
Equation (l28T l has nontrivial solutions for densities x < Xfc if the corresponding Landau amplitude, which is density- 
dependent, is positive and sufficiently large for the potential energy to be higher than the kinetic energy. For instance, 
such a state is realized in a low-density electron liquid. The transformation of the Fermi step function n(p) - 6(pf - p) 



into a smooth function determined by Eq. (l28T l then becomes possible i4lll42l 

It follows from Eq. ( l28l l that the quasiparticles of FC form a collective state, because their state is determined by the 
macroscopic number of quasiparticles with momenta pi < p < pj. The shape of the single-particle spectrum related to 
FC is independent of the Landau interaction, which is in general determined by the properties of the system as a whole, 
including the collective states, irregularities of structure, the presence of impurities, and composition. The length of the 
interval from p, to where FC exists is the only characteristic determined by the Landau interaction; of course, the 
interaction must be strong enough for FCQPT to occur Therefore, we conclude that spectra related to FC have a universal 
shape. In Sections [4. Sl and lSTI we show that these spectra are dependent on the temperature and the superconducting gap 
and that this dependence is also universal. The existence of such spectra can be considered a characteristic feature of a 
"quantum protectorate", in which t he prop erties of the material, including the thermodynamic properties, are determined 



by a certain fundamental principle [84|, |85|] . In our case, the state of matter with FC is also a quantum protectorate, since 
the new type of quasiparticles of this state determines the special universal thermodynamic and transport properties of 
Fermi liquids with FC . 

4.3. The influence of FCQPT at flnite temperatures 

According to Eq. ([TJ, the single-particle energy e(p, T) is linear in T for T ^ Tf within the interval (pf - pd lH; 
Expanding ln((l - n(p))/n(p)) in a series in n(p) at /? ^ pf, we can write the expression 



e(p, T)-^l{T) ^ 1 - «(p) _ 1 - 2«(p) 



n(p) «(p) 

where T j is the temperature above which the effect of FC is insignificant 115 ill : 



(29) 

P~PF 



with Qfc being the volume occupied by FC, Sf being the Fermi energy, and being the volume of the Fermi sphere. 
We note that for T <k T f, the occupation numbers «(p) obtained from Eq. (ISTT i are almost perfectly independent of 
T l^^iS^. At finite temperatures, according to Eq. ( |29] l, the dispersionless plateau £(p) - y. shown in Fig. |5]is 
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slightly rotated counterclockwise in relation to ju. As a result, the plateau is slightly tilted and rounded off at its end points. 
According to Eqs. ^ and ( |29] l. the effective mass M^^ that refers to the FC quasiparticles is given by 



Kc-PF^^^- (31) 



In deriving ( ISTT i. we approximated the derivative as dn{p)/dp ^ -l/(pf - pi). Equation ( ISTT i clearly shows that for 
< T Tf, the electron liquid with FC behaves as if it were placed at a quantum critical point, since the electron 
effective mass diverges as T — > 0. Actually, as we shall see in Subsection 14.41 the system is at a quantum critical line, 
because critical behavior is observed behind QCP with x = Xfc of FCQPT as T 0. In Sections |7]and[T0] we show that 
the behavior of such a system differs dramatically from that of a system at a quantum critical point. 
Upon using Eqs. (l30t and (ISTT i. we estimate the effective mass M^^-, as 

M ~ A^o(O) ~ r ' 

where A^o(O) is the density of states of a noninteracting electron gas and A'(O) is the density of states on the Fermi surface. 
Equations (ISTT i and (l32l l yield the temperature dependence of M^^. 

Multiplying both sides of Eq. ( ISTT i by (/?/ - p,), we obtain an expression for the characteristic energy, 

£0 - 4r, (33) 

which determines the momentum interval (pf- pi) with the low-energy quasiparticles characterized by the energy |e(p) - 
H\ < Eo/2 and the effective mass M^^^. The quasiparticles that do not belong to this momentum interval have an energy 
|e(p) - fi\ > Eo/2 and an effective mass M*^ that is weakly temperature-dependent |74, 75, 87]. Equation ( l33T l shows that 
Eq is independent of the condensate volume. We conclude from Eqs. ( ISTT i and ( |33] ) that for T <K Tf, the single-electron 
spectrum of FC quasiparticles has a universal shape and has the features of a quantum protectorate. 

Thus, a system with FC is characterized by two effective masses, M^^ and M^. This fact manifests itself in a break or 
an abrupt change in the quasiparticle dispersion law, which for quasiparticles with energies e(p) < p can be approximated 
by two straight lines intersecting at £0/2 ^ 2T. Figure |5] shows that at T = 0, the straight hnes intersect at p = This 
break also occurs when the system is in its superconducting state at temperatures Tc < T «i Tf, where Tc is the critical 
temperature of the superconducting phase transition, which agrees with the experimental data in Issll and, as we will see 
in Section |5] this behavior agrees with the experimental data at T < T^. At T > Tc, the quasiparticles are well-defined. 



because their width y is small compared to their energy and is proportional to the temperature, y ~ rj|34l 15111 . The 
quasiparticle excitation curve (see Section [6]l can be approximately described by a simple Lorentzian |87j, which also 
agrees with the experimental data 1 88 , 8^ 93 91 ] . 



We estimate the density xpc at which FCQPT occurs. We show in Section|8]that an unlimited increase in the effective 
mass precedes the appearance of a density wave or a charge density wave formed in electron systems at = rcdw, where 
rs = ro/aB, ro is the average distance between electrons, and ag is the Bohr radius. Hence, FCQPT certainly occurs at 
T = when rj reaches its critical value rfc corresponding to Xfc, with rpc < ^cdw iItoIi . We note that the increase in the 
effective mass as the electron number density decreases was observed in experiments, see Figs. |3]and|4] 

Thus, the formation of FC can be considered a general property of different strongly correlated systems rather than 
an exotic phenomenon corresponding to the anomalous solution of Eq. (ISTT i. Beyond FCQPT, the condensate volume is 
proportional to (r, - rfc), with Tf/sp ~ (rs - rfc)/rfc, at least when (r, - rFc)/rFC ^ 1- This implies that 

rs - rpc ^ Pf - Pi ^ Xfc-x ^^^^ 
rpc Pf Xfc 

Because a state of a system with FC is highly degenerate, FCQPT serves as a stimulator of phase transitions that could lift 
the degeneracy of the spectrum. For instance, FC can stimulate the formation of spin density waves, antiferromagnetic 
state and ferromagnetic state etc., thus strongly stimulating the competition between phase transitions eliminating the 
degeneracy. The presence of FC strongly facilitates a transition to the superconducting state, because both phases have 
the same order parameter 
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4.4. Phase diagram of Fermi system with FCQPT 

At r = 0, a quantum phase transition is driven by a nonthermal control parameter, e.g. the number density x. As 
we have seen, at QCP, x - Xfc, the effective mass diverges. It follows from Eq. ( fT9l ) that beyond QCP, the effective 
mass becomes negative. To avoid an unstable and physically meaningless state with a negative effective mass, the system 
undergoes FCQPT leading to the formation of FC. 




Figure 6: Schematic phase diagram of system with FC. The number density x is taken as the control parameter and depicted as x/xfc- The dashed line 
shows M*(A-/.Vf(;) as the system approaches QCP,x/xfc = 1 , of FCQPT which is denoted by the arrow. Atx/xfc > 1 and sufficiently low temperatures, 
the system is in the LFL state as shown by the shadow area. At T = and beyond the critical point, x/xfc < 1, the system is at the quantum critical 
line depicted by the dashed line and shown by the vertical arrow. The critical line is characterized by the FC state with finite superconducting order 
parameter At any finite low temperature T > = 0, k is destroyed, the system undergoes the first order phase transition, possesses finite entropy So 
and exhibits the NFL behavior at any finite temperatures T < Tf. 



A schematic phase diagram of the system which is driven to the FC state by variation of x is reported in Fig. |6] 
Upon approaching the critical density Xfc the system remains in the LFL region at sufficiently low temperatures as it 
is shown by the shadow area. The temperature range of the shadow area shrinks as the system approaches QCP, and 
M*{x/xfc) diverges as shown by the dashed line and Eq. (fT9] l. At QCP Xfc shown by the arrow in Fig. |6] the system 
demonstrates the NFL behavior down to the lowest temperatures. Beyond the critical point at finite temperatures the 



behavior remains the NFL and is determined by the temperature-independent entropy So la [79|]. In that case at T — > 0, 
the system is approaching a quantum critical line (shown by the vertical arrow and the dashed line in Fig. |6]l rather than a 
quantum critical point. Upon reaching the quantum critical line from the above at T — > the system undergoes the first 
order quantum phase transition, which is FCQPT taking place at Tc - 0. While at diminishing temperatures, the systems 
located before QCP do not undergo a phase transition and their behavior transits from NFL to LFL. 

It is seen from Fig. |6]that at finite temperatures there is no boundary (or phase transition) between the states of 
systems located before or behind QCP shown by the arrow. Therefore, at elevated temperatures the properties of systems 
with x/xfc < 1 or with x/xfc > 1 become indistinguishable. On the other hand, at T > the NFL state above the critical 
line and in the vicinity of QCP is strongly degenerate, therefore the degeneracy stimulates the emergence of different 
phase transitions lifting it and the NFL state can be captured by the other states such as superconducting (for example, 
by the superconducting state (SC) in CeCoIng li63[ l79ll ) or by antiferromagnetic (AF) state (e.g. AF one in YbRh2Si2 
OSll ) etc. The diversity of phase transitions occurring at low temperatures is one of the most spectacular features of the 
physics of many HE metals. Within the scenario of ordinary quantum phase transitions, it is hard to understand why 
these transitions are so different from one another and their critical temperatures are so extremely small. However, such 
diversity is endemic to systems with a FC [23]. 

Upon using nonthermal tuning parameters like the number density, pressure or magnetic field, the NFL behavior is 
destroyed and the LFL one is restored as we shall see in Sections|9]and[T0] For example, the application of magnetic field 
B > Bed drives a system to QCP and destroys the AF state restoring the LFL behavior. Here, Bco is a critical magnetic 
field, such that atB > Bco the system is driven towards its LFL state. In some cases as in the HE metal CeRu2Si2, BcO = 0, 
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see e.g. while in YbRhaSii, B^o ^ 0.06 T yj. 



5. The superconducting state with FC 

In this section we discuss the superconducting state of a 2D hquid of heavy electrons, since high-Tc superconductors 
are represented mainly by 2D structures. On the other hand, our study can easily be generalized to the 3D case. To show 
that there is no fundamental difference between the 2D and 3D cases, we derive Green's functions for the 3D case in 
SectionlS^ 



5.1. The superconducting state at T — 

As we have seen in Subsection l4. II the ground-state energy Egs{K{p), n(p)) of a 2D electron liquid is a functional of the 
superconducting state order parameter A-(p) and of the quasiparticle occupation numbers n(p). This energy is determined 
by the well-known Bardeen-Cooper-SchriefFer (BCS) equations and in the weak-coupling superconductivity theory is 



given by |54l78l|93|] 



£j,Wp),n(p)) = E(n(p)) + Ao j VipuPi) 

X 4PiK(P2)^^. (35) 

It is assumed that the constant Aq, which determines the magnitude of the pairing interaction /loV(pi,p2), is small. We 
define the superconducting gap as 

A(p) = -^o J y(p,piMPi)^- (36) 
Minimizing Egs in v(p) and using (l36l l. we arrive at equations that relate the single-particle energy e(p) to A(p) and ^(p) 

1 - 2v\p) A(p) 

Here the single-particle energy e(p) is determined by Eq. (O, and 

E(p) = ^^2(p) + A2(p), (38) 

with ^(p) - e(p) - ju. Substituting the expression for k(p) from ( l37l i in Eq. ( l36b . we obtain the well-known equation of 
the BCS theory for A(p) 

As Aq 0, the maximum value of Ai of the superconducting gap A(p) tends to zero and each equation in (|37] ) reduces to 
Eq. (EB 

5£'(n(p)) , , „ 
^ , s = e(p) -H = 0, (40) 

if < n(p) < 1, or k(p) + 0, in the interval p, < p < pj. Equation ( l40l l shows that the function no(p) is determined 
from the solution to the standard problem of finding the minimum of the functional £'(n(p)) [41, 51, 52]. Equation ( l40t 
specifies the quasiparticle distribution function no(p) that ensures the minimum of the ground-state energy E(k(p), «(p)). 
We can now study the relation between the state specified by Eq. (|40] | or Eq. (1211 1 and the superconducting state. 

At r = 0, Eq. ( l40l i determines the specific state of a Fermi liquid with FC, the state for which the absolute value of 
the order parameter |/f(p)| is finite in the momentum interval p, < p < pf as Ai — » 0. Such a state can be considered 
superconducting with an infinitely small value of Ai. Hence, the entropy of this state at T = is zero. Solutions «o(p) 
of Eq. (l40l i constitute a new class of solutions of both the BCS equations and the Landau Fermi-liquid equations. In 
contrast to the ordinary solutions of the BCS equations iFtsIi . the new solutions are characterized by an infinitely small 
superconducting gap Ai — > 0, with the order parameter A-(p) remaining finite. On the other hand, in contrast to the standard 
solution of the Landau Fermi-liquid theory, the new solutions no(p) determine the state of a heavy-electron liquid with a 
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finite entropy Sq as T — > (see Eq. (|27])). We arrive at an important conclusion tliat the solutions of Eq. ( l40l i can be 
interpreted as the general solutions of the BCS equations and the Landau Fermi-liquid theory equations, while Eq. ( l40t 
itself can be derived either from the BCS theory or from the Landau Fermi-liquid theory. Thus, as shown in Subsection l4. II 
both states of the system coexist as T — > 0. As the system passes into a state with the order parameter k{j^), the entropy 
suddenly vanishes, with the system undergoing the first-order transition near which the critical quantum and thermal 
fluctuations are suppressed and the quasiparticles are well- defined excitations (see also Section [TOb. It follows from 
Eq. (|22] ) that FCQPT is related to a change in the topological structure of the Green's function and belongs to Lifshitz's 
topological phase transitions, which occur at absolute zero [73)]. This fact establishes a relation between FCQPT and 
quantum phase transitions under which the Fermi sphere splits into a sequence of Fermi layers 11941 19511 (see Sections [T] 
andfTSl). We note that in the state with the order parameter k(p), the system entropy 5=0 and the Nernst theorem holds 
in systems with FC. 

If 0, the gap Ai becomes finite, leading to a finite value of the effective mass M^^, which may be obtained from 
Eq. ( l37T i by taking the derivative with respect to the momentum p of both sides and using Eq. (|6| ll74l75ll87ll : 



Mpc - Pf 



Pf - Pi 
2Ai ■ 



(41) 



It follows from Eq. ( |4TI) that in the superconducting state the effective mass is always finite. As regards the energy scale, 
it is determined by the parameter Eq: 



Eo = e(p/) - e(p/) - PF^-^hr^ - 2Ai. 



(42) 



5.2. Green 's function of the superconducting state with FC at T — 

We write two equations for the 3D case, the Gor'kov equations ll96ll . which determine the Green's functions F^(p, to) 
and G(p, w) of a superconductor (e.g., see Ref. |20]): 



G 

The gap A and the function E are given by 



{(jj - £(p) + i 0)(w + £(p) - / 0) 

«'(P) ^ v'(p) 
Lo-E(p) + iO oj + E(p)-iO' 



A^Aom, /H= r rF-(p,^)^ 



dcodp 



(43) 



(44) 



We recall that the function F^(p, a>) has the meaning of the wave function of Cooper pairs and S is the wave function of 
the motion of these pairs as a whole and is just a constant in a homogeneous system |20]. It follows from Eqs. ( l37b and 
(l44l) fliat 



Taking Eqs. ( l44b and ( l37b into account, we can write Eqs. ( l43l l as 



(45) 



G 



cj-E(p) + iO aj + E(p)-iO' 
«'(P) vHp) 



oj-Eip) + iO w + £'(p)-/0 
As Aq 0, the gap A — > 0, but H and /c(p) remain finite if the spectrum becomes fiat, ^(p) = 0, and Eqs. ( l46b become 

1 1 



(46) 



F+(p,w) = -k(p) 

.,2 



0) + iO (jj - iO 

2, 



G(p,a>) 



«'(P) ^ v^(p) 



w -I- / ^ (D - iO 



(47) 



in the interval pi < p < pf. The parameters v(p) and u(p) are determined by the condition that the spectrum be flat: 
e(p) = iJ. If we take the Landau equation Q into account, this condition again reduces to Eqs. ( I^TT i and (l40l i for 
determining the minimum of the functional £'(n(p)). 

We construct the functions F^{p, o)) and G(p, <y) in the case where the constant /Iq is finite but small, such that v(p) 
and /(-(p) can be found on the basis of the FC solutions of Eq. (ISTT i. Then H, ^(p) and A are given by Eqs. ( l45T l. (l44l i. and 
dJTl l respectively. Substituting the functions constructed in this manner into (ij&b . we obtain F^(p, w) and G(p, oS) [QT]]- 
We note that Eqs. ( l44l i imply that the gap A is a linear function of Aq under the adopted conditions. As we shall see in 
Subsection l5.3l this gives rise to high-Tf at common values of the superconducting coupling constant 



5.3. The superconducting state at finite temperatures 

We assume that the region occupied by FC is small: {pf - pd/pp ^ 1 and Ai <sc Tf. Then, the order parameter /<-(p) 
is determined primarily by FC, i.e., the distribution function no(p) 11741 17511 . To be able to solve Eq. (l39l l analytically, we 
adopt the BCS approximation for the interaction [TgJ: /loV(p, Pi) - -Aq if |e(p) - /i| < cl)d and the interaction is zero 
outside this region, with ljd being a certain characteristic energy. As a result, the superconducting gap depends only on 
the temperature, A(p) = Ai(r), and Eq. (|39] l becomes 

Eo/2 

r 

1 = NfcIo ' ^ 



r ^ 



+ NlAo j , (48) 

£„/2 



where we introduced the notation ^ = e(p) - yu and the density of states Nfc in the interval (pf - pd or in the Eo-energy 
interval. It follows from Eq. (1411 1 that Nfc = (pf - Pf)Pf li^^^i)- Within the energy interval {(jJd - Eo/2), the density 
of states Nl has the standard form Nl - Myin. As £0 — > 0, Eq. (l48l l becomes the BCS equation. On the other hand, 
assuming that Eq < Iojd and discarding the second integral on the right-hand side of Eq. ( l48T l. we obtain 

AKO) = -^o^^(^/-^^) i„(i^V2) 

= 2/3gf ^^~^^ ln(l + V2), (49) 

where Sf - p^p/2M*^ is the Fermi energy and /3 - AoMyin is the dimensionless coupling constant. Using the standard 
value of p for ordinary superconductors, e.g., /3 ^ 0.3, and assuming that (pf - PfMpf - 0.2, we obtain a large value 
Ai(0) ~ O.lsf from Eq. (|49|i; for ordinary superconductors, this gap has a much smaller value: Ai(0) ~ lO^^e^. With the 
integral discarded earlier taken into account, we find that 

Ai(0) ^ ipsF ^^''^'' ln(l + V2) 
Pf ^ ' 

On the right-hand side of Eq. ( l50l l, the value of Ai is given by ( |49] l. As £0 ^ and Pf Pf, the first term on the 
right-hand side of Eq. ( l48T l is zero, and we obtain the ordinary BCS result with Ai oc exp(-l//lo). The correction related 
to the second integral in ( l48T l is small because the second term on the right-hand side of Eq. dSOt contains the additional 
factor /3. In what follows, we show that 2Tc - Ai(0). The isotopic effect is small in this case, because Tc depends on ojd 
logarithmical, but the eff'ect is restored as £0 — » 0. 

At r ^ Tc, Eqs. (gg and ^ are replaced by Eqs. ([SB and (O, which also hold for T^ < T ^ Tf. 



Kc - PF^^^^ " if r - r,; (5i) 
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Mp(. ^ Pf \J , Eo - AT, at r < T,.. (52) 



Pf-P 
AT 

Equation ( l48b is replaced by its standard generalization valid for finite temperatures 



1 ^f^Fch) , tanh — 



7 

+ NlAo tanh-!!— . (53) 

Because Ai(r — > T^) — > 0, Eq. ( l53T l implies a relation that closely resembles the BCS result jstl, 

2r, - Ai(0), (54) 

where Ai(r = 0) is found from Eq. ( l50t . Comparing ( |4TI ) and ( l42l i with ( ISTI ) and ( l52l i. we see that both M^^ and £0 are 
temperature-independent for r < Tc- 

5.4. Bogoliubov quasiparticles 

Equation ( [39l ) shows that the superconducting gap depends on the single-particle spectrum e(p). On the other hand, 
it follows from Eq. ( l37T i that e(p) depends on e(p) if Eq. ( l40t has a solution that determines the existence of EC as 
Aq — > 0. We assume that Aq is so small that the pairing interaction /loV(p, pi) leads only to a small perturbation of 
the order parameter /<'(p). Equation (l4Tl i implies that the effective mass and the density of states A^(0) oc M*^^ oc 1/Ai 
are finite. Thus, in contrast to the spectrum in the standard superconductivity theory, the single-particle spectrum e(p) 
depends strongly on the superconducting gap, and Eqs. (O and ( [39l ) must be solved by a self-consistent method. 

We assume that Eqs. O) and (l39b have been solved and the effective mass M^^ has been found. This means that 
we can find the quasiparticle dispersion law e(p) by choosing the effective mass M* equal to the obtained value of ML, 



and then solve Eq. ( l39b without taking (O into account, as is done in the standard BCS superconductivity theory 17811 
Hence, the superconducting state with EC is characterized by Bogoliubov quasiparticles |98] with dispersion ( |38] ) and 
the normalization condition v^(p) + M^(p) - 1 for the coefficients v(p) and m(p). Moreover, quasiparticle excitations of 
the superconducting state in the presence of EC coincide with the Bogoliubov quasiparticles characteristic of the BCS 
theory, and superconductivity with EC resembling the BCS superconductivity, which points to the applicability of the 
BCS formalism to the description of the high-T^ superconducting state l99ll . At the same time, the maximum value of the 
superconducting gap set by Eq. (l50b and other exotic properties are determined by the presence of EC. These results are 
in good agreement with the experimental facts obtained for the high-T^ superconductor Bi2Sr2Ca2Cu30io+tf 1100]. 

In constructing the superconducting state with EC, we returned to the foundations of the LEL theory, from which the 
high-energy degrees of freedom had been eliminated by the introduction of quasiparticles. The main difference between 
the LEL, which forms the basis for constructing the superconducting state, and the Eermi liquid with EC is that in the latter 
case we must increase the number of low-energy degrees of freedom by introducing the new type of quasiparticle with the 
effective mass M^^ and the characteristic energy £0 given by Eq. ( l42l i. Hence, the dispersion law e(p) is characterized 
by two types of quasiparticles with the effective masses M*^ and M^^^ and the scale Eq. The extended paradigm and 
new quasipa rticles determine the properties of the superconductor, including the lineshape of quasiparticle excitations 



I 74I75LIIOIII . while the dispersion of the Bogoliubov quasiparticles has the standard form. 

We note that for T < Tc, the effective mass M*^^ and the scale £0 are temperature-independent OlOlll . Eor T > Tc, 
the effective mass M^^^ and the scale Eq are given by Eqs. dSTT l and ( l33l l. Obviously, we cannot directly relate these new 
quasiparticles (excitations) of the Eermi liquid with EC to excitations (quasiparticles) of an ideal Eermi gas, as is done 
in the standard LEL theory, because the system is beyond ECQPT. The properties and dynamics of quasiparticles are 
given by the extended paradigm and closely related to the properties of the superconducting state and are of a collective 
nature, formed by ECQPT and determined by the macroscopic number of EC quasiparticles with momenta in the interval 
(pf- Pi). Such a system cannot be perturbed by scattering on impurities and lattice defects and, therefore, has the features 
of a quantum protectorate and demonstrates universal behavior [|74l I75[l84l85ll . 
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Several remarks concerning the quantum protectorate and the universal behavior of superconductors with FC are in 
order Similarly to the Landau Fermi liquid theory, the theory of high-T^ superconductivity based on FCQPT deals with 
quasiparticles that are elementary low-energy excitations. The theory provides a qualitative general description of the 
superconducting and the normal states of superconductors and HF metals. Of course, with phenomenological parameters 
(e.g., the pairing coupling constant) chosen, we can obtain a quantitative description of superconductivity, in the same way 
as this can be done in the Landau theory when describing a normal Fermi liquid, e.g., ''He. Hence, any theory capable of 
describing FC and compatible with the BCS theory gives the same qualitative picture of the superconducting and normal 
states as the picture based on FCQPT. Obviously, both approaches may be coordinated on the level of numerical results 
by choosing the appropriate parameters. For instance, because the formation of FC is possible in the Hubbard model 
118311 . it allows reproducing the results of the theory based on FCQPT. It is appropriate to note here that the correspo nding 



description restricted to the case of T = has been obtained within the framework of the Hubbard model Ill02lll03 l 



5.5. The pseudogap 

We now discuss some features of the superconducting state with FC 11571 l87l Il04ll considering two possible types 
of the superconducting gap A(p) determined by Eq. ( |39T l and the interaction /loV(p, pi). If the interaction is caused by 
attraction, occurring, for instance, as a result of an exchange of phonons or magnetic excitations, the solution of Eq. (39[ 
with an i-wave or 5 -i- li-mixed waves has the lowest energy. If the pairing interaction /toy(pi , pa) is a combi nation of a n 
attractive interaction and a strongly repulsive interaction, li-wave superconductivity may occur (e.g., see Refs ri05','l06'l). 
However, both the s- and tZ-wave symmetries lead to approximately the same result for the size of the gap Ai in Eq. (l50l l. 
Hence, d-wa\e superconductivity is not a universal and n ecessary property of high - Tc superconductors. This conclusion 



agrees with the experimental evidence described in Refs. 1 107 , 108 , 109 , 1 lOi 111 1 



We can define the critical temperature T* as the temperature at which Ai(r*) s 0. For T > T*, Eq. (153b has 
only the trivial solution Ai = 0. On the other hand, the critical temperature can be defined as the temperature at 
which superconductivity disappears and the gap occupies only a part of the Fermi surface. Thus, there are two different 
temperatures and T*, which may not coincide in the case of the li-wave symmetry of the gap. As shown in Refs [57, StI], 
in the presence of FC , Eq. (|53] | has nontrivial solutions atT^ < T < T*, when the pairing interaction /loy(pi , pa) consists 
of attraction and strong repulsion, which leads to d-wave superconductivity. In this case, the gap A(p) as a function of the 
angle (/>, or A(p) = A(pf, (p), has new nodes at J > T„ode, as shown in Fig. 17) 187.1 . 
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Figure 7: The gap A{pf,(p) as a function of (p calculated for three values of the temperature expressed in units of T„„jg ^ T^. The solid curve (a) 
represents the function A{pf,<l>) calculated for the temperature 0.9r„orf(,. The dashed curves (b) represents the same function at T = T„oji,, and the 
dotted curve (c) depicts the function calculated at 1.2 T„ocie- The an'ows indicate the region 6^ limited by the two new zeros that emerge at T > T„ode- 



Figure |2] shows the ratio A{pf,(p)IT* calculated for three temperatures: 0.9 Tnode, Tnode and 1.2 T„ode- In contrast to 
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curve (a), curves (b) and (c) have flat sections. Clearly, the flattening occurs because of the two new zeros that emerge 
at r = Tnode- As the temperature increases, the region 6c between the zeros (indicate d by arrow s in Fig. |7]) increases in 
size. It is also clear that the gap A is very small within the interval 9c. It was found in 111 12l ll that the magnetism and 
the superconductivity affect each other, which leads to suppression of the magnetism at temperatures below Tc- In view 
of this, we can expect suppression of superconductivity due to magnetism. 

Thus, we may conclude that the gap in the vicinity of Tc can be destroyed by strong antiferromagnetic correlations 
(or spin density waves), impurities, and sizable inhomogeneities existing in high-Tc superconductors 111 1411 . Because the 
superconducting gap is destroyed in a macroscopic region of the phase space, 9c, superconductivity is also destroyed, and 
therefore Tc - Tnode- The exact value of Tc is determined by the competition between the antiferromagnetic state (or 
spin density waves) and the superconductivity in the interval 9c. The behavior and the shape of the pseudogap closely 
resemble the similar characteristics of the superconducting gap, as Fig. |2]shows. The main difference is that the pseudogap 
disappears in the segment 9c of the Fermi surface, while the gap disappears at isolated nodes of the d-wave. This result 
is in accord with observations 111 1511 . Our estimates show that for small values of the angle i/r, the function 0t (i/') rapidly 
increases, 0c(i/') - y/tp. These estimates agree with the results of numerical calculations of the function 9c{[T - Tc]/Tc), 
(Fig. |8]l. Hence, we may conclude that Tc is close to T„ode- Thus, the pseudogap state appears at T > Tc - T„ode and 
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Figure 8: . The result of a numerical calculation of the angle 8^ separating two zeros as a function of {T - Tc)/Tc. 



disappears at temperatures T > T* at which Eq. (I53I I has only the trivial solution Ai s 0. Obviously, Ai determines T* 
and not Tc, with the result that Eq. ( |54] | should be rewritten as 

IT* ^ Ai(0). (55) 

The temperature T* has the physical meaning of the temperature of the BCS transition between the state with the order 
parameter k i^O and the normal state. 

At temperatures below T < Tc, the quasiparticle excitations of the superconducting state are characterized by the pres- 
ence of sharp peaks. When the temperature becomes high (T > Tc) and A(9) = in the interval 9c, normal quasiparticle 
excitations with a width y appear in the segments 9c of the Fermi surface. A pseudogap exists outside the segments 9c, 
and the Fermi surface is occupied by excitations of the BCS type in this region. Excitations of both types have widths 
of the same order of magnitude, transferring their energy and momenta into excitations of normal quasiparticles. These 
results are in accord wi th strong indications of the pa iring or the formation of preformed pairs in the pseudogap regime at 
temperatures above Tc ful [Tm [TH [TTl [TTl [l20ll . 



We now estimate the value of y. If the entire Fermi surface were occupied by the normal state, the width y would be 
y ~ N{QifT^ I Sd{Tf wifli the density of states A^(0) ~ M*(T) ~ l/T [see Eq. The dielecti-ic constant ed(T) ~ A^(0) 
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and, hence, y ~ T 11511 15211 . However, only a part of the Fermi surface within 6c is occupied by normal excitations in 
our case. Therefore, the number of states accessible for quasiparticles and quasiholes is proportional to 9c, and the factor 

is replaced by the factor T^O^. Taking all this into account yields y ~ O^T ~ T{T - Tc)/Tc ~ {T - Tc). Here, we 
ignored the small contribution provided by excitations of the BCS type. It is precisely for this reason that the width y 
vanishes at T = Tc- Moreover, the resistivity of the normal state p(T) ccy cc (T - Tc), because y ~ (T - Tc). Obviously, at 
temperatures T > T* , the relation p{T) cc y cc T remains valid up to T ~ Tf, and Tj- may be as high as the Fermi energy 
if FC occupies a significant part of the Fermi volume. 

The temperature T„ode is determined mainly by the repulsive interaction, which is part of the pairing interaction 
'ioV(Pi,P2)- The value of the repulsive interaction, in turn, may be determined by the properties of the materials, such 
as composition or doping. Because superconductivity is destroyed a t Tc - T„ode, the ratio 2Ai/Tc may vary within broad 
limits and strongly depends on the properties of the material 1571 l87L 1 10411 . For instance, in the case of Bi2Sr2CaCu2 Oe+a 



it is assumed that superconductivity and the pseudogap are of common origin: 2Ai/Tc — 28, while lAjJT* ^ 4, which 



agrees with the experimental data obtained in measurements involving other high-T^ superconductors 1I105I1 . 

We note that Eq. ( |55] ) also provides a good description of the maximum value of the gap A] in the case of d- wave 
superconductivity, because different regions with the maximum density of states may be considered unrelated flOel]. We 
may also conclude that without a strong repulsion, with which i-wave pairing is possible, there can be no pseudogap. 
Thus, the transition from the superconducting gap to the pseudogap may proceed only in the case of t/-wave pairing, when 
superconductivity is destroyed at Tc - T„ode and the superconducting gap gradually transforms into a pseudogap, which 

The fact that there is no pseudogap in the case of i-wave pairing 



closes at a certain temperature T* > Tc [57, 87, 104l 



agrees with the experimental data (e.g., see Ref. Ill 1 111 ') 



5.6. Dependence of the critical temperature Tc of the superconducting phase transition on doping 

We examine the maximum value of the superconducting gap Ai as a function of the number density x of mobile charge 
carriers, which is proportional to the degree of doping. Using Eq. 041 ). we can rewrite Eq. (l49b as 

A,^^Oc^c^_ (56) 
sf xfc 

where we took into account that the Fermi level sf and that the number density x ~ pj^/(2M*), with the result 

that Sf X. It is realistic to assume that Tc cc Ai, because the curve Tdx) obtained in experiments with high-Tc 
super conductors must be a smooth function of x. Hence, we can approximate Tc(x) by a smooth bell-shape function 



11211: 

Tc(x) cc j3(xFc - x)x. (57) 

To illustrate the application of the above analysis, we examine the main features of a superconductor that can hy- 
pothetically exist at room temperature. Such a superconductor must be a two-dimensional structure, just as high-Tc 
superconducting cuprates are. Equation (|49] | implies that Ai ~ /3sf Bearing in mind that FCQPT occurs at 

rs ~ 20 in 3D systems and at r, ~ 8 in 2D systems |70], we can expect that in 3D systems Ai amounts to 10% of the 
maximum size o f the superconducting gap in 2D systems, which in our case amounts to 60 mV for lightly doped cuprates 
with Tc = 70 K fl23\. On the other hand, Eq. i4% implies that Ai may be even larger. A] ~ 75 mV. We can expect that 
Tc ~ 300 K in the case of i-wave pairing, as the simple relation 2Tc - Ai implies. Indeed, we can take sf ~ 500 mV, 
y6 ~ 0.3, and (pf- pd/pF ~ 0.5. 

Thus, the hypothetical superconductor at room temperature must be an s-wave superconductor in order to eliminate 
the pseudogap effect, which dramatically decreases the temperature Tc at which superconductivity is destroyed. We note 
that the number density x of mobile charge carriers must satisfy the condition x < xfc and must be varied to reach the 
optimum degree of doping x„p, ^ xfcI'^. 

5.7. The gap and heat capacity near Tc 

We now calculate the gap and heat capacity at temperatures T Tc. Our analysis is valid if T* ^ Tc, since otherwise 
the discontinuities in the heat capacity considered below are smeared over the temperature interval between T* and r^.. 
To simplify matters, we calculate the leading contribution to the gap and heat capacity related to FC. We use Eq. (l53l l to 
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find the function Ai(r — > Tc) simply by expanding the first integral on its righ t-han d side in powers of Ai and dropping 
the contribution from the second integral. This procedure leads to the equation OlOlll 



Ai(T) ^ 3 AT, 



T 



(58) 



Therefore, the gap in the spectrum of single-particle excitations behaves in the ordinary manner. 
To calculate the heat capacity, we can use the standard expression for the entropy S 117 811 : 



S(T) 



J [/(p)ln/(p) + (l 
dp 



/(P)) 



X ln(l-/(p))] 



(271) 



2' 



(59) 



where 



The heat capacity C is given by 



/(P) -- 
£(P) -- 

C(T) = T 



(l+exp[£(p)/r])- 



■ju)2 + A?(r). 



(60) 



J/(£)(l -/(£)) 



t/Ai(r)] 



+ TAi(T)- 



dT 



\d^ 



|| J f(E)(l-f(E)) 

dAiiry 



E^ + TAi(T)- 



dT 



d^. 



(61) 



In deriving Eq. (|6TI ). we again used the variable ^, the above notation for the density of states, Nfc and Nl, and the 
notation E - + A^{T). Equation ( ICTb describes a jump in heat capacity, 6C{T) - Cs{T) - C„{T), where CsiT) and 
C„{T) are respectively the heat capacities of the superconducting and normal states at T,; the jump is determined by the 
last two terms in the square brackets on the right-hand side of this equation. Using Eq. ( l58b to calculate the first term on 
the right-hand side of Eq. dMT l, we find lIlOlll 



(62) 



where n = 1 in the 2D case and n - 2 m the 3D case. This result differs from the ordinary BCS result, according to 
which the discontinuity in the heat capacity is a linear function of T,- The jump 6C{Te) is independent of because, 
as Eq. ( |52] | shows, the density of state varies in inverse proportion to T,. We note that in deriving Eq. ( l62b we took the 
leading contribution coming from EC into account. This contribution disappears as Eq — > 0, and the second integral on 
the right-hand side of Eq. dFTT l yields the standard result. 

As we will show in Section [10] [see Eq. ( I134b l. the heat capacity of a system with EC behaves as C„(T) cc ^jT/Tf. 
The jump in the heat capacity given by Eq. (l62b is temperature-independent. As a result, we find that 




(63) 



C„{Tc) \ Tc Pf 

In contrast to the case of normal superconductors, in which 6C(T,)/C„(T,) = 1.43 jEOl, in our case Eq. ( l63b implies that 
the ratio 6C(Tc)/Cn(Tc) is not constant and may be very large when Tf/T, » 1 1791 101 1. It is instructive to apply this 
analysis to CeCoIns, where -2.3 K |79]. In this material lll23ll . SC/C,, — 4.5 is substantially higher than the BCS 
value, in agreement with Eq. (|63]) . 
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6. The dispersion law and lineshape of single-particle excitations 



The recently discovered break in the dispersion of quasiparticles at energies between 40 and 70 mV, resulting in a 
change in the quasiparticle speed at this energy 19111 . can hardly be explained by the marginal Fermi-liquid 



theory, because this theory contains no additional energy scales or parameters that would allow taking the break into 
account 1.76., 77i1 . We could assume that the break, which leads to a new energy scale, occurs because of the interaction of 
electrons and collective excitations, but then we would have to discard the idea of a quantum protectorate, which would 
contradict the experimental data llsii 851. 



As shown in Sections|4]and|5] a system with FC has two effective masses: M^^, which determines the single-particle 
spectrum at low energies, and M^, which determines the spectrum at high energies. The fact that there are two effective 
masses manifests itself in the form of a break in the quasiparticle dispersion law. The dispersion law can be approximated 
by two straight lines intersecting at a binding energy £'0/2 [see Eqs. (l33T l and (l42T i1. The break in the dispersion law occurs 
at temperatures much lower than T <k Tj, when the system is in the superconducting or normal state. Such behavior is in 
good agreement with the experimental data fs§\. It is pertinent to note that at temperatures below T < Tc, the effective 
mass is independent of the momenta pf, Pf, and p,, as shown by Eqs. Wli and (|49] l: 

~ — . (64) 
This formula implies that M^^ is only weakly depende nt on jc if a dependence of /Iq on x is allowed. This result is in good 



Y depenae r 

agreement with the experimental facts Hl24i 125ill26ll . The same is true of the dependence of the Fermi velocity 
PfIM*p^ o n X because the Fermi momentum pf ~ Vn is weakly dependent on the electron number density n - «o(l - x) 



1 124112511 : here, «o is the single-particle electron number density at half-filling. 



Because is the coupling constant that determines the magnitude of the pairing interaction, e.g., the electron-phonon 
interaction, we can expect the break in the quasiparticle dispersion law to be caused by the electron-phonon interaction. 
The phonon scenario could explain the constancy of the break at T > Tc because phonons are temperature independent. 
On the other hand, it was found that the quasiparticle dispersion law distorted by the interaction with phonons has a 
tendency to restore itself to the ordinary single particle dispersion law when the quasiparticle energy becomes higher than 



the phonon energy 1112711 . However, there is no experimental evidence that such restoration of the dispersion law actually 
takes place llssll . 

The quasiparticle excitation curve Uq, a)) is a function of two variables. Measurements at a constant energy lj = wq, 
where (dq is the single particle excitation energy, determine the curve L{q, co - wq) as a function of the momentum q. 
We established above that M^^, is finite and constant at temperatures not exceeding Tc- Hence, at excitation energies 
u) < En, the system behaves as an ordinary superconducting Fermi liquid with the effective mass determined by Eq. iAli 
HIzlliTtl- At r, < r, the effective mass Mpi^ is also finite and is given by Eq. Olb . In other words, at w < E^), 
the system behaves as a Fermi liquid whose single-particle spectrum is well defined and the width of the single-particle 
excitations is of the order of T fl^Jl^Jsl^. Such behavior has been observed in experiments in measuring the quasiparticle 
excitation curve at a fixed energy [34, 90, 128]. 

The quasiparticle excitation curve can also be described as a function of a>, at a constant momentum q = ^o- For small 
values of oj, the behavior of this function is similar to that described above, with L{q - qo, a)) having a characteristic 
maximum and width. For o) > Eq, the contribution provided by quasiparticles of mass becomes significant and leads 
to an increase in the function L{q - qo, oS). Thus, L{q - qo, a>) has a certain structure of maxima and minima il29ll 



directly determined by the existence of two effective masses, M^^ and M* 1174 175L 18711 . We conclude that, in contrast to 
Landau quasiparticles, these quasiparticles have a more complicated spectral lineshape. 

We use the Kramers-Kronig transformation to calculate the imaginary part Im l)(p, g ) of the self-energy part 2(p, e). 
But we begin with the real part Re2;(p, e), which determines the effective mass M* lll30ll . 

\m Pf dp )j \ ds ) 



M* 

The corresponding momenta p and energies e satisfy the inequalities \p - PfMpf ^ 1, and e/e/r <K 1. We take ReE(p,e) 
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in the simplest form possible that ensures the variation of the effective mass at the energy £0/2, 



M* 

ReS(p, E) = -e— ^ + 



FC 



■M* 



m \ Z / m 



(66) 



where 6(s) is the step function. To ensure a smooth transition from the single-particle spectrum characterized by M^^ to 
the spectrum characterized by M*, we must replace the step function with a smoother function. Substituting Eq. ( |66] ) in 
Eq. ( |65] ). we see that M* ^ M^^-, within the interval (-£'0/2, Eo/2), while M* ^ M*^ outside this inte rval. Applying the 
Kramers-Kronig transformation to ReZ(p, s), we express the imaginary part of the self-energy as OlOlll 



ImS(p, e) ~ e 



2 Kc 
EfM 



■Ml 



filn 



2s + £■() 



2s - Eq 



Eo, 



4s^ 



-^0 



^0 



(67) 



Clearly, with s/Eq <k 1, the imaginary part is proportional to e^; at 2s/E() ^ 1, we have ImS ~ s, and for Eq/s «: 1, the 
main contribution to the imaginary part is approximately constant. 

It follows from Eq. (l67T i that as £0 — > 0, the second term on its right-hand side tends to zero and the single-particle 
excitations become well-defined, which resembles the situation with a normal Fermi liquid, while the pattern ofminima 
and maxima eventually disappears. Now the quasiparticle renormalization factor zip) is given by the equation il30ll 



1 

i(p) 



1 



gRe£(p,g) 



(68) 



Consequently, from Eqs. ( |67] ) and ( |68] ) for T < Tc, the amplitude of a quasiparticle on the Fermi surface increases as 
the characteristic energy £0 decreases. Equations (l42l i and (ISTT i imply that £0 ~ (xpc - x)/xfc- When T > Tc, it follows 
from (l66T l and ( |68] | that the quasiparticle amplitude increases as the effective mass M^^ decreases. So, from Eqs. dSTT i and 
(l34ll M^.^ ~ ipf- Pi)/ Pf ~ {xfc — x)Ixfc- As a result, we conclude that the amplitude increases with the doping level and 
the single-particle excitations are better defined in heavily doped samples. As x — > Xfc, the characteristic energy £0 — » 
and the quasiparticles become normal excitations of LFL. We note that such behavior has been observed in e xperiments 
with heavily doped Bi2212, which demonstrates high-r^ superconductivity with a gap of about 10 mV ill 3 ill . The size 
of the gap suggests that the region occupied by FC is small because £o/2 ^ A]. For x > Xfc and low temperatures, the 
HE liquid behaves as LFL (see Fig. |6|and Section |9]l. Experimental data show that, as expected, the LFL state exists in 
super-heavily doped nonsuperconductingLai.7Sro.3Cu04 lll32ill33ll . 



7. Electron liquid with FC in magnetic fields 

In this Section, we discuss the behavior of HE liquid with FC in magnetic field. We assume that the coupling constant 
is nonzero, Aq + 0, but is infinitely small. We found in Section |5] that at T = the superconducting order parameter 
/(■(p) is finite in the region occupied by FC and that the maximum value of the superconducting gap Ai oc is infinitely 
small. Hence, any weak magnetic field B is critical and destroys A:(p) and FC. Simple energy arguments suffice to 
determine the type of rearrangement of the FC state. On the one hand, because the FC state is destroyed, the gain in energy 
A£b 0^ tends to zero as B — > 0. On the other hand, the function no(p), which occupies the finite interval {pj - pi) in 
the momentum space and is specified by Eq. (1211 1 or (l42l l. leads to a finite gain in the ground-state energy compared to the 
ground-state energy of a normal Fermi liquid 141] . Thus, the distribution function is to be reconstructed so that the order 
parameter is to vanish while a new distribution function is to deliver the same ground state energy. 

7.1. Phase diagram of electron liquid in magnetic field 

Thus, in weak magnetic fields, the new ground state without FC must have almost the same energy as the state with 
FC. As shown in Section [15] such a state is formed by multiply connected Fermi spheres resembhng an onion, in which 



26 



a s mooth distribution function of quasiparticles, no(p), is replaced in the interval {pf - pd with the distribution function 
I9i[l3l 



V(P) = ^ - P2k-l)9{P2k - P)- 



(69) 



where the parameters pi < pi < pi < ■ ■ ■ < Pin ^ Pf ^"e chosen such that they satisfy the normalization condition and the 
condition needed for the conservation of the number of particles: 



Figure |9] shows the corresponding multiply connected distribution. For definiteness, we present the most interesting case 
of a three-dimensional system. The two-dimensional case can be examined similarly. We note that the possibility of the 
existence of multiply connected Fermi spheres was studied in e.g. 1.23, . 135., 136.. 137.1 . 



V(p) 







p. Pa 



2n-1 
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Figure 9: The function y(p) for the multiply connected distribution that replaces the function «o(p) in the region (pf - pi) occupied by FC. The momenta 
satisfy the inequalities p, < pf < pf. The outer Fermi surface at p ^ p2„ - Pf has the shape of a Fermi step, and therefore the system behaves like LFL 
at sufficiently low temperatures. 



We assume that the thickness of each inner slice of the Fermi sphere, dp - p2k+i - Pik, is determined by the magnetic 
field B. Using the well-known rule for estimating errors in calculating integrals, we find that the minimum loss of the 
ground-state energy due to slice formation is approximately (Sp)'*. This becomes especially clear if we account for 
the fact that the continuous FC functions «o(p) ensure the minimum value of the energy functional £'[n(p)], while the 
approximation of v(p) by steps of width dp leads to a minimal error of the order of (dp)*. Recalling that the gain due to 
the magnetic field is proportional to and equating the two contributions, we obtain 



6p cc Vb. 



(70) 



Therefore, as T — > 0, with B — > 0, the slice thickness 6p also tends to zero and the behavior of a Fermi liquid with FC 
is replaced with that of LFL with the Fermi momentum pf. Equation ( l40l i implies that pj > pf and the electron number 
density x remains constant, with the Fermi momentum of the multiply connected Fermi sphere p2n - Pf > Pf (see Fig. 
|9]l- We see in what follows that these observations play an important role in studying the behavior of the Hall coefficients 
Rh{B) as a function of B in heavy-fermion metals at low temperatures. 

To calculate the effective mass M*{B) as a function of the applied magnetic field B, we first note that at T = the 
field B splits the FC state into Landau levels, suppresses the superconducting order parameter /c(p), and destroys FC, 
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which leads to restoration of LFL 138 1. The Landau levels near the Fermi surface can be approximated by separate 



slices whose thickness in momentum space is 6p. Approximating the quasiparticle dispersion law within a single slice, 
e{p) - jj. ~ ip - Pf + Sp)(p - pf)/M*, we find the effective mass M*(B) ~ M* /(6p/pf). The energy increment AEfc 
caused by the transformation of the PC state can be estimated based on the Landau formula lEoll 

AEfc^ [(Mp)-f^)Sn{p)-^. (71) 

The region occupied by the variation (Jn(p) has the thickness 6p, with (e(p) - p) ~ {p - pf)pf/M*(B) ~ 6ppf/M*(B). As 
a result, we find that AEfc ~ p^jdp^ /M*{B). On the other hand, there is the addition AEb ~ iBfiB)^M*(B)pf caused by 
the applied magnetic field, which decreases the energy and is related to the Zeeman splitting. Equating AEb and AEfc 
and recalling that M*(B) oc 1 fSp in this case, we obtain the chain of relations 

6p^ 1 



M'{E) (M*(B))^ 
which implies that the effective mass M*{B) diverges as 



oc B'-M^B). (72) 



M'iB) oc ^ (73) 
where Bcq is the critical magnetic field, which places HF metal at the magnetic-field-tuned quantum critical point and 



nullifies the respective Neel temperature, Tnl{Bm) - 114911 . In our simple model of HF liquid, the quantity Bco is a 



parameter determined by the properties of the specific metal with heavy fermions. We note that in some cases B^q - 0, 
e.g., the HF metal CeRuiSia has no magnetic order, exhibits no superconductivity, and does not behave like a Landau 
Fermi liquid even at the lowest reached temperatures ll92ll . 

Formula (l73T l and Fig. |9]show that the application of a magnetic field B > B^ brings the FC system back to the LFL 
state with the effective mass M*(B) that depends on the magnetic field. This means that the following characteristic of 
LFL are restored: C/T - joiB) oc M*(B) for the heat capacity andxoiB) oc M*(B) for the magnetic susceptibility. The 
coefficient A(B) determines the temperature-dependent part of the resistivity, p{T) - po + Ap, where po is the residual 
resistivity and Ap - A{B)T^. Because this coefficient is directly determined by the effective mass, A{B) oc {M*{B)f l,139ll . 
Eq. (l73T l yields 

A(B)cx-i— . (74) 

Thus, the empirical Kadowaki- Woods relation ll27ll K — A/yf^ ^ const is valid in our case lll39|] . Furthermore, K may 
depend on the degree of degeneracy of the quasiparticles. With this degenera tion, the K adowaki- Woods relation provides 
a good description of the experimental data for a broad class of HF metals 1 14d 141 1. In the simplest case, where HF 



liquid is formed by_^spin-l/2 quasiparticles with the degeneracy degree 2, the value of K turns out to be close to the 
empirical value 1 139] known as the Kadowaki- Woods ratio |27]. Hence, under a magnetic field, the system returns to the 
state of LFL and the constancy of the Kadowaki-Woods relation holds. 

At finite temperatures, the system remains in the LFL state, but when T > T*{B), the NFL behavior is restored. As 
regards finding the function T*(B), we note that the effective mass M* characterizing the single-particle spectrum cannot 
change at T*{B) because no phase transition occurs at this temperature. To calculate M*{T), we equate the effective mass 
M*(r) in Eq. (O to M*(B) in M*(T) ~ M*(B), 

^ T'iB) 0^ tTTT;- 0^ - Bco, (75) 



M*{T) M*(B) 

whence 

r{B) oc ^B - B,„. (76) 

At temperatures T > T*(B), the system returns to the NFL behavior and the effective mass M* specified by Eq. ( ISTT l. 
Thus, expression ( |76] | determines the line in the T - B phase diagram that separates the region where the effective mass 
depends on B and the heavy Fermi liquid behaves like a Landau Fermi liquid from the region where the effective mass 
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Figure 10: Schematic T -B phase diagram of heavy electron liquid. B^o denotes the magnetic field at which the eft'ective mass diverges as given by j73t . 
The horizontal an'ow illustrates the system moving in the NFL-LFL direction along B at fixed temperature. As shown by the dashed curve, at B < 5^0 
the system can be in its ferromagnetic (FM), antiferromagnetic (AFM) or superconducting (SC) states. The NFL state is characterized by the entropy 
So given by Eq. j27t . The solid curve T'(B) separates the NFL state and the weakly polarized LFL one and represents the transition regime. 



is temperature-dependent. At T*{B), the temperature dependence of the resistivity ceases to be quadratic and becomes 
Hnear 

A schematic T - B phase diagram of HF Hquid with FC in magnetic field is shown in Fig. [10] At magnetic field 
B < Bc{) the FC state can captured by ferromagnetic (FM), antiferromagnetic (AFM) and superconducting (SC) states 
lifting the degeneracy of the FC state. It follows from ( |76] | that at a certain temperature T*{B) <sc T j, the heavy-electron 
liquid transits from its NFL state to LFL one acquiring the properties of LFL at {B - Bco) oc {T*{Byf'. At temperatures 
below T*(B), as shown by the horizontal arrow in Fig. [TO] the heavy-electron liquid demonstrates an increasingly metallic 
behavior as the magnetic field B increases, because the effective mass decreases [see Eq. (l73Tl1. Such behavior of the 
effective mass can be observed, for instance, in measurements of the heat capacity, magnetic susceptibility, resistivity, 
and Shubnikov-de Haas oscillations. From the T - B phase diagram shown in Fig. [10] and constructed in this manner, 
it follows that a unique possibility emerges where a magnetic field can be used to control the variations in the physical 
nature and type of behavior of the electron liquid with FC. 

We briefly discuss the case where the system is extremely close to FCQPT on the ordered size of this transition, and 
hence dppc - (Pf - pd/PF 1- Because 6p oc M*{B), it follows from Eqs. (TTOt and (TTSl l that 



5p B-B,o 

r , (77) 

Pf V Bm 

where Qc is a constant of the order of unity, ~ 1. As the magnetic field B increases, Sp/pF becomes comparable to 
Spfc, and the distribution function v(p) disappears, being absorbed by the ordinary Zeeman splitting. As a result, we 
are dealing with HF liquid located on the disordered side of FCQPT. We show in Section [TO] that the behavior of such a 
system differs markedly from that of a system with FC. Equation (fTTT i impUes that the relatively weak magnetic field Bcr, 

Bred = ^-^^ = (6pFcf ~ B,r, (78) 
BcO 

where Bred is the reduced field, takes the system from the ordered side of the phase transition to the disordered if Spfc <K 1 ■ 

7.2. Dependence of effective mass on magnetic fields in HF metals and high-Tc superconductors 

Observations have shown that in the normal state obtained by applying a magnetic field whose strength is higher 
than the maximum critical field Bc2 that destroys superconductivity, the heavily doped cuprate (Tl2Ba2Cu06+a) 12 61 and 
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the optimally doped cuprate (Bi2Sr2Cu06+5) 112911 exhibit no significant violations of the Wiedemann-Franz law. Studies 
of the electron-doped superconductor Pro.ciiLaCe() ()9Cu04_y {Tc-2A K), revealed that when a magnetic field destroyed 
superconductivity in this material, the spin-lattice relaxation constant l/Ti obeyed the relation T\T - const, known as 
the Korringa law, down to temperatures about T ^ 0.2 K [142, 143]. At higher temperatures and in magnetic fields 
up to 15.3 T perpendicular to the Cu02 plane, the ratio 1/TiT remains constant as a function of T f or T < 55 K. In 
the temperature range from 50 to 300 K, the ratio 1/TiT decreases as the temperature increases Measurements 
involving the heavily doped nonsuperconducting material Lai 7Sro.3Cu04 have shown that the resistivity p varies with T 
as and that the Wiedemann-Franz law holds il32ill33ll . 

Because the Korringa and Wiedemann-Franz laws strongly indicate the presence of the LFL state, experiments show 
that the observed elementary excitations cannot be distinguished from Landau quasiparticles in high-Tc superconductors. 
This places seve re restricti ons on models describing hole- or electron-doped high-Tc superconductors. For instance, for a 
Luttinger liquid i l44 , 145 1. for spin-charge separation |8], and in some t- J models 1 146], a violation of the Wiedemann- 



Franz law was predicted, which contradicts experimental evidence and points to the limited applicability of these models. 

If the constant Aq is finite, then a HF liquid with FC is in the superconducting state. We examine the behavior of 
the system in magnetic fields B > Bc2- In this case, the system becomes LFL induced by the magnetic field, and the 
elementary excitations become quasiparticles that cannot be distinguished from Landau quasiparticles, with the effective 
mass M*(B) given by Eq. ( |73] ). As a result, the Wiedemann-Franz law holds as T — > 0, which agrees with the experimental 
facts 1 2^ 29 1 . The low-temperature properties of the system depend on the effective mass; in particular, the resistivity 
p(T) behaves as given by Eq. (fTSl ), with A{B) oc (M*(B))^. Assuming that the critical field B = Bco in the case of high-Tc 
superconductors, we deduce from Eq. ( l73l l that 

70 - BcO = const. (79) 

Taking Eqs. ( |74] | and ( |79] | into account, we find that 

ro ~ A{B) VB - B,o. (80) 

At finite temperatures, the system remains LFL, but for T > T*{B) the effective mass becomes temperature-dependent, 
M* oc 1/r, and the resistivity becomes a linear function of the temperature, p(T) oc T |147ll . Such behavior of the 
resistivity has been observed in the high-r^ superconductor Tl2Ba2Cu06+5 (Tc < 15 K) 114811 . At < 10 T, the resistivity 
is a linear function of the temperature in the range from 120 mK to 1.2 K, and at B = 10 T the t emperature dependence 



of the resistivity can be written in the form p(T) oc AT^ in the same temperature range Ill48lll49[l . clearly demonstrating 
that the LFL state is restored under the application of magnetic fields. 

In LFL, the spin-lattice relaxation parameter 1/Ti is determined by the quasi particles near the Fermi level, whose 



population is proportional to M*T, whence l/TiT cc M*, and is a constant quantity II 1421114311 . When the superconducting 



ipa 



state disappears as a magnetic field is applied, the ground state can be regarded as a field-induced LFL with the field- 
dependent effective mass. As a result, T\T - const, which implies that the Korringa law holds. According to Eq. (l73T l. 
the ratio 1/TiT oc M*{B) decreases as the magnetic field increases at T < T*(B), whereas in the case of a Landau Fermi 
liquid it remains constant, as noted above. On the other hand, at T > T*(B), the ratio l/TiT is a decreasing function of 



the temperature, 1 /TiT oc M*{T). These results are in good agreement with the experimental facts B 14311 . Because T*{B) 
is an increasing function of the magnetic field [see Eq. (f76]ll, the Korringa law remains valid even at higher temperatures 
and in stronger magnetic fields. Hence, at Tq < T*{Bq) and high magnetic fields B > Bo, the system demonstrates distinct 
metallic behavior, because the effective mass decreases as B increases, see Eq. ( l73T l. 

The existence of FCQPT can also be verified in experiments, because at number densities x > Xfc or beyond the 
FCQPT point, the system must become LFL at sufficiently low temperatures 1 1381 . Experiments have shown that such 



a liquid indeed exists in the heavily doped non-superconducting compound Lai.7Sro.3Cu04 01321 Il33ll . It is remarkable 



that fo r r < 55 K, the resistivity exhibits a T^-behavior without an additional linear term and the Wiedemann-Franz law 



holds Ill32lll33|]. At t emperatures above 55 K, experimenters have observed significant deviations from the LFL behavior 
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Observations |0, 134ll are in accord with these experimental findings showing that the system can again be returned to the 



LFL state by applying sufficiently strong magnetic fields (also see Section|9]l. 

7.2.7. Common QCP in the high-Tc Tl2Ba2Cu06+x ond the HF metal YbRh2Si2 

Under the application of magnetic fields B > Bc2 > BeO and at T < T*{B), a high-Tc superconductor or HF metal can 
be driven to the LFL state with its resistivity given by Eq. ( fTSl l. In that case measurements of the coefficient A produce 
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information on its dependence on the applied field. We note that relationships between critical magnetic fields Bc2 and 
BcQ are clarified in Subsection |9.9l 




B[T] 



B[T] 



Figur e 11: The charge tansport coefficient A(B) as a function of magnetic field B obtained in measurements on YbRh2Si2 fl?] and Tl2Ba2Cu06+x 
llSOll . The different field scales are clearly seen. The solid curves represent our fit by using Eq. )8U 

Precise measurements of the coeflicient A{B) on high-Tc Tl2Ba2Cu06+x lll50ll allow us to establish relationships 
between the physics of both high-Tc superconductors and HF metals and clarify the role of the extended quasiparticle 
paradigm. The A(B) coefficient, being proportional to the quasiparticle-quasiparticle scattering cross-section, is found to 
be A oc (M*(B))2 dliH. With respect to Eq. this implies that 



A{B) - Ao + 



D 



B-B, 



(81) 



cO 



where Aq and D are fitting parameters. 

Figure [TT| reports the fit of our theoretical dependence dSTT l to the experimental data for the measurements of the 
coefficient A(B) for two different classes of substances: HF metal YbRh2Si2 (with Bco - 0.06 T, left panel) (l^ and high- 
Tc Tl2Ba2Cu06+x (with Bco = 5.8 T, right panel) [ISO). In Fig. (TT] left panel, A{B) is shown as a function of magnetic 
field B, applied boffi along and perpendicular to the c axis. For ffie latter the B values have been multiplied by a factor of 
1 1 lIlSll . The different scales of field Bco are clearly seen and demonstrate that Bco has to be taken as an input parameter. 
Indeed, the critical field of Tl2Ba2Cu06+x with Bco = 5.8 T is 2 orders of magnitude larger than that of YbRh2Si2 with 
Bm = 0.06 T. 

Figure [TT] displays good coincidence of the theoretical dependence ( l74l i with the experimental facts 1 150l 151 1. This 
means that the physics underlying the field-induced reentrance of LFL behavior, is the same for both classes of substances. 
To further corroborate this point, we rewrite Eq. ( ISTI ) in ffie reduced variables A/Aq and B / Bco- Such rewriting immediately 
reveals the scaling nature of the behavior of these two substances - both of them are driven to common QCP related to 
FCQPT and induced by the application of magnetic field. As a result, Eq. dSTT i takes the form 



MB) 



^ 1 -I- 



B/B,o - 1 ' 



(82) 



where = D/(AoBco) is a constant. From Eq. ( l82l i it is seen that upon applying the scaling to both coefficients A{B) 
for Tl2Ba2Cu06+x and A(B) for YbRh2Si2 they are reduced to a function depending on the single variable B/Bco thus 
demonstrating universal behavior To support Eq. ( l82l i. we replot both dependencies in reduced variables A/Ao and B/BcO 
in Fig. [121 Such replotting immediately reveals the universal scaling nature of the behavior of these two substances. It is 
seen from Fig. [T2]that close to ffie magnetic induced QCP ffiere are no "external" physical scales revealing the scaling. 
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Figure 12: Normalized coefficient A(B)/A(i ^ 1 + D^liy - 1) given by Eq. {82) as a function of normalized magnetic field y = B/Bco shown by squai'es 
for YbRh2Si2 and by circles for high-r^ Tl2Ba2Cu06-i-x- is the only fitting parameter. 



Therefore the normalization by the scales Aq and Bcq immediately reveals the common physical nature of these substances, 
allowing us to get rid of the specific properties of the system that define the values of Aq and Bco- 

Based on the above analysis of the A coefficients, we conclude that there is at least one q uantum phase transition inside 
the superconducting dome of high-Tf superconductors, and this transition is FCQPT il20ll . 

8. Appearance of FCQPT in Fermi systems 

We say that Fermi systems that approach QCP from a disordered state are highly correlated systems in order to 
distinguish them from strongly correlated systems (or liquids) that are already beyond FCQPT placed at the quantum 
critical line as shown in Fig. |6l A detailed description of the properties of highly correlated systems are given in Section 
|9] and the properties of strongly correlated systems are discussed in Section [TO] In the present section, we discuss the 
behavior of the effective mass M* as a function of the densit y x of the s ystem as x — > Xfc- 



The experimental facts for high-density 2D ^He |64, 72, 152 , 153 ^ show that the effective mass becomes divergent 



when the value of the density at which the 2D liquid He begins to solidify is reached uM- Also observed was a sharp 
increase in the effective mass in the metallic 2D electron system as the density x decreases and tends to the critical density 
of the metal-insulator transition fisl]. We note that there is no ferromagnetic instability in the Fermi systems under 
consideration and the corresponding Landau amplitude > -1 li65> , ,7Zl . which agrees with the model of nearly localized 
fermionsiHliSlill]. 



We examine the divergence of the effective mass in 2D and 3D highly correlated Fermi liquids at T = as the density 
X — > Xfc approaching FCQPT from the disordered phase. We begin by calculating M* as a function of the difference 
(x - Xfc) for a 2D Fermi liquid. For this, we use the equation for M* derived in ItoIi . where the divergence of M* related 
to the generation of density wave in various Fermi liquids was predicted. As jc — » Xfc, the effective mass M* can be 
approximately written as 



M* m 47r2 j j JT^ 



1 So 

Vi-r 

(83) 
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Here we use the notation pf ^jlil^^y) - q(y), where q(y) is the momentum transfer, v(y) is the pair interaction, the 
integral with respect to the coupUng constant g is taken from zero to the actual value ^o, Xq(Q^ ^) is the linear response 
function for the noninteracting Fermi liquid, and R(q, o)) is the effective interaction, with both functions taken at w = 0. 
The quantities R axidxa determine the response function for the system, 

Xiq, 0), g) = — . (84) 

1 - R{q, CD, g)xo{q. (^) 

Near the instability related to the generation of density wave at the density Xcdw, the singular part of the response function 
X has the well-known form, see e.g. ^ 

X'Uq, w, g) - a{XcdK - x) + b(q - qcf + c(gQ - g), (85) 



where a, b, and c are constants and qc ^ Ipfis the density-wave momentum. Substituting Eq. ( |85] | in SS3i and integrating, 
we can represent the equation for the effective mass M* as 

' ' ' (86) 



M*(x) m yjx - Xcdw ' 

where c is a positive constant. It follows from Eq. ( |86] | that M*(x) diverges as a function of the difference (x - x/rc) and 
M*(x) ^ oo as X ^ x/rc I62l l63 1 

- «! + , (87) 

m X — xpc 

where ai and 02 are constants. We note that Eqs. ( l86b and ( |87] ) do not explicitly contain the interaction v{q), although v(^) 
affects fli, fl2 and Xfc- This result agrees with Eq. ( fT9b , which determines the same universal type of divergence (i.e., a 
divergence that is independent of the interaction). Hence, both Eqs. ( fT9] l and dSTl i can be applied to 2D ^He, the electron 
liquid, and other Fermi liquids. We also see that FCQPT precedes the formation of density waves (or charge-density 
waves) in Fermi systems. 

We note that the difference (x - Xfc) must be positive in both cases, since the density x approaches xpc when the 
system is on the disordered side of FCQPT with the finite effective mass M*(x) > 0. In the case of ^He, FCQPT occurs 
as the density increases, when the potential energy begins to dominate the ground-state energy due to the strong repulsive 
short ranged part of the inter-particle interaction. Thus, for the 2D ^He liquid, the difference (x - Xfc) on the right hand 
side of Eq. ( |87] ) must be replaced with {xfc - x). Experiments have shown that the effective mass indeed diverges at high 



densities in the case of 2D ■'He and at low densities in the case of 2D electron systems ll65Ll72ll . 

In Fig. [13] we report the experimental values of the effective mass M* (z) obtained by the measurements on ^He 



monolayer |72]. These measurements, in coincidence with those from Ref. 111531 115411 . show the divergence of the 
effective mass at x = Xfc- To show, that our FCQPT approach is able to describe the above data, we represent the 
fit of M*(z) by the fractional expression M*{z)/M cc bi + i>2/(l - z) and the reciprocal effective mass by th e linear fi t 



m/M*(z) oc bi z- We note here, that the linear fit has been used to describe the experimental data for bilayer ^He II153LI155 1 



and we use this function here for the sake of illustration. It is seen from Fig. [I3]that the data of Ref. 015311 (^He bilayer) 



can be equally well approximated by both linear and fractional functions, while the data in Ref. 117211 cannot. For instance, 
both fitting functions give for the critical density in bilayer xpc ~ 9.8 nm"^, while for monolayer f/S] these values are 
different - Xfc = 5.56 nm"^ for linear fit and xpc = 5.15 nm"^ for fractional fit. It is seen from Fig. [T3] that linear fit is 
unable to properly describe the experiment ll72ll at small 1 - z (i.e. near x = Xpc)-, while the fractional fit describes the 
experiment very well. This means that more detailed measurements are necessary in the vicinity x — Xfc 

The effective mass as a function of the electron density x in a silicon MOSFET is shown in Fig. [3] We see that 
Eq. ( l86b provides a good description of the e xperi mental results. The divergence of the effective mass M*(x) discovered 



in measurements involving 2D ^He 1164 l72l 115311 is illustrated by Figs. |4]and[T3] Figures [3l l4l and [T3] show that the 



description provided by Eqs. ( fT9] l. ( |86] | and (l87b is in good agreement with the experimental data. 
In the case of 3D systems, as x — > xpc, the effective mass is given by the expression jTOll 

I go 

M* m "^n^ } [\-R{q(y),g)x^{qiy))f' 
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Figure 13: The dependence of the effective mass M'(z) on dimensionless density I - z = 1 - x/Xfc- Experimental data from Ref. |72'] are shown by 
circles and squares and those from Ref. [153] are shown by triangles. The effective mass is fitted as M'(z)/in bi + bil(} - z) (also see Eq. M9\ ). while 
the reciprocal one sls m/M' (z) ^ z, where i , ^2 and are constants . 

Comparison of Eqs. ( l83T l and dSSl l shows that there is no essential difference between them, although they describe different 
cases, 2D and 3D. In the 3D case, we can derive equations similar to ( l86l l and (l87T i just as we did in the 2D case, but the 
numerical coefficients are different, because they depend on the number of dimensions. The only difference between 2D 
and 3D electron systems is that FCQPT occurs in 3D systems at densities much lower than those corresponding to 2D 
systems. No such transition occurs in massive 3D ^He because the transition is absorbed by the first-order liquid-solid 
phase transition 1,64. 72,1 . 

9. A highly correlated Fermi liquid in HF metals 

As noted in the Introduction, the challenge for the theories is to explain the scaling behavior of the normalized effective 
mass M^{y) displayed in Fig. |2j the theories analyzing only the critical exponents that characterize M*^{y) at y » 1 
consider a part of the problem. In this Section we analyze and derive the scaling behavior of the normalized effective mass 
near QCP as depicted in Fig. |2]and show that numerous facts collected in measurements of the thermodynamic, transport 
and relaxation properties carried out at the transition regime on HF metals can be explained within the framework of the 
extended quasiparticle paradigm describing the scaling behavior 

9.1. Dependence of the effective mass M* on magnetic field 

When the system approaches FCQPT from the disordered side, at sufficiently low temperatures as shown in Fig. |6] it 
remains LFL with the effective mass M* that strongly depends on the distance r - {x-Xfc)/xfc and magnetic field B. This 
state of the system with M* that strongly depends on r and B resembles the state of strongly correlated liquid described 
in Sections |4] and [10] But in contrast to a strongly correlated liquid, the system in question involves no temperature 
independent entropy specified by Eq. ( |27] | and at low temperatures becomes LFL with effective mass M* oc 1/r [see 
Eqs. ( fT9] l and dSTlll. Such a liquid can be called a highly correlated liquid; as we see shortly, its effective mass exhibits 
the scaling behavior. We study this behavior when the system transits from its LFL to NFL states. 

We use the Landau equation to study the behavior of the effective mass M*{T,B) as a function of the temperature 
and the magnetic field. For the model of homogeneous HF liquid at finite temperatures and magnetic fields, this equation 
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acquires the form ll20ll 



1 1 v r pfp 

(T,B) m pI 



M*(T,B) m p^^ 

dn^,(p,T,B) dp 



dp (27r)3 



(89) 



where fa-,cri(PF, p) is the Landau amplitude dependent on momenta Pf, p and spin cr. For the sake of definiteness, we 
assume that the HF Uquid is 3D Hquid. As seen in Section[8] the scaHng behavior calculated within the model of HF liquid 
does not depend on dimensionality and on the inter-particle interaction, while the values of scales like M*^ and Tm do 
depend. To simplify matters, we ignore the spin dependence of the effective mass, because M*{T, B) is nearly independent 
of the spin in weak fields. The quasiparticle distribution function can be expressed as 



n^(p, r) = < 1 + exp 



(e(p, T) - fia-) 



(90) 



where e(p, T) is determined by Q. In our case, the single-particle spectrum depends on the spin only weakly, but 
the chemical potential may depend on the spin due to the Zeeman splitting. When this is important, we specifically 
indicate the spin dependence of physical quantities. We write the quasiparticle distribution function as no-(p, T, B) = 
5no-(p, r, B) + no-(p, T - Q,B - Q). Equation ( l89b then becomes 



m 

+ 



^ y r PfPi 

pI^J pf 



M*iT,B) M*(x) pI 

d6n^S9uT,B) dpi 

X /^^,^,(Pf,Pi) 7 ■prii- ^91) 

dpi {2ny 

We assume that the highly correlated HF liquid is close to FCQPT and the distance r — > 0, and therefore M/M*{x) — > 0, 
as follows from Eq. ( fT9l ). For normal metals, where the electron liquid behaves like LFL with the effective mass of several 
bare electron masses M* /m ~ 1, at temperatures even near 1000 K, the second term on the right hand side of Eq. ( |9TT l 
is of the order of T^/i-i^ and is much smaller than the first term. The same is true, as can be verified, when a magnetic 
field B ~ 100 T is applied. Thus, the system behaves Uke LFL with the effective mass that is actually independent of the 
temperature or magnetic field, while p(T) oc AT^. This means that the corrections to the effective mass determined by the 
second term on the right-hand side of Eq. ( |9TT i are proportional to (T/p.)^ or (psB/fif'. 

Near QCP Xfc, with m/M*{x — > Xfc) — » 0, the behavior of the effective mass changes dramatically because the first 
term on the right-hand side of Eq. ( IgTT i vanishes, the second term becomes dominant, and the effective mass is determined 
by the homogeneous version of Eq. (|9TJ as a function of B and T. As a result, the LFL state vanishes and the system 
demonstrates the NFL behavior down to lowest temperatures. 

We now qualitatively analyze the solutions of Eq. ( |9l1 i at x ^ Xfc and 7=0. Application of magnetic field leads 
to Zeeman splitting of the Fermi surface, and the distance Sp between the Fermi surfaces with spin up and spin down 
becomes 6p - p^. - Pp ~ PbBM*(B)I pf. We note that the second terr n on t he right-hand side of Eq. W\\ is proportional 
to {5pf oc (pBBM*(B)/pFf, and therefore Eq. (|9B reduces to ll49i Isl [Til 

{pBBM'{B)f 

+ c 7 , (92) 



M*{B) M'ix) p4 

where c is a constant. We also note that M*{B) depends on x and that this dependence disappears at x = x fc- A t this 



point, the term m/M*{x) vanishes and Eq. ( |92] i becomes homogeneous and can be solved analytically 05311631113811 



^*^^^-(^3^- (93) 



where B^o is the critical magnetic field, regarded as a parameter (see remarks to Fig. [TTT i. 



35 



Equation ( |93] ) specifies the universal power-law behavior of the effective mass, irrespective of the interaction type and 
is valid in 3D and 2D cases. For densities x > Xfc, the effective mass M*(x) is finite and we deal with the LFL state if the 
magnetic field is so weak that 

^ « 1, (94) 
M*iB) 

with M*{B) given by Eq. ( l93T l. The second term on the right-hand side of Eq. (|92] |. which is proportional to {BM*(x))^, is 
only a small correction. In the opposite case, at T/T*{B) 1, where 

^ » 1, (95) 

the electron liquid behaves as if it were at the quantum critical point. In the LFL state, the main thermodynamic and 
transport properties of the system are determined by the effective mass. It therefore follows from Eq. ( |93] l that we 
have the unique possibility of controlling the magnetoresistance, resistivity, heat capacity, magnetization, thermal bulk 
expansion, etc by varying the magnetic field B. It must be noted that a large effective mass leads to a high density of 
states, which causes the emergence of a large number of competing states and phase transitions. We believe that such 
states can be suppressed by applying an external magnetic field, and we examine the thermodynamic properties of the 
system without considering such competition. 

9.2. Dependence of the effective mass M* on temperature and the damping of quasiparticles 

For a qualitative examination of the behavior of M*(T,B,x) as the temperature increases, we simplify Eq. ( l9Tb 
by dropping the variable B and by imitating the effect of an external magnetic field by a finite effective mass in the 
denominator of the first term on the right hand side of Eq. ( |9TT l. Then the effective mass becomes a function of the distance 
r, M*(r), determined also by both the magnitude of the magnetic field B and x. In a zero magnetic field, r = (x- Xfc)/xfc, 
We integrate the second term on the right-hand side of Eq. (|9TT l with respect to the angular variables, then integrate by 
parts with respect to p, and replace p with z = {s{p) - fj)/T. In the case of a flat and narrow band, we can use the 
approximation where {s{p) - jj) - prip - Pf)/M*(T). The result is 

M m F(pf,pf(l+az))dz 



Jo 



M*{T) M*ir) Jo 1 + e< 



Jo 



Fipf,pFil-az))- :. (96) 

1 4- e- 



where we use the notation 



d(F^p^) TM*(T) TM'(T) 
m- — — 



dp ' pI iT,MHr)y 

Tk - pl/M*{r), and the Fermi momentum is defined by the condition s(pf) = yU. 

We first consider the case where a <K 1. Then, discarding terms of the order exp(-l/Q'), we can set the upper limit 
in the second integral on the right hand side of Eq. ( l96b to infinity, with the result that the sum of the second and third 
terms is an even function of a. The resulting integrals are typical expre ssion s involving the Fermi-Dirac function in the 
integrand and can be evaluated by a standard procedure (e.g., see Ref. il56ll V Because we need only an estimate of the 
integrals, we write Eq. ( |96] ) as 



+ ai 



TM*(T) 



M*{T) M*(r) \TkM*{r) 

TM'iTW 



+ 02 —\ +-, (97) 

where a\ and 02 are constants of the order of unity. 

Equation Wh is a typical equation of the LFL theory. The only exception is the effective mass M*{r), which depends 
strongly on the distance r and diverges as r — > 0. Nevertheless, Eq. ( |97] ) implies that as T — > 0, the corrections to M*(r) 
begin with terms of the order if 

'tm*{tW t'- 



» — - —7, (98) 
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and the system behaves like LFL. Condition ( |98T l impUes that the behavior inherent in LFL disappears as r — > and 
M*(r) — > oo. Then the free term on the right-hand side of Eq. ( |96l ) is negligible, m/M*{r) — > 0, and Eq. ( |96] l becomes 
homogeneous and determines the universal behavior of the effective mass M*{T). At a certain temperature T* <K Tk, 
the value of the sum on the right-hand side of Eq. ( l97T i is determined by the second term and relation ( l98l l becomes 
invalid. Keeping only the second term in Eq. (|97] |. we arrive at an equation for determining M*{T) in the transition region 



I53LI157I1: 

^'^^)'^^- ^99) 

As regards an estimate of the transition temperature T*(B) at which the effective mass becomes temperature- 
dependent, we note that the effective mass is a continuous function of the temperature and the magnetic field; M*(B) ~ 
M*(r*). With Eqs. (|93]l and we obtain 

T*(B) ^ fiB(B - Bm). (100) 

As the temperature increases, the system transfers into another mode. The coefficient a is then of the order of unity, a ~ 1, 
the upper limit in the second integral in Eq. (|96] l cannot be set to infinity, and odd terms begin to play a significant role. 
As a result, Eq. (|97j breaks down and the sum of the first and second integrals on the right-hand side of Eq. ( |96] l becomes 
proportional to M*(T)T. Ignoring the first term m/M*{r) and approximating the sum of integrals by M*{T)T, we obtain 
from (|96]l that 

M*(T) oc-^. (101) 

We note that M*{T) is also given by Eq. ( llOll l if the Landau amplitude F{p) is determined by a nonanalytic function, that 
is the function cannot be expanded in Tailor series at /? — > 0, see SectionfTSl 

We therefore conclude that as the temperature increases and the condition x ^ Xfc is satisfied, the system demonstrates 
regimes of three types: (i) the behavior of the Landau Fermi liquid ata I, when Eq. ( |98] l is valid and the behavior of the 
eff'ective mass is described by Eq. ( |93] l; (ii) the behavior defined by Eq. ( |99] l. when M*{T) cc 7" 2/3 ^j^^j g ^j-^ M*(T)T oc 
r^/^; and (iii) the behavior at a ~ 1, when Eq. ( fTOTT l is valid, M'(T) oc 1/ Vf, the entropy S(T) cc M*{T)T oc Vf, and 
the heat capacity C{T) = T(dS(T)/dT) oc ^/f. 




Figure 14: The entropy S (7") of a highly correlated Fermi liquid at the critical point of FCQPT. The solid line represents the function S (T) oc VT, and 
the squares mark the results of calculations. 
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We il lustrate the behavior of S (T) when Eq. dlOlb is valid using calculations based on the model Landau functional 
I4l[l58l] 



with the nonanalytic Landau amplitude 



n(Pi)«(P2)^^, (102) 



V(p) = .o^^e^. (103) 
IPI 



We normalize the effective mass to m, i.e., M* - M* /m, and the temperature To to the Fermi energy e*^, T - Tq/s'-^, and 
use the dimensionless coupling constant g = (gom)/(2n^) and also /3 = (^oPf- FCQPT occurs when these parameters reach 
the critical values /? = and g - gc- On the other hand, a transition of this kind occurs as M* — > oo. This condition 



allows deriving a relation between be and gc 11421 115811 : 



^{\+bc) exp(-Z7 J [be coMbc) - sinh(fec)] = 1 . ( 1 04) 

be 

This relation implies that the critical point of FCQPT can be reached by varying go if A) ^"d pf are fixed, by varying 
Pf if /3o and go are fixed, etc. For definiteness, we vary g to reach FCQPT or to study the properties of the system 
beyond the critical point. Calculations of M*{T), S (T), and C{T) based on the model functional ( 1102b with the parameters 
g - gc - 6.7167 and /5 - be - 3 show that M*(T) oc 1/ Vt, S(T) oc y/T, and C(T) oc ^/T. The temperature dependence 
of the entropy in this case is depicted in Fig. [14] 

We now estimate the quasiparticle damping yiT). In the Landau Fermi-liquid theory, y(T) is given by 1I20I1 

\3t2 



y ~ \rnMyT\ (105) 



where T is the particle-hole amplitude. In the case of highly correlated HF system with a high density of states caused by 
the enormous effective mass, F cannot be approximated by the "bare" interaction between particles but can be estimated 
within the "ladder" approximation, which yields |r| ~ I /{pfM*{T)) 115 U 15211 . As a result, we find that yiT) oc T^ 



in the Landau Fermi-liquid regime since M* is temperature-independent. Then, y{T) oc T'^^^ in the T ^^^-regime, and 
y{T) oc r3/2 in the 1/ Vf -regime. We note that in all these cases, the width is small compared to the characteristic 
quasiparticle energy, which is assumed to be of the order of T, and hence the quasiparticle concept is meaningful. 

The conclusion that can be drawn here is that when the HF liquid is localized near QCP of FCQPT and is on the 
disordered side, its low-energy excitations are quasiparticles with the effective mass M*{B, T).We note that at FCQPT, 
the quasiparticle renormalization z-factor remains approximatel y constan t and the divergence of the effective mass that 
follows from Eq. ( fT9b is not related to the fact that z — > |23, 24, 159ll . Therefore, the quasiparticle concept remains 
valid and it is these quasiparticles that constitute the extended paradigm and determine the transport, relaxation and 
thermodynamic properties of HF liquid. 



9.3. Scaling behavior of the effective mass 

As was mentioned in the Introduction, to avoid difficulties associated with the anisotropy generated by the crystal 
lattice of HF metals, we study the universal behavior of HF metals using the model of the homogeneous HF (electron) 
liquid. The model is quite meaningful because we consider the scaling behavior exhibited by the effective mass at low 
temperatures. The scaling behavior of the effective mass is determined by energy and momentum transfers that are small 
compared to the Debye characteristic temperature and momenta of the order of the reciprocal lattice cell length a ^ 
Therefore quasiparticles are influenced by the crystal lattice averaged over large distances compared to the length a. 
Thus, we can use the well-known jelly model. We note that the values of such scales as M^, Tm, Bcq and Z?c2 etc depend 
on the properties of a HF metal, its lattice, composition etc. For example, the critical magnetic field Beo depends even on 
its orientation with respect to the lattice. 

To explore the scaling behavior of M*, we start with qualitative analysis of Eq. ( l89b . At FCQPT the effective mass 
M* diverges and Eq. ( [89] l becomes homogeneous determining M* as a function of temperature as given by Eq. ( |99] ). If 
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the system is located before FCQPT, M* is finite, and at low temperatures the system demonstrates the LFL behavior, that 
is M*{T) ^ M* + aiT^. As we have seen in Subsection 19.21 the LFL behavior takes place when the second term on the 
right hand side of Eq. ( [89] l is small in comparison with the first one. Then, at increasing temperatures the system enters 
the transition regime: M* grows, reaching its maximum M^, ?AT — Tm, with subsequent diminishing. Near temperatures 
T > Tm the last "traces" of LFL regime disappear, the second term starts to dominate, and again Eq. ( |89] | becomes 
homogeneous, and the NFL behavior restores, manifesting itself in decreasing of M* as T^^^^. When the system is near 
FCQPT, it turns out that the solution of Eq. ( [89l ) M*(T) can be well approximated by a simple universal interpolating 
function Il60ll . The interpolation occurs between the LFL (M* ^ M* + a^T^) and NFL (M* oc T'^^^) regimes, thus 
describing the above crossover. Introducing the dimensionless variable y -T^ - T/Tm, we obtain the desired expression 

M*^(y) ^ co-i^^. (106) 

1 + C2y'^'^ 

Here M*^ - M* jM*^ is the normalized effective mass, cq = (1 +C2)/(1 +ci), c\ and C2 are fitting parameters, parameterizing 
the Landau amplitude. 

It follows from Eq. ( |93] l, that the application of magnetic field restores the LFL behavior, so that M^, depends on B as 

Ml,^{B-B,^,)-^l\ (107) 

while 

Tm^MB-B,o). (108) 

Employing Eqs. (1107b and dlOSb to calculate M*^ and Tm, we conclude that Eq. ( 1106b is valid to describe the normalized 
effective mass in external fixed magnetic fields with y = T/{B - Bco). On the other hand, Eq. (1106b is valid when the 
applied magnetic field becomes a variable, while temperature is fixed at T = Tf. In that case, it is convenient to represent 
the variable as y = (B - Bco)/Tf. 



9.3.1. Schematic phase diagram ofHF metal 

The schematic phase diagram of HF metal is reported in Fig. [15] panel a. Magnetic field B is taken as the control 
parameter In fact, the control parameter can be pressure P or doping (the number density) x etc as well. At B = Bcq, 
FCQPT takes place leading to a strongly degenerate state, where Bcq is a critical magnetic field, such that at B > Bed the 
system is driven towards the LFL state. We recall, that in our simple model Bcq is a parameter. The FC state is captured 
by the superconducting (SC), ferromagnetic (FM), antiferromagnetic (AFM) etc. states lifting the degeneracy. Below in 
Subsection|93]we consider the HF metal YbRh2Si2. In that case, Bcq ^ 0.06 T (B±c) and at T = and B < B^ the AFM 
state takes place with temperature dependent resistivity p{T) oc T^ 11511 . At elevated temperatures and fixed magnetic 
fields the NFL regime occurs, while rising B again drives the system from the NFL state to the LFL one as shown by the 
dash-dot horizontal arrow in Fig. [15] We consider the transition region when the system moves from the NFL state to LFL 
one along the horizontal arrow and also moves from LFL state to NFL one along the vertical arrow as shown in Fig. [15] 
The inset to Fig. [T5]demonstrates the scaling behavior of the normalized effective mass M]^ - M* jM*^ versus normalized 
temperature Tt^ - T/Tm, where is the maximum value that M* reaches at T = Tm- The T^^^^ regime is marked as 
NFL since the effective mass depends strongly on temperature. The temperature region T ^ Tm signifies the crossover 
(or the transition region) between the LFL state with almost constant effective mass and NFL behavior, given by T^^^^ 
dependence. Thus temperatures T ~ Tm can be regarded as the crossover region between the LFL and NFL states. 

The transition (crossover) temperature T*(B) is not really the temperature of a phase transition. It is necessarily broad, 
very much depending on the criteria for determination of the point of such a crossover, as it is seen from the inset to Fig. 
T5]a. see e.g. Refs. 1.15^,150.1 . As usual, the temperature T*{B) is extracted from the field dependence of charge transport, 
for example from the resistivity p{T) given by Eq. (ITsT l. The LFL state is characterized by the T"'" dependence of the 
resistivity with - 2, see Subsection 19.51 The crossover (that is the transition regime shown by the hatched area both 
in the panel a of Fig. [15] and in its inset) takes place at temperatures where the resistance starts to deviate from the LFL 
behavior with - 2 so that the exponent becomes 1 < q-r < 2, see Subsection [93] As it will be shown in Subsection 
931 in the NFL state or - I. 

The panel b of Fig. [TS] represents the experimental T — B phase diagram of the exponent a^iT, B) as a function of 
temperature T versus magnetic field B ||7|l. The evolution of aniT, B) is shown by the color: the yellow color corresponds 
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Magnetic field B ^^^^ 

Figure 15: The panel a represents a schematic phase diagram of HF metals. B^o is magnetic field at which the effective mass divergences. SC, FM, AFM 
denote the superconducting (SC), ferromagnetic (FM) and antiferromagnetic (AFM) states, respectively. At B < B^o the system can be in SC, FM or 
AFM states. The vertical arrow shows the transition from the LFL to the NFL state at fixed B along T with M* depending on T. The dash-dot horizontal 
arrow illustrates the system moving from the NFL to LFL state along B at fixed T. The exponent a/t determines the temperature dependent part of the 
resistivity, see Eq. )18t . In the LFL state the exponent = 2 and in the NFL o-j; = 1. In the transition regime the exponent evolves from its LFL value 
to the NFL one. The inset shows a schematic plot of the normalized eft'ective mass versus the normahzed temperature. Transition regime, where 
reaches its maximum value at 7 = Tm, is shown by the hatched area both in the panel a and in the inset. The arrows mark the position of inflection 
point in Mj^ and the transition region. The panel b shows the experimental T - B phase diagram of the exponent ai;(T, B) as a function of temperature 
T versus magnetic field B |7]. The evolution of (1-r(7',B) is shown by the color: the yellow color coiresponds to ai((T,B) = 1 (the NFL state) and the 
blue color corresponds to aR(T,B) = 2 (the LFL state). The NFL behavior occurs at the lowest temperatures right at QCP tuned by magnetic field. At 
rising magnetic fields B > B^o and T ~ T'(B), the broad transition regime from the NFL state to the field-induced LFL state occurs. As in the panel a, 
the both transitions from the LFL to the NFL state and from the NFL to LFL state are shown by the corresponding arrows. 



to afi(r, B) - I and the blue color corresponds to a^iT, B) - 2. It is seen from the panel that at the critical field 
BcQ - 0.66 T {B\\c) the NFL behavior occurs down to the lowest temperatures. While YbRhaSia transits from the NFL to 
LFL behavior under the application of magnetic field. It is worthy to note that the phase diagram displayed in Fig. [T5](the 
panel a) coincides with that of shown in the panel b. 

A few remarks are in order here. As we shall see, the magnetic field dependence of the effective mass or of other 
observable like the longitudinal magnetoresistance do not have "peculiar points" like maxima. The normalization is to be 
performed in the other points like the inflection point alT - Ti„f (or at B = shown in the inset to Fig. [T5]by the 
arrow. Such a normalization is possible since it is established on the scales, r,„y cx Tm <x (B — Bco)- As a result, we obtain 

fiB(Binf-B,.o)^Tf. (109) 

It follows from Eq. (1106b that in contrast to the Landau paradigm of quasiparticles the effective mass strongly depends on 
T and B. This dependence leads to the extended quasiparticle paradigm and forms the NFL behavior. Also from Eq. ( 1106b 
the scaling behavior of M* near QCP is revealed by the application of appropriate physical scales to measure the effective 
mass, magnetic field and temperature. At fixed magnetic fields, the characteristic scales of temperature and of the function 
M*{T, B) are defined by both Tm and M*^ respectively. At fixed temperatures, the characteristic scales are (Bm - Bco) and 
M^. From Eqs. ( 1107b and (1108b it is seen that at fixed magnetic fields, Tm — > 0, and — > oo, and the width of the 
transition region shrinks to zero as B ^ Bco when these are measured in "external" scales like T in K, magnetic field B in 
T etc. In the same way, it follows from Eqs. ( |99] l and ( |109b that at fixed temperature Tf, {Bi„f - Bco) — > 0, and M^^ — > oo, 
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and the width of the transition region shrinks to zero as Tf — > 0. Thus, the application of the external scales obscure the 
scaling behavior of the effective mass and the thermodynamic and transport properties. 




Figure 16: The normalized entropy SaiCB/B,,,/) versus y = B/Bi„f and the normalized entropy S n{T /Ti„f) versus v = T/Ti„f calculated at fixed 
temperature and magnetic field, correspondingly, are represented by the solid lines and shown by the arrows. The inflection point is depicted by the 
dash-dot arrow. 



In what follows, we compute the effective mass using Eq. ( l89l l and employ Eq. ( 11061) for estimations of the considered 
values. To compute the effective mass M*{T,B), we solve Eq. ( l89b with a quite general form of Landau interaction 
amplitude ifssll . Choice of the amplitude is dictated by the fact that the system has to be at QCP, which means that the first 
two /7-derivatives of the single-particle spectrum e(p) should equal zero. Since the first derivative is proportional to the 
reciprocal quasiparticle effective mass l/M*, its zero just signifies QCP of FCQPT. The second derivative must vanish; 
otherwise sjp) - has the same sign below and above the Fermi surface, and the Landau state becomes unstable before 
r 112311 1 371] . Zeros of these two subsequent derivatives mean that the spectrum e(p) has an inflection point at pf so 
that the lowest term of its Taylor expansion is proportional to (p - ppf ■ After solution of Eq. ( |89] |. the obtained spectrum 
has been used to calculate the entropy 5 (B, T), which, in turn, has been used to recalculate the effective mass M*{T, B) by 
virtue of the well-known LFL relation M*(T, B) - S(T, B)/T. Our calculations of the normalized entropy as a function 
of the normalized magnetic field B/Bi„f - y and as a function of the normalized temperature y = T/Ti„f are reported in 



Fig. 



Here T,,,/ and B,„/ are the corresponding inflection points in the function S . We normalize the entropy by its 



value at the inflection point S N(y) - S(y)/S(l). As seen from Fig. [16] our calculations corroborate the scaling behavior 
of the normalized entropy, that is the curves at different temperatures and magnetic fields merge into a single one in terms 
of the variable y. The inflection point T,,,/ in S (T) makes M*(T, B) have its maximum as a function of T, while M*(T, B) 
versus B has no maximum. We note that our calculations of the entropy confirm the validity of Eq. (1106b and the scaling 
behavior of the normalized effective mass shown in Fig. [15] 



9.4. Non-Fermi liquid behavior in YbRh2Si2 

In this Subsection, we analyze the transition regime and the NFL behavior of the HF metal YbRh2Si2. We demonstrate 
that the NFL behavior observed in the thermodynamic and transport properties of YbRh2Si2 can be described in terms 
of the scaling behavior of the normalized effective mass. This allows us to construct the scaled thermodynamic and 
transport properties extracted from experimental facts in a wide range of the variation of scaled variable and conclude that 
the extended quasiparticles paradigm is strongly valid. We show that "peculiar points" of the normalized effective mass 
give rise to the energy scales observed in the thermodynamic and transport properties of HF metals. Our calculations 
of the thermodynamic and transport properties are in good agreement with the heat capacity, magnetization, longitudinal 
magnetoresistance and magnetic entropy obtained in remarkable measurements on the heavy fermion metal YbRh2Si2 
1 1^ [iTl |36l llZll • For YbRh2Si2 the constructed thermodynamic and transport functions extracted from experimental facts 
show the scahng over three decades in the variable. The energy scales in these functions are also explained BSll . 
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9.4.1. Heat capacity and the Sommerfeld coefficient 

Exciting measurements of C IT oc M* on samples of the new generation of YbRh2Si2 in different magnetic fields B up 
to 1.5 T [36] allow us to identify the scaling behavior of the effective mass M* and observe the different regimes of M* 
behavior such as the LFL regime, transition region from LFL to NFL regimes, and the NFL regime itself. A maximum 
structure in C/J oc M^^ at Tm appears under the application of magnetic field B and Tm shifts to higher T as B is increased. 
The value of CjT - yo is saturated towards lower temperatures decreasing at elevated magnetic fields. 

The transition region corresponds to the temperatures where the vertical arrow in the main panel a of Fig. [T5] crosses 
the hatched area. The width of the region, being proportional to Tm oc (B - Bco) shrinks, Tm moves to zero temperature 
and 70 M* increases as Z? — » Bcq. These observations are in accord with the experimental facts ISf 
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Figure 17: The normalized eft'ective mass extracted from the measurements of the specific heat C/T on YbRh2Si2 in magnetic fields B listed 
in the legend. Our calculations are depicted by the solid curve tracing the scaling behavior of . 

To obtain the normalized effective mass M*^, the maximum structure in C/T was used to normalize C/T, and T was 
normalized by Tm- In Fig. [T7] Mt as a function of normaUzed temperature T^ is shown by geometrical figures, our 
calculations are shown by the solid line. Figure [T7]reveals the scaling behavior of the normalized experimental curves - 
the scaled curves at different magnetic fields B merge into a single one in terms of the normalized variable y = T /Tm- As 
seen, the normalized mass extracted from the measurements is not a constant, as would be for LFL. The two regimes 
(the LFL regime and NFL one) separated by the transition region, as depicted by the hatched area in the inset to Fig. [15] 
a, are clearly seen in Fig. [T7] displaying good agreement between the theory and experimental facts. It is worthy to note 
that the normalization procedure allows us to construct the scaled function C/T extracted from the experimental facts in 
wide range variation of the normalized temperature. Indeed, it integrates measurements of C/T taken at the application of 
different magnetic fields into unique function of the normalized temperature which demonstrates the scaling behavior over 
three decades in the normalized temperature as seen from Fig. [17] As seen from Figs. [T] the NFL behavior extends at least 
to temperatures up to few Kelvins. Thus, we conclude that the extended quasiparticle paradigm does take into account 
the remarkably large temperature ranges over which the NFL behavior is observed. We note that at these temperatures the 
contribution coming from phonons is still small. 

9.4.2. Magnetization 

Consider now the magnetization M as a function of magnetic field B at fixed temperature T - Tj 

M{B,T)^ f xib,T)db, (110) 
Jo 
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where the magnetic susceptibihty x is given by ll20ll 



X(B,T) 



PM*(B, T) 

1 + F" 
i -I- r 



(111) 



Here, is a constant and is the Landau amplitude related to the exchange interaction. In the case of strongly correlated 
systems > -0.9 16111 . Therefore, as seen from Eq. ( Ill II ). due to the normalization the coefficients p and (1 +Fg) 

drops out from the result, andx °^ M* . One might assume that F^ can strongly depend on B. This is not the case [i38i[39fl . 
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Figure 18: The field dependencies of the normalized magnetization M collected at dift'erent temperatures shown at light bottom corner are extracted 
from measurements collected on YbRu2Si2 flTllSTll . The kink (shown by the arrow) is clearly seen at the normalized field Bn = BjBt — 1. The solid 
curve represents our calculations. 



since the Kadowaki-Woods ratio is conserved lfl4ll27 , 141 1. A(B)/yf^(B) oc A(B)/x^(B) oc const, we have yo M* oc )(. 



Note that the Sommerfeld coefficient does not depend on F^. 

Our calculations show that the magnetization exhibits a kink at some magnetic field B - B^. The experimental 
magnetization demonstrates the same behavior We use Bk and MiBk) to normalize B and M respectively. 

The normalized magnetization M(B)/M{Bk) extracted from experimental facts depicted by the geometrical figures and 
calculated magnetization shown by the solid line are reported in Fig. [18] As seen, the scaled data at different Tf merge 
into a single one in terms of the normalized variable y - B/Tt- It is also seen, that these exhibit energy scales separated 
by kink at the normalized magnetic field B^ - BjBk - 1. The kink is a crossover point from the fast to slow growth 
of M at rising magnetic field. Figure [18] shows that our calculations are in good agreement with the experimental facts, 
and all the data exhibit the kink (shown by the arrow) ?A B^ ^ 1 taking place as soon as the system enters the transition 
region corresponding to the magnetic fields where the horizontal dash-dot arrow in the main panel a of Fig. [lacrosses the 
hatched area. Indeed, as seen from Fig. [18] at lower magnetic fields M is a linear function of B since M* is approximately 
independent of B. Then, Eqs. (1106b and ( 11071 ) show that at elevated magnetic fields M* becomes a diminishing function 



of B and generates the kink in M{B) separating the energy scales discovered in Refs. Il7il37ll . It is seen from Eq. ( |109l l 
that the magnetic field Bk at which the kink appears, B^ - Bm oc Tf, shifts to lower B as Tf is decreased. This observation 
is in accord with experimental facts 

Consider now an "average" magnetization M = Bx + M as a function of the magnetic field B at fixed temperature 
T - Tf\\l]. We again use Bk and M_{Bit) to normalize B and M_ respectively. The normalized M_ vs the normalized field 
Bff - B/Bk are shown in Fig. [19] Our calculations are depicted by the solid line. The stars trace out our calculations of 
M with M*{y) extracted from the data C/T shown in Fig. [17] It is seen from Fig. [T9]that our calculations are in good 
agreement with the experimental facts, and all the data exhibit the kink (shown by arrow) at B^ - 1 taking place as soon 
as the system enters the transition region corresponding to the magnetic fields where the horizontal dash-dot arrow in the 
main panel a of Fig. [lacrosses the hatched area. Indeed, as seen from Fig. [19] at lower magnetic fields M is a hnear 
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Figure 19: The field dependence of the normalized "average" magnetization M = M+Bx is shown by squares and has been extracted from measurements 
collected on YbRuTSi? 1 17]. The kink (shown by the arrow) is clearly seen at the normalized field Bn = B/Bt ^ 1. The solid curve and stars (see text) 
represent our calculations. 



function of B since M* is approximately independent of B. It follows from Eq. (1107b that at elevated magnetic fields 
M* becomes a diminishing function of B and generates the kink in M_{B) separating the energy scales discovered in Ref. 
UtIi . Then, it seen from Eq. il09\ that the magnetic field Bk — Bm at which the kink appears shifts to lower B as Tf is 
decreased. 

9.4.3. Longitudinal magnetoresistance 

Consider a longitudinal magnetoresistance (LMR) p(B, T) - pQ+ AT^ as a function of B at fixed Tf. In that case, the 
classical contribution to LMR due to o rbital motion of carriers induced by the Lorentz force is small, while the Kadowaki- 
Woods relation lH l27l [l39[ [l40[ 141 ^. K - A/Jq cc A/x^ - const, allows us to employ M* to construct the coefficient A, 



since yo cc ^ cc M*. Omitting the classical contribution to LMR, we obtain that p{B, T) - po cc (M*)^. Fig. |20]reports the 
normalized magnetoresistance 

p^iy)^ ^^^^'^' ^(M*^(y)f (112) 

Pinf 

versus normalized magnetic field y = B/Bi„f at different temperatures, shown in the legend. Here p,„/ and B,„/ are LMR 
and magnetic field respectively taken at the inflection point marked by the arrow in Fig. |20] Both theoretical (shown by the 
solid line) and experimental (marked by the geometrical symbols) curves have been normalized by their inflection points, 
which also reveal the scaling behavior - the scaled curves at different temperatures merge into a single one as a function 
of the variable y and show the scaling behavior over three decades in the normalized magnetic field. The transition region 
at which LMR starts to decrease is shown in the inset to Fig. \T5\a by the hatched area. Obviously, as seen from Eq. (1109b . 
the width of the transition region being proportional to Bm - Bj„f o^Tj decreases as the temperature T j is lowered. In the 
same way, the inflection point of LMR, generated by the inflection point of M* shown in the inset to Fig. [T5]by the arrow. 



shifts to lower B as is decreased. All these observations are in good agreement with the experimental facts II17U37I1 . 



9.4.4. Magnetic entropy 

The evolution of the derivative of magnetic entropy dS(B, T)/dB as a function of magnetic field B at fixed temper- 
ature Tf is of great importance since it allows us to study the scaling behavior of the derivative of the effective mass 
TdM*{B, T)/dB cc dS(B, T)/dB. While the scaling properties of the effective mass M*{B, T) can be analyzed via LMR, 
see Fig. |20] As seen from Eqs. ( 1106b and (1 109b . at y < 1 the derivative 

dM^iy) 
oc y 

dy 
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Normalized magnetic field 

Figure 20: Magnetic field dependence of the normalized magnetoresistance pN versus normalized magnetic field, pjv was extracted from LMR of 
YbRh2Si2 at dift'erent temperatures fnilBTll Hsted in the legend. The inflection point is shown by the arrow, and the sohd line represents our calculations. 




Figure 21: Normalized magnetization difl'erence divided by temperature increment (AM/AT)n vs n ormalized magnetic field at fixed temperatures 
(listed in the legend in the upper left corner) is extracted from the facts collected on YbRh2Si2 II6III . Our calculations of the normalized deiivative 
(dS /dB)fj ^ {AM/AT)fj vs normalized magnetic field are given at fixed dimensionless temperatures T/fi (listed in the legend in the upper right comer). 
All the data are shown by the geometrical figures depicted in the legend at the upper left comer. 



with y = (B - Bco)/(Bj„f - Bcq) oc (B - Bco)/Tf. We note that the effective mass as a function of B does not have the 
maximum. At elevated y the derivative -dMN(y)/dy possesses a maximum at the inflection point and then becomes a 
diminishing function of y. Upon using the variable y = (B - BcQ)/Tf, we conclude that at decreasing temperatures, the 
leading edge of the function -dS jdB oc -TdM* /dB becomes steeper and its maximum at (B,„/ - Bco) oz Tf is higher. 
These observations are in quantitative agreement with striking measurements of the magnetization difference divid ed b y 
temperature increment, -AM /AT, as a function of magnetic field at fixed temperatures Tf collected on YbRh2Si2 lll6l[l . 
We note that according to the well-known thermodynamic equality dM/dT = dS /dB, and AM /AT ^ dS /dB. To carry 
out a quantitative analysis of the scaling behavior of -dM*(B, T)/dB, we calculate as described above the entropy S{B, T) 
shown in Fig. [T6]as a function of B at fixed dimensionless temperatures Tj/^ shown in the upper right corner of Fig. [21] 
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This figure reports tlie normalized (dS /dB)^ as a function of the normalized magnetic field. The function (dS /dB)i^ is 
obtained by normalizing ( -dS / dB) by its maximum taking place at B m, and the field B is scaled by Bm- The measurements 
of -AM /AT are normalized in the same way and depicted in Fig. |2T|as (AM/AT)^ versus normalized field. It is seen 
from Fig. [2T|that our calculations are in good agreement with the experimental facts and both the experimental functions 
{AM/AT)n and the calculated (dS /dB)^ show the scaling behavior over three decades in the normalized magnetic field. 

9.4.5. Energy scales 




Figure 22: Temperature versus magnetic field T - B phase diagram for YbRh2Si2- Solid circles represent the boundary between AF and NFL states. 
The solid squares denote the boundary of the NFL and LFL regime IT^ ITtIIstIi shown by the dashed line which is approximated by " ^cO idl- 
Diamonds mark the maximums Tm of C/T 1, 1 89.1 shown in Fi g. 1171 T he dash-dot line is approximated by Tm a(B - Bco), a is a fitting parameter, see 
Eq. 1108) . Triangles along the solid line denote Ti„f in LMR I17l[37ll sown in Fig. [21] and the solid line represents the function Ti„f oc h(B - B^o), b is 
a fitting parameter, see Eq. (T09). 



Fig. |22]reports T,-,,/ and Tm versus B depicted by the solid and dash-dotted lines, respectively. The boundary between 
the NFL and LFL regimes is shown by the dashed line, and AF marks the antiferromagnetic state. The corresponding data 
are taken from Ref. 1117113711 16 ill . It is seen that our calculations are in good agreement with the experimental facts. In 
Fig. |22] the solid and dash-dotted lines corresponding to the functions Tinf an d Tm, respectively, represent the positions of 
the kinks separating the energy scales in C and M reported in Ref. lll7lll6lll . Furthermore, our calculations are in accord 
with experimental facts, and we conclude that the energy scales are reproduced by Eqs. ( llOSI l and (|109l l and related to the 
peculiar points T,,,/ and Tm of the normalized effective mass M]^ which are shown by the arrows in the inset to Fig. [15] 

At B — > BcQ both Tj„f — » and Tm — » 0, thus the LFL and the transition regimes of both C/T and M as well as 
those of LMR and the magnetic entropy are shifted to very low temperatures. Therefore due to experimental difficulties 
these regimes cannot often be observed in experiments on HF metals. As it is seen from Figs. [T7] [18] [20] |2T] and 
|22] the normalization allows us to construct the unique scaled thermodynamic and transport functions extracted from 
the experimental facts in a wide range of the variation of the scaled variable y. As seen from the mentioned Figures, 
the constructed normalized thermodynamic and transport functions show the scaling behavior over three decades in the 
normalized variable. 



9.5. Electric resistivity of HF metals 

The electric resistivity of strongly correlated Fermi systems, p{T) - po + Ap\{B, T), is determined by the effective 
mass , because o f the Kadowaki- Woods relation Api(B, T) = A{B, T)T^ oc {M*{B, T)T)^, see Subsection |9.4.3| and Refs. 

IM [M [mH, 

and therefore the temperature dependence of the effective mass discussed above can be observed in 
measurements of the resistivity of HF metals. At temperatures T <K T*(B), the system is in the LFL state, the behavior of 
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the effective mass as x ^ xfc is described by Eq. ( |93] ). and the coefficient A(B) can be represented as 



A(B) cc 



1 



(113) 



In this regime, the resistivity behaves as Api - ciT^/{B - Bcof^^ cc T^. The second regime, a highly correlated Fermi 
liquid determined by Eq. (l99b. is characterized bv the resistivitv dependence Api = C2T^ /(T^^^)^ cc T^l^. The third regime 
at r > T*{B) is determined by Eq. ( IIOII ). In that case we obtain Api - cj,T^ /(T^^^)^ oc T. If the system is above the 
quantum critical line as shown in Fig. |6] the dependence of the effective mass on temperature is given by Eq. ( |3TI ). so we 
obtain from Eq. ( llOSI l that the q uasiparticle damping y{T) cc T IsTll . As a result, we see that the resistivity dependence on 
temperature is Api = C4T 114711 . Here, ci, C2, C3 and C4 are constants. If the system at the transition regime, as shown by 
the arrows in Fig. [15] the dependence of the effective mass on temperature cannot be characterized be a single exponent 
as it is clearly seen from the inset to Fig. [T5]a. So we have that Api oc T"" with 1 < qr < 2. We note that all temperature 
dependencies corresponding to these regimes have been observed in measurements involving the heavy-fermion metals 
CeCoIns, YbRh2Si2 and YbAgGe 1 15, 3d.l3l [l62i[l63ll . 



9.6. Magnetic susceptibility and magnetization measured on CeRu2Si2 

Experimental investigations of the magnetic properties of CeRu2Si2 down to the lowest temperatures (down to 170 
mK) and ultrasmall magnetic fields (down to 0.21 mT) have shown neither evidence of the magnetic ordering, super- 
conductivity nor conventional LFL behavior foS]. These results imply a magnetic quantum critical point in CeRu2Si2 is 
absent and the critical field Bco = 0. Even if the magnetic quantum critical point were there it should maintain the NFL 
behavior over four decades in temperature. Such a strong influence can hardly exist within the framework of conventional 
quantum phase transitions. 

Temperature dependence on a logarithmic scale of the normalized AC susceptibility xiB, T) is shown at different 
applied magnetic fields B as indicated in the left panel of Fig. |23] versus normalized temperature. The right panel of 
the Figure shows the normalized static magnetization Mb{B, T) (DC susceptibility) in the same normalized temperature 
range. The temperature is normalized to Tm (the temperature at which the susceptibility reaches its peak value), the 
susceptibility is normalized to the peak value x{B, Tm), and the magnetization is normalized to Mb{B, T 0), for each 
value of the field [92]. If we use Eq. (II 10b and the definition of susceptibility (II 1 lb . we conclude that the susceptibility 
and magnetization also demonstrate the scaling behavior and can be represented by the universal function (1106b of the 
single variable y, if they are respectively normalized as discussed above. We see from Fig. |23]that at finite field strengths 




Normalized temperature 



Figure 23: The normalized magnetic susceptibility ^(B, T)/x{B, Tm) (the left panel) and normalized magnetization Mb(B, T)/Mb(B, Tm) (DC suscep- 
tibility, the right panel) for CeRu2Si2 in magnetic fields 0.20 mT (squares), 0.39 mT (triangles), and 0.94 mT (circles) as functions of the normalized 
temperature T/Tm t92ll . The solid lines depict the calculated scaling behavior i53h as described in Subsection l9.3.1l 



B, the curve descri bing x(B, T)/x(B, Tm) has a peak at a certain temperature Tm, whil e Mb (B, T)/Mb{B, Tm) has no such 
peak ll49[ Issi 157ll . This behavior agrees well with the experimental results ll49i Issl Il57ll obtained in measurements on 
CeRu2Si2 ll92ll . We note that such behavior of the susceptibility is not typical of ordinary metals and cannot be explained 
within the scope of theories that take only ordinary quantum phase transitions into account Ii92il . 
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Figure 24: Temperature dependence of the AC susceptibility XAC for CeRu2Si2. The solid curve is a fit for the data shown by the triangles at B = 0.02 
mT I92I1 and represented by the function xC^) = 'j/ given by Eq. jlOlt with a being a fitting parameter. Inset shows the normalized effective mass 
versus normalized temperature T^j extracted from^^c measured at different fields as indicated in the inset |92]. The solid curve traces the universal 
behavior of M'^{Tfi) determined by Eq. )106t . Parameters ci and C2 are adjusted to fit the average behavior of the normalized effective mass M^. 



To verify Eq. ( llOlb and illustrate the transition from LFL behavior to NFL one, we use measurements of Xac(T) 
in CeRu2Si2 at magnetic field B = 0.02 mT at which this HP metal demonstrates the NFL behavior down to lowest 
temperatures I92i1 . Indeed, in this case we expect that LFL regime to start to form at temperatures lower than Tm ~ fJ^gB ~ 
0.01 mK as it follows from Eq. dlOSb . It is seen from Fig. |24]that Eq. dlOll ) gives good description of the experimental 
facts in the extremely wide range of temperatures: the susceptibility ;\fAC as a function of T, is not a constant upon cooling, 
as would be for a Fermi liquid, but shows a 1/ Vt divergence over almost four decades in temperature. The inset of Fig. 
|24]exhibits a fit for M*^ extracted from measurements of xac(T) at different magnetic fields, clearly indicating the change 
from LFL behavior at Tai < 1 to NFL one at Tat > 1 when the system moves along the vertical arrow as shown in Fig. 
[TS] It seen from Figs. |23] and |24] that the function given by Eq. ( IIO6I 1 represents a good approximation for M*^ within the 
extended paradigm. In Subsection l9.4l we have seen that the same is true in the case of YbRhiSii with the AF quantum 
critical point. We conclude that both alloys CeRu2Si2 and YbRh2Si2 demonstrate the universal NFL thermodynamic 
behavior, independent of the details of the HE metals such as their lattice structure, composition and magnetic ground 
state. This conclusion implies also that numerous QCPs related to conventional quantum phase transitions assumed to 
be responsible for the NFL behavior of different HF metals can be well reduced to a single QCP related to FCQPT and 
accounted for within the extended quasiparticle paradigm 1.164,1 . 



9.7. Transverse magnetoresistance in the HF metal CeCoIns 

Our comprehensive theoretical study of both the longitudinal and transverse magnetoresistance (MR) shows that it is 
(similar to other thermodynamic characteristics like magnetic susceptibility, specific heat, etc) governed by the scaling 
behavior of the quasiparticle effective mass. The crossover from negative to positive MR occurs at elevated temperatures 
and fixed magnetic fields when the system transits from the LFL behavior to NFL one and can be well captured by this 
scaling behavior. 

By definition, MR is given by 

p{B,T)-piO,T) 

pMT)^ ' ^114) 
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We apply Eq. ( II 14l i to study MR of strongly correlated electron liquid versus temperature T as a function of magnetic 
field B. The resistivity p(B, T) is 

p(B, T) = po + Ap(B, T) + ApdB, T), (115) 
where po is a residual resistance, Ap = ciAT^, ci is a constant. The classical con tribution ApiiB, T) to MR due to 



orbital motion of carriers induced by the Lorentz force obeys the Kohler's rule 1116511 . We note that Api{B) <s: p(0, T) 



as it is assumed in the weak-field approximation. To calculate A, we again use the quantities yo - C/T cc M* and/or 
X M* as well as employ the fact that the Kadowaki-Woods ratio K = A/y^ cc A/x^ = const. As a result, we obtain 
A cc (M*)^, so that Ap{B, T) = c(M*(B, T))^T^ and c is a constant. Suppos e tha t the temperature is not very low, so that 
po < Ap{B = 0, T), and B > B^. Substituting (fTTsT i into (fTT4l i. we find that fl6^ 

po + ApdB,T) .{M*{B,T)f-{M'{0,T)f 

Pmr - JTTTf: + cT — — . (116) 

p(0, T) p(0, T) 



Consider the qualitative behavior of MR described by Eq. ( II 16b as a function of B at a certain temperature T - Tq. 
In weak magnetic fields, when the system exhibits NFL (see Fig. [TSl l. the main contribution to MR is made by the term 
Api{B), because the eff'ective mass is independent of the applied magnetic field. Hence, \M*{B, T)-M*{Q, T)\/M*{0, T) <K 
1 and the leading contribution is made by Api{B). As a result, MR is an increasing function of B. When B becomes 
so high that T*{B) ~ hb{B - Bm) ~ To, the difference (M*(B, T) - M*(0, T)) becomes negative because M*(B, T) is 



now the diminishing function of B given by Eq. ( 1107b . Thus, MR as a function of B reaches its maximal value at 
T*{B) ~ Tn(B) ~ To. At further increase of magnetic field, when Tm{B) > To, the effective mass M*{B, T) becomes a 
decreasing function of B. As B increases, 

(M*(B,T)-M*(0,T)) ^ ^ 

and the magnetoresistance, being a decreasing function of B, can reach its negative values. 

Now we study the behavior of MR as a function of T at fixed value Bo of magnetic field. At low temperatures T <s: 
T*(Bo), it follows from Eqs. (fT06l l and (|93]l that M*(Bo, r)/M*(0, T) « 1, and it is seen from Eq. ( fTTTl l that p,„ABo, T) ~ 
-1, because Api(Bo, T)/p(0, T) 1. We note that Bo must be relatively high to guarantee that M*(Bo, T)/M*(0, T) <K 1. 
As the temperature increases, MR increases, remaining negative. At T ^ T*(Bo), MR is approximately zero, because 
p(Bo, T) ^ p(0, T) at this point. This allows us to conclude that the change of the temperature dependence of resistivity 
p(Bo, T) from quadratic to linear manifests itself in the transition from negative to positive MR. One can also say that 
the transition takes place when the system goes from the LFL behavior to the NFL one. At T > T*{Bo), the leading 
contribution to MR is made by Api{Bo, T) and MR reaches its maximum. At Tm{Bo) T, MR is a decreasing function 
of the temperature, because 

\M^iB,T)-M^(0,T)\^ (118) 

M*(o,r) ^ ^ 

and p,„r{Bo, T) <K 1 . Both transitions (from positive to negative MR with increasing B at fixed temperature T and from 
negative to positive MR with increasing T at fixed B value) have been detected in measurements of the resistivity of 
CeColns in a magnetic field [30]. 

Let us turn to quantitative analysis of MR 1 166]. As it was mentioned above, we can safely assume that the classical 
contribution ApiiB, T) to MR is small as compared to Ap(B, T). Omission of Api(B, T) allows us to make our analysis 
and results transparent and simple while the behavior of Apiifio, T) is not known in the case of HE metals. Consider 
the ratio = p(B, T)/p(Q, T) and assume for a while that the residual resistance po is small in comparison with the 
temperature dependent terms. Taking into account Eq. (II 15b and p(0, T) cc T, we obtain from Eq. (II 16b that 

RP = p„r + 1 = T(M*(B, T))\ (1 19) 

P(0, T) 

and consequently, from Eqs. ( 1106b and (II 19b that the ratio R^ reaches its maximal value at some temperature T^^ ~ 
Tm- If the ratio is measured in units of its maximal value and T is measured in units of rR„, ~ Tm then it is seen from 
Eqs. ( 1106b and ( |119b that the normalized MR 



(120) 




Normalized temperature 

Figure 25: The normalized magnetoresistance fi^(v) given by Eq. jllOj ^ versus normalized temperature v = r/TR,,,. ^C*') was extracted from MR 
shown in Fig. I27l and collected on CeCoIns at fixed magnetic fields B ISO] listed in the right upper comer. The starred line represents our calculations 
based on Eqs. 11061 and iilO) with the parameters extracted from AC susceptibility of CeRu2Si2 (see the caption to Fig. 1241 . The solid Hne displays 
our calculations based on Eqs. jl2U and U20t : only one parameter was used to fit the data, while the other were extracted from the AC susceptibility 
measured on CeRu2Si2. 



becomes a function of the only variable y = T/T^^. To verify Eq. ( I120I I. we use MR obtained in measurements on 
CeCoIns, see Fig. 1(b) of Ref. OOfl . The results of the normalization procedure of MR are reported in Fig. |25] It is clearly 
seen that the data collapse into the same curve, indicating that the normalized magnetoresistance obeys the scaling 



behavior well given by Eq. ( I120l i. This scaling behavior obtained directly from the experimental facts is a vivid evidence 
that MR behavior is predominantly governed by the effective mass M*(B, T). 

Now we are in position to calculate R^^^iy) given by Eq. ( 11201 ). Using Eq. ( IIO6I 1 to parameterize M*^(y), we extract 



parameters ci and C2 from measurements of the magnetic AC susceptibility ;t' on CeRu2Si2 ll92ll and apply Eq. ( 11201 ) to 



calculate the normalized ratio. It is seen that the calculations shown by the starred line in Fig. |25] start to deviate from 



experimental points at elevated temperatures. To improve the coincidence, we employ Eq. ( llOlb which describes the 
behavior of the effective mass at elevated temperatures and ensures that at these temperatures the resistance behaves as 
p(T) oc T. To correct the behavior of M*^{y) at rising temperatures M* ~ T"'^^, we add a term to Eq. ( 11061 ) and obtain 



1 + ci/ _^ exp(-l/y) 
1 + C2y^/^ ^ y/y 



(121) 



where C3 is a parameter The last term on the right hand side of Eq. (11211 ) makes satisfy Eq. (llOll l at temperatures 
T/Tm > 2. In Fig. |25] the fit of ^^(y) by Eq. (1121b is shown by the solid line. Constant C3 is taken as a fitting parameter, 
while the other were extracted from AC susceptibility of CeRu2Si2 as described in the caption to Fig. l24l 

Before discussing the magnetoresistance Pmr{B, T) given by Eq. ( II 14b . we consider the magnetic field dependence of 
both the MR peak value RmnxiB) and the corresponding peak temperature Tr^{B). It is possible to use Eq. (11 19b which 
relates the position and value of the peak with the function M*{B,T). Since Tf^^ oc fiB{B - Bco), B enters Eq. dl 19b 
only as tuning parameter of QCP, as both Api and po were omitted. At B — > Bco and T <K Tji,m(B), this omission is not 
correct since Api and po become comparable with Ap(B, T). Therefore, both Rm-dxiB) and 7'Rm(B) are not characterized by 
any critical field, being continuous functions at the quantum critical field Bco, in contrast to M*(B, T) whose peak value 
diverges and the peak temperature tends to zero at Bco as seen from Eqs. (11071 ) and ( 1108b . Thus, we have to take into 
account ApiiB, T) and po which prevent rRm(B) from vanishing and make Rm-dxiB) finite at B — > Bco- As a result, we have 
to replace B^o by some effective field Bgff < Bco and take B^jf as a parameter which imitates the contributions coming 
from both ApiiB, T) and po. Upon modifying Eq. (II 191 ) by taking into account Api(B, T) and po, we obtain 
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Figure 27: MR versus temperature T as a function of magnetic field B. The experimental data on MR were collected on CeCoIns at fixed magnetic field 
B 13011 shown in the right bottom comer of the Figure. The solid lines represent our calculations, Eq. 1106) is used to fit the eftective mass entering Eq. 

{So). 



njB) ^ bi(B - B,ff), (122) 

b2(B - BeffT^I^ - 1 
i?max(B) - ■ ' "\ I , . (123) 

biiB - Beff)-^ + 1 

Here b\, bi, bj and Bg/f are the fitting parameters. It is pertinent to note that while deriving Eq. ( 1123b we use Eq. 
( II22I 1 with substitution (B - Be//) for T. Then, Eqs. (1122b and ( 1123b are not valid at B < Bm- In Fig. |26] we show 
the field dependence of both TRn, and /?niax, extracted from measurements of MR ll30ll . Clearly both and Rm^x are 
well described by Eqs. ( 1122b and (1123b with Bgff =3.8 T. We note that this value of B^ff is in good agreement with 
observations obtained from the T - B phase diagram of CeCoIns, see the position of the MR maximum shown by the 
filled circles in Fig. 3 of Ref. Q. 
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To calculate pmi (B, T), we apply Eq. (11201 ) to describe its universal behavior, Eq. ( IIO6I 1 for the effective mass along 
with Eqs. ( 1122b and ( I123l l for MR parameters. Figure |27] shows the calculated MR versus temperature as a function of 
magnetic field B together with the experimental points from Ref. |30]. We recall that the contributions coming from 
ApiiB, T) and po were omitted. As seen from Fig. |27] our description of the experiment is good. 

9.8. Magnetic-field-induced reentrance of Fermi-liquid behavior and spin-lattice relaxation rates in YbCu5_xAUx 

One of the most interesting and puzzling issues in the research on HE metals is their anomalous dynamic and relaxation 
properties. It is important to verify whether quasiparticles with effective mass M* still exist and determine the physical 
properties of the muon and ^^Cu nuclear spin-lattice relaxation rates l/Ti in HE metals throughout their temperature - 
magnetic field phase diagram, see Eig. [15] This phase diagram comprises both LEL and NEL regions as well as NEL- 
LEL transition or the crossover region, where magnetic-field-induced LEL reentrance occurs. Measurements of the muon 
and ^^Cu nuclear spin-lattice relaxation rates 1 IT\ in YbCu4 4AU0.6 have shown that it differs substantially from ordinary 
Eermi liquids obeying the Korringa law 1 167]. Namely, it was reported that for T — > reciprocal relaxati on ti me diverges 



as 1 / Ti r oc r following the behavior predicted by the self-consistent renormalization (SCR) theory 16811 . The static 



unifo rm susceptibility x diverges as oc r so that l/TiT scales with x^- Latter result is at variance with SCR theory 



1 167I1 . Moreover, the application of magnetic field B restores the LEL behavior from initial the NEL one, significantly 



reducing l/Ti. These experimental fi ndings are hard to explain within both the conventional LEL approach and in terms 



of other approaches like SCR theory Ill67lll68ll . 

In this Subsection we show that the above anomalies along with magnetic-field-induced reentrance of LEL properties 
are indeed determined by the dependence of the quasiparticle effective mass M* on magnetic field B and temperature T 
and demonstrate that violations of the Korringa law also come from M*(B, T) dependence. Our theoretical analysis of 
experimental data on the base of ECQPT approach permits not only to explain the above two experimental facts in a unified 
manner, but to unveil their universal properties, relating the peculiar features of both longitudinal magnetoresistance and 
specific heat in YbRh2Si2 to the behavior of spin-lattice relaxation rates. 

To discuss the deviations from the Korringa law in light of NEL properties of YbCu4 4AU0.6, we notice that in LEL 
theory the spin-lattice relaxation rate l/T[ is determined by the quasiparticles near the Fermi level. The above relaxation 
rate is related to the decay amplitude of the quasiparticles, which in turn is proportional to the density of states at the 
Eermi level N{Ef). Formally, the spin-lattice relaxation rate is determined by the imaginary paitx" of the low-frequency 
dynamical magnetic susceptibility ;k'(q, w — » 0), averaged over momentum q 



where Aq is the hyperfine coupling constant of the muon (or nuclei) with the spin excitations at wave vector q il68ll . If 
Aq = A() is independent of q, then standard LEL theory yields the relation 

= 7tAIn\Ef). (125) 

Equation ( I125I I can be viewed as Korringa law. Since in our ECQPT approach the physical properties of the system under 
consideration are determined by the effective mass M*{T, B, x), we express 1 /Ti T in Eq. (1125b via it. This is accomplished 
with the standard expression 1I20I1 NiEp) = M* ppjn^, rendering Eq. (1125b to the form 

= ^M*2 = 1^ [M'{T, B, x)f , (126) 

where rj - (Aq/?^)/:/!^ =const. The empirical expression 

T^^xHn (127) 



extracted from experimental data in YbCu5_xAux lll67ll . follows expUcitly from Eq. ( 1126b and well-known LEL relations 

M* ccxocC/T. 
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In what follows, we compu te the effective mass as it was explained in Subsection 19.3.11 and employ Eq. ( 1106b for 
estimations of obtained values [169]. The decay law given by Eq. ( |99] l along with Eq. ( 1126b permits to express the 
relaxation rate in this temperature range as 

^^a,+a2T-'l' <^x\T), (128) 

where ai and a2 are fitting parameters. The dependence (1128b is reported in Fig. l28l a long with experimental points for the 
muon and nuclear spin-lattice relaxation rates in YbCu4 4Auo.6 at zero magnetic field il67ll . It is seen from Fig. |28]that Eq. 
( 1128b gives good description of the experiment in the extremely wide temperature range. This means that the extended 
paradigm is valid and quasiparticles survive in close vicinity of FCQPT, while the observed violation of Korringa law 
comes from the dependence of the effective mass on temperature. 
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Figure 28: Temperature dependence of muon (squares) and nuclear (circles) spin-lattice relaxation rates (divided by temperature) for YbCu4 4Auo.6 at 
zero magnetic field 1167]. The solid curve represents our calculations based on Eq. (128). 

Figure |29] displays magnetic field dependence of normalized muon spin-lattice relaxation rate 1/7"^^ in YbCu5_xAux 
(x=0.6) along with our theoretical B-dependence. To obtain the latter theoretical curve we (for fixed temperature and in 
magnetic field B) employ Eq. (1126b and solve the Landau integral equation to calculate M*(T,B) as it was described in 
Subsection |9.3.1| We note that the normalized effective mass M*^(y) was obtained by normalizing M*{T, B) at its infection 
point shown in the inset to Fig. [15] 

It is instructive to compare the LMR analyzed in Subsection l9.4.3l and l/T^. LMR p(B, T) - po + Pb + MB, T)T^ is 
as a function of B at fixed T, where po is the residual resistance, ps is the contribution to LMR due to orbital motion of 
carriers induced by the Lorentz force, and A is the coefficient. As we see in Subsection |9.4.3l ps is small and we omit 
this contribution. The Kadowaki-Woods relation allows us to employ M* to calculate A(B, T). As a result, p(B, T)- po cc 
(M*f, and 1/r^^ oc (M'f as seen from Eq. (1126b . As a result, we see that that LMR and the magnetic field dependence 
of normalized muon spin-lattice relaxation rate 1 /7"^^ can be evaluated from the same equation 

^(3') = = ^ - (M^Miy))'- (129) 

Inset to Fig. |29]reports the normalized LMR vs normalized magnetic field y - BIBm at different temperatures, shown 
in the legend. Here pM and Bint are respectively LMR and magnetic field taken at the inflection point. The inflection 
points of both LMR and 1 /Tiai are generated by the inflection point of M* shown in the inset to Fig. [15] a by the arrow. 
The transition region where LMR starts to decrease is shown in the inset by the hatched area and takes place when the 
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Figure 29: Magnetic field dependence of normalized at the inflection point muon spin-lattice relaxation rate 1 IT'^^j extracted from measurements Il67ll 
on YbCu44Auo.6 along with our calculations of B-dependence of the quasiparticle effective mass. Inset shows the normalized LMR f^Cv) versus 
normalized magnetic field. S^Cv) was extracted from LMR of YbRh2Si2 at dift'erent temperatures fnll Usted in the legend. The solid curves represent 
our calculations. 



system moves along the horizontal dash-dot arrow. We note that the same normalized effective mass has been used to 
calculate both l/T^^ in YbCu44Auo.6 and the normalized LMR in YbRh2Si2. Thus, Eq. ( I129l l determines the close 
relationship between the quite different dynamic properties, showing the validity of the quasiparticle extended paradigm. 
In Fig. |29] both theoretical and experimental curves have been normalized by their inflection points, which also reveals 
the scaling behavior - the curves at different temperatures merge into a single one in terms of the scaled variable y. Figure 
|29l shows clearly that both normalized magnetoresistance 7?^ and reciprocal spin-lattice relaxation time obey well the 
scaling behavior given by Eq. (1129b . This fact obtained dkectly from the experimental findings is vivid evidence that the 
behavior of both the above quantities is predominantly governed by the field and temperature dependence of the effective 
mass. 

We remark that the same normalized effective mass determines the behavior of the thermodynamic and transport 
properties in YbRh2Si2, see Subsection 19.41 It is seen from the Figures presented in Subsection 19.41 that our calculations 
of the effective mass offer good descriptions of such different quantities as the relaxation rates (l/TiT) and the transport 
(LMR) and thermodynamic properties in such different HE metals asYbCu5_xAux and YbRh2Si2. It is pertinent to note 
that the obtained good description makes an impressive case in favor of the reliability of the quasiparticle extended 
paradigm. 



9.9. Relationships between critical magnetic fields Bm and Bc2 in HF metals and high-Tc superconductors 

Recently, in high-Tc superconductors, exciting measurements revealing their physics have been performed. One type 
of th e measurements demonstrate the existence of Bogoliubov quasiparticles (BQ) in the superconducting state 1 100.,, 



type 

11 g 



1 ITll . While in the pseudogap regime at temperatures above Tc when the superconductivity vanishes, a strong indication 



of the pairing of electrons or the formation of preformed pairs of electrons was observed, while the gap continues to 
follow the simple d-wave form 1 1 16., 1 171 . Another type of the measurement explored the normal state induced by the 
application of magnetic field, when the transition from the NFL behavior to LFL one occurs lll50ll . As we have mentioned 
in Subsection l7.2.1l there are the experimental relationships between the critical fields Bc2 > Bco, where Bc2 is the field 
destroying the superconducting state, and Bm is the critical field at which the magnetic field induced QCP takes place. Now 
we show that Bc2 > BcO- We note that to study the aforementioned transition experimentally in high-Tc superconductors, 
strong magnetic fields of B > Bc2 are required; earlier, such investigation was technically inaccessible. An attempt to 
study the transition experimentally had already been made I148ll . 
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Figure 30: T ~ B phase diagram of the superconductor Tl2Ba2Cu06+j:. The crossover (from LFL to NFL regime) line T'(B) is given by the Eq. i fTOSt . 
Open squares and solid circles ai'e experimental values |150]. Thick line represents the boundary between the superconducting and normal phases. 
Arrows near the bottom left corner indicate the critical magnetic field Ed destroying the superconductivity and the critical field B^o- Inset displays the 
peak temperatures rn,nx(B)> extracted from measurements of CjT and^^c on YbRh2(Sio.95Geo.o5)2 iTL IlVOll and approximated by straight lines Eq. 
(108). The lines intersect at B ^ 0.03 T. 



Let us now consider the T - B phase diagram of the high-Tc superconductor Tl2Ba2Cu06+jc shown in Fig. |30] The 
substance is a superconductor with Tc from 15 K to 93 K, being controlled by oxygen content fl50^. In Fig. 1301 open 
squares and solid circles show the experimental values of the crossover temperature from the LFL to NFL regimes [15^- 
The solid line given by Eq. (llOOI l shows our fit w ith Bco - 6T that is in good agreement with Bco = 5.8 T obtained from 
the field dependence of the charge transport fl5(f\ . 

As it is seen from Fig. |30] the linear behavior agrees well with the experimental data 1 150i 166ll . The peak temper ature s 
Tmax shown in the inset to Fig. 1301 depict the maxima of C{T)/T and xac(T) measured on YbRh2(Sio.95Geo.o5)2 
From Fig. [30l r^ax is seen to shift to higher values with increase of the applied magnetic field and both functions can be 
represented by straight lines intersecting at B ^ 0.03 T. This observation is in good agreement with experiments |7, 17(J]. 
Clearly from Fig. [30]the critical field Bc2 = 8 T destroying the superconductivity is close to Bco = 6 T. We now show that 
this is more than a simple coincidence, and Bc2 ^ BcO- Indeed, at B > B^o and low temperatures T < T*(B), the system is 
in its LFL state. The superconductivity is then destroyed since the superconducting gap is exponentially small as we have 
seen in Subsection l5.3l At the same time, there is the FC state at B < Bco and this low-field phase has large prerequisites 
towards superconductivity as in this case the gap is a linear function of the superconducting coupling constant /lo as it was 
shown in Subsection [53] We note that this is exactly the case in CeCoIns where Bco - Bc2 - 5 T |30] as seen from Fig. 
[3n while the application of pressure makes Bc2 > BcO ll33ll . However, if the superconducting coupling constant is rather 
weak then antiferromagnetic order wins the competition. As a result, Bc2 - 0, while B^ can be finite as in YbRhiSii and 
YbRh2(Sio.95Geo.o5)2 115, 170]. 

Comparing the phase diagram of Tl2Ba2Cu06+.i with that of CeCoIns shown in Figs. |30] and |3T] respectively, it is 
possible to conclude that they are similar in many respects. Further, we note that the superconducting boundar y lin e 
Bc2iT) at decreasing temperatures acquires a step, i.e. the corresponding phase transition becomes first order ll97iTl7l. 
This leads us to speculate that the same may be true for Tl2Ba2Cu06+.v We expect that in the NFL state the tunneling 
conductivity is asymmetrical function of the applied voltage, while it becomes symmetrical at the a pplic ation of increased 
magnetic fields when Tl2Ba2Cu06+.v transits to the LFL behavior, as it predicted to be in CeCoIns M172ll . 

It follows from Eq. (ISTb that it is impossible to observe the relatively high values of A{B) since in our case Bc2 > BeO- 
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Figure 31: T - B phase diagram of tlie CeCoInj heavy femiion metal. The interface between the superconducting and normal phases is shown by the 
solid line to the squai'e where the phase transition becomes a first-order phase transition. AtT < Tq, the phase transition is a first-order phase transition 
I171II . The interface between the superconducting and normal phases is shown by the dashed line. The solid straight line represented by Eq. U08t with 
the experimental points |31] shown by squares is the interface between the LFL and NFL states. 



We note that Eq. (ISTT i is valid when the superconductivity is destroyed by the application of magnetic field, otherwise 
the efifective mass is also finite being given by Eq. (HTt . Therefore, as was mentioned above, in high-Tc QCP is poorly 
accessible to experimental observations being "hidden in superconductivity". Nonetheless, thanks to the experimental 
facts 1 150], we have seen in Subsection 17 .2. 1 1 that it is possible to study QCP of high-Tc 1,120.1 . As seen from Fig. [I2l 
the facts give evidences that the physics underlying the field-induced reentrance of LFL behavior, is the same for both HF 
metals and hight-Je superconductors. 

9.70. Scaling behavior of the HF CePdi^yJRh^ferrotnagnet 

QCP can arise by suppressing the transition temperature T^l of a ferromagnetic (FM) (or antiferromagnetic (AFM)) 
phase to zero by tuning some control parameter ^ othe r than tem perature, such as pressure P, magnetic field B, or doping 
X as it takes place in the HF ferromagnet CePdi_xRhx inlliTl or the HF metal Celnj-xSn^ [ 175]. 



The HF metal CePdi_xRhx evolves from ferromagnetism at x = to a non-magnetic state at some critical concen- 
tration Xfc- Utilizing the extended quasiparticle paradigm picture and the concept of FCQPT, we address the question 
about the NFL behavior of the ferromagnet CePdi_xRhx and show that it coincides with that of the antiferromagnets 
YbRh2(Si().95Geo.o5)2 and YbRh2Si2, and paramagnets CeRu2Si2 and CeNi2Ge2. We again conclude that the NFL behav- 
ior, being independent of the peculiarities of a specific alloy, is universal. Incidentally, numerous quantum critical points 
assumed to be responsible for the NFL behavior of different HF metals can be well reduced to the only quantum critical 
point related to FCQPT il76ifl77ll . 

As we have seen above, the efifective mass M*{T,B) can be measured in experiments on HF metals. For example, 
M*{T,B) DC C(T)/T oc a{T)/T and M*{T,B) Xac{T) where XAciT) is ac magnetic susceptibility. If the corresponding 
measurements are carried out at fixed magnetic field B (or at fixed both the concentration x and B) then the effective mass 
reaches its maximum at some temperature Tm- Upon normalizing both the effective mass by its peak value at each field B 
and the temperature by Tm, we observe that all the curves merge into a single one, given by Eq. (1106b . thus demonstrating 
a scaling behavior. 

It is seen from Fig. |32] that the behavior of the normalized ac susceptibility ;^'^(-,(3') - xac(T /Tm, B)/xac(^, B) - 
M*j^{Tn) obtained in measurements on the HF paramagnet CeRu2Si2 [92.] agrees with both the approximation given by 
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Figure 32: Normalized magnetic susceptibility Xn(Tn, B) = XAciT /Tm, B)/xac(K B) = M'^(Tn) for CeRu2Si2 in magnetic fields 0.20 mT (squares), 
0.39 mT (upright triangles) and 0.94 mT (circles) versus normalized temperature = T/Tu 1 92]. The susceptibiHty reaches its maximum ;(^/ic(7'i^, B) 
at r = Tm- The normalized specific heat (C(Tn)ITn)IC(\) of the HF ferromagnet CePdi-^Rhx with .v = 0.8 versus Tn is shown by downright triangles 
Il74ll . Here Tu is the temperature at the peak of C(T)/T. The solid curve traces the universal behavior of the normalized effective mass determined by 
Eq. jl06t . Parameters ci and C2 are adjusted for;^'jv(7'jv, B)at B = 0.94 mT. 



Eq. ( IIO6I 1 and the normalized specific heat (C(Tn)/Tn)/C(1) - M'^^iT^) obtained in measurements on CePdi_xRhx lll74ll . 
Also, from Fig. |32] we see that the cur ve gi ven by Eq. (1106b agrees perfectly with the measurements on CeRu2Si2 whose 
electronic system is placed at FCQPT 1116411 . 




Figure 33: The normalized effective mass at elev ated magnetic fields as a function of y = T/Tm- The mass taken from the specific heat C/T of the 
HF ferromagnet CePdi_xRhx with x = 0.8 (Ref. Il74ll ) is shown at different magnetic fields B depicted at the right upper corner. At B > IT, M^(3') 
coincides with that of CeRu2Si2 (solid curve, see the caption to Fig. 132) . The normalized specific heat C(>')/C(7'm) of CePdi-xRh.i at different magnetic 
fields B is shown in the inset. The kink in the specific heat is clearly seen at =: 2. The soUd curve represents the function yM^(y) with parameters ci 
and C2 adjusted for the magnetic susceptibility of CeRu2Si2 at B = 0.94 mT. 
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Now we consider the beha vior o f M'^{T), extracted from measurements of the specific heat on CePdi_xRhx under the 
application of magnetic field lll74ll and shown in Fig. |33] It is seen from Fig. |33]that for B > 1 T describes the 
normalized specific heat almost perfectly, coinciding with that of CeRu2Si2 and is in accord with the universal behavior of 
the normalized effective mass given by Eq. ( 1106b . Thus, we conclude that the thermodynamic properties of CePdi_xRhx 
with X - 0.8 are determined by quasiparticles rather than by the critical magnetic fluctuations. On the other hand, one 
could expect the growth of the critical fluctuations contribution as jc ^ Xfc so that the behavior of the normalized effective 
mass would deviate from that given by Eq. ( 11061 ). This is not the case as observed from Fig. |33] It is also seen that at 
increasing magnetic fields B all the curves corresponding to the normalized effective masses extracted from CePdi_xRhx 
with X - 0.8 merge into a single one, thus demonstrating a scaling behavior in accord with equation (IIO6I 1. We note 
that existing theories based on the quantum and thermal fluctuations predict that ma gnetic and thermal p roperties of the 
ferromagnet CePdi_xRhx differ from those of the paramagnet CeRuzSia, Refs. iilH [Hi [Hi \vm . Clearly, from 
the inset of Fig. |33] there is the kink in the temperature dependence of the normalized specific heat C{Tn)/C{Tm) of 
CePdi_xRhx appearing at - 2. In the inset, the solid line depicts the function TnM*^{Tn) with parameters ci and C2 
which are adjusted for the magnetic susceptibility at B = 0.94 mT. Since the function T^M'^^iTfi) describes the normalized 
specific heat very well and its bend (or kink) comes from the crossover from the LFL regime to the NFL one, we safely 
conclude that the kink emerges at temperatures when the system transits from the LFL behavior to the NFL one. As shown 
in Subsection 19.71 the magnetoresistance changes from positive values to negative ones at the same temperatures. One 
may speculate that there is an energy scale which could make the kink coming from fluctuations of the order parameter 
ilTll . In that case we must to concede that such different HF metals as CePdi_xRhx, CeRu2Si2 and CeCoIns with different 
magnetic ground states have the same fluctuations which exert coherent influence on the heat capacity, susceptibility and 
transport properties. Indeed, as we have seen above and will also see below in this Subsection, that Eq. (1106b allows us to 
describe quantitatively all the mentioned quantities. 
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Figure 34: Same as in Fig. I33l but x = 0.85 [1741 . At B > IT, M*^(Tn) demonstrates tlie universal behavior (solid curve, see the caption to Fig. I32t . 



In Fig. [34l the effective mass M*^(Tn) at fixed B's is shown. Since the curve shown by circles and extracted from 
measurements at Z? = does not exhibit any maximum down to 0.08 K 1, 1741 . we conclude that in this case x is very close 
to Xfc and the maximum is shifted to very low temperatures. As seen from Fig. [34l the application of magnetic field 
restores the scaling behavior given by Eq. ( 1106b . Again, this permits us to conclude that the thermodynamic properties of 
CePdi-xRhx with x = 0.85 are determined by quasiparticles rather than by the cr itica l magnetic fluctuations. 

The thermal expansion coefficient a{T) is given by a{T) ^ M*T / (pj,K) il56l l. The compressibility K(p) is not 
expected to be singular at FCQPT and is approximately constant [180]. Taking into account Eq. ( llOlb . we find that 
a(T) oc Vr and the specific heat C(T) - TM* oc ^/T. Measurements of the specific heat CiT) on CePdi_xRhx with 
X - 0.9 show a power-law temperature dependence. It is described by the expression C{T)/T - AT^'' with q ^ 0.5 and 
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Figure 35: Thermal expansion coef ficien t a(T) as a function of temperature in tlie interval 0.1 < T < 6 K. The experimental values for doping levels 
X = 0.90,0.87 are taken from Ref. |173|| . The solid lines represent approximations of the experimental values of a{T) = C[ V?, where C[ is a fitting 
parameter. 



A^const 1173 1. 

Figure [35] shows that at the critical point x = 0.90 at which the critical temperature of the ferromagnetic phase transi- 
tion vanishes, the thermal expansion coefficient is well approximated by the dependence a(T) oc ^/T as the temperature 
varies by almost two orders of magnitude. However, even a small deviation of the system from the critical point destroys 
the correspondence between this approximation and the experimental data. We note that it is possible to describe the 
critical behavior of two entirely different heavy-fermion metals (one is a paramagnet and the other a ferromagnet) by the 
function aiT) - ci Vt with only one fitting parameter ci. This fact vividly shows that fluctuations do not determine the 
behavior of a(T). Heat-capacity measurements for CePdi_xRhx with x - 0.90 have shown that C{T) oc y/T IJ73il . Thus, 
the electron systems of both metals can be interpreted as being highly correlated electron liquids. Hence, we conclude 
that the behavior of the effective mass given by Eq. (llOlb agrees with experimental facts. 

Measurements of a(T)/T on both CePdi_xRhx with x - 0.9 1 173] and CeNi2Ge2 1 14] are shown in Fig. |36] It is seen 
that the approximation a(T) — C3 Vt is in good agreement with the results of measurements of a{ T) in CePdi_xRhx and 
CeNi2Ge2 over two decades in T^- It is noted that measurements on Celn3_xSnx with x - 0.65 117511 demonstrate the same 
behavior a(T) oc yjT (not shown in Fig. l36T l. As a result, we suggest that Celn3_xSnx with x — .65, CePd i-xRhx with 
X ^ 0.9, and CeNi2Ge2 are located at FCQPT; recall that CePdi_xRhx is a three dimensional FM ifnlflTl . CeNi2Ge2 
exhibits a paramagnetic ground state 114] and Celn3_xSnx is AFM cubic metal il75ll . 

The normalized effective mass M*^(Tn) extracted from measurements on the HF metals YbRh2(Sio.95Geo.()5)2, 
CeRu2Si2, CePdi_xRhx and CeNi2Ge2 is reported in Fig. [37] Clearly, the scaling behavior of the effective mass given 
by Eq. (1106b is in accord with the experimental facts and M*^(Tn), shown by inverted triangles and collected on the 
AFM phase of YbRh2(Sia95Geo.o5)2 1 170], coincides with that collected on the PM phase (shown by upright triangles) of 
YbRh2(Sio.95Geo.o5)2 lll70ll . We note that in the case of LFL theory the corresponding normalized effective mass M^^ ^ 1 
is independent of both T and B as shown in Fig. |2] 

The peak temperatures Tmax, where the maxima of C(T)/T, XAciT) and a{T)/T occur, shift to higher values with 
increase of the applied magnetic field. In Fig. |38] T^^^{B) are shown for C/T and^AC, measured on YbRh2(Sio. 95Geo o5)2 . 
It is seen that both functions can be represented by straight lines intersecting at B ^ 0.03 T. This observationj 161 , 170ll 



as well as the measurements on CePdi_xRhx, CeNi2Ge2 and CeRu2Si2 demonstrate similar behavior 1114 192L 117411 which 
is well described by Eq. (1108b . 
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Figure 36: The normalized thermal expansion coefficient (a{TN)ITN)la(\) = MIj(Tn) for CeNi2Ge2 Ql] and for CePdi_xRhx with x = 0.90 fT?! 
versus Tn = T/Tm- Data obtained in measurements on CePdj-xRhx at B = are multiplied by some factor to adjust them at one point to the data for 
CeNi2Ge2. Dashed line is a fit to the data shown by the circles and pentagons at B = 0; it is represented by the function a{T) = cj, Vt with C3 being a 
fitting parameter The solid curve traces ffie universal behavior of the normalized eff'ective mass determined by Eq. )106t . see the caption to Fig. 1321 
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Figure 37: The universal behavior of M'j^(Tr^), extracted from xac(X,B)Ixac(Tm,B ) for both YbRh2(Sio.95Geo.o5)2 and CeRu2Si2 f9?, 'l70p, 
(C{T)IT)I(C(T„)IT„) for both YbRh2(Sio.95Geo.o5)2 and CePdi_xRhx with x = 0.80 (nlinSl. and (a{T)IT)l{a(T„)IT„) for CeNi2Ge2 uJ. 
All the measurements were performed under the application of magnetic field as shown in the insets. The solid curve gives the universal behavior of 
M*^ determined by Eq. )106t . see the caption to Fig. 1321 



We conclude, that subjecting the different experimental data (like C(T)/T,xac{T), a(T)/T etc) collected in measure- 
ments on different HF metals (YbRh2(Sio.95Geo.o5)2, CeRu2Si2, CePd i_xR hx, CeIn3_xSnx and CeNi2Ge2) to the above 
normalized form immediately reveals their universal scaling behavior \ VTf\. This is because all the above experimental 
quantities are indeed proportional to the normalized effective mass exhibiting the scaling behavior. Since the effective 
mass determines the thermodynamic properties, we further conclude that the above HF metals demonstrate the same scal- 
ing beha vior, independent of the details of HF metals such as their lattice structure, magnetic ground states, dimensionality 
etc il64l[l77ll . 
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10. Metals with a strongly correlated electron liquid 

For T <s: Tf, the function «()(p) given by Eq. (I2TI 1 determines the entropy S nfl{T) given by Eq. (|4|l of the HE liquid 
located above the quantum critical line shown in Eig. |6] Erom Eqs. ^ and dZTI l. the entropy contains a temperature- 
independent contribution, 

P f - Pi 

So ~ -^—^ ~ (130) 
Pf 

where r = (x - Xfc)/xfc- Another specific contribution is related to the spectrum e(p), which ensures a link between 
the dispersionless region (pf - pi) occupied by EC and the normal quasiparticles in the regions p < pi and p > pj. This 
spectrum has the form £;(p) oc (p - p/)^ ~ {pi - p)^. Such a shape of the spectrum, corroborated by exactly solvable 
models for systems with EC , leads to a contribut ion t o the heat capacity C ~ sjT/Tf ll4lll . Therefore, for < T ^ Tf, 
the entropy can be approximated by the function jlSlll 

SNFL{T)^So + aJ^+b^, (131) 

where a and b are constants. The third term on the right-hand side of Eq. (I131l l. which emerges because of the contribution 
of the temperature-independent part of the spectrum e(p), yields a relatively small addition to the entropy. As we will 
see shortly, the temperature-independent term So determines the universal transport and thermodynamic properties of 
the heavy-electron liquid with EC, which we call a strongly correlated Eermi liquid. The properties of this liquid differ 
dramatically from those of highly correlated Eermi liquid that at T — > becomes LEL Uquid. As a result, we can 
think of QCP of ECQPT as the phase transition that separates highly correlated and strongly correlated Eermi liquids. 
Because the highly correlated liquid behaves like LEL as T — > 0, QCP separates LEL from a strongly correlated Eermi 
liquid. On the other hand, as was shown in Subsection l4.4l at elevated temperatures the properties of both liquids become 
indistinguishable. Thus, as shall be seen below, both systems can be discriminated at diminishing temperatures when the 
impact of both QCP and the quantum critical line on the properties become more vivid. 

Eigure |39] shows the temperature dependence of S{T) calculated on the basis of the model functional (I102l i. The 
calculations were done with g = 7, 8, 12 and fi - be - 3. We recall that the critical value of g\s gc - 6.1161 . We see in 
Eig. [39]that in accord with Eq. (1130b So increases as the system moves away from QCP along the quantum critical line, see 
Eig. |6] Obviously, the term 5o on the right-hand side of Eq. ( 11311 ). which is temperature-independent, contributes nothing 
to the heat capacity; the second term in ( 11311 ) makes a contribution so that the heat capacity behaves as C{T) oc Vt. 
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Figure 39: Entropy S (T) as a function of temperature. The lines represent the approximation for S(T) based on Eq. )13U . the symbols mark the results 
of calculations based on )102t . 



10.1. Entropy, linear expansion, and Griineisen's law 

The unusual temperature dependence of the entropy of a strongly correlated electron liquid specified by Eq. dOlb 
determines the unusual behavior of the liquid. The existence of a temperature-independent term So can be illustrated by 
calculating the thermal expansion coefficient a{T) 1 181, 182], which is given by Ii20il 

-H'-^] ^ (132) 



3\ dT Ip 3V\ dP 



where P is the pressure and V is the volume. We note that the comp ressibility K - dfi/cKVx) does not develop a s ingularity 



at FCQPT and is approximately constant in systems with FC OlSOll . Substituting (I131I I in Eq. (I132l l. we find that lll8llll82 1 

aFc(T) flo ^FC 



T 



(133) 



where aq ~ dS^/dPis temperature-independent. In (1133b . we took only the leading contribution related to 5 o into account. 
We recall that 

55(r)^a /y (134) 
^ dT 2^Tf ' 

and obtain from Eqs. (1133b and (1134b that the Gruneisen ratio r(r) diverges as 

r(r).^.2^J5, (135) 

C[T) a y T 

from which we conclude that Griineisen's law does not hold in strongly correlated Fermi systems. 

Measurements that have been conducted with YbRh2(Sio.95Geo.o5)2 show that a/T cc l/T and that the Gruneisen ratio 
diverges as r(T) ^ J"^, q ^ 0.33, which allows classifying the electron system of this compound as strongly correlated 
liquid il4i] . Our estimate q - 0.5, which follows from Eq. (1135b . is in satisfactory agreement with this experimental 
value. The behavior of a(T)/T given by Eq. (1133b contradicts the LFL theory, according to which the thermal expansion 
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coefficient a{T)/T - M* - const as T — > 0. The l/T-dependence of the ration ajT predicted in Ref. il82ll is in good 
agreement with facts collected on YbRh2(Sio.95Geo.o5)2 [11411 . 

Equation (ISTT i implies that M*(T — > 0) ^ oo and that the strongly correlated electron system behaves as if it were 
placed at the quantum critical point. Actually, as we have seen in Subsection 14.41 the system is at the quantum critical 
Une xlxpc < 1, and critical behavior is observed for all x < Xfc as T — > 0. It was shown in Section |5] that as T — > 0, 
the strongly correlated electron liquid undergoes the first-order quantum phase transition, because the entropy becomes a 
discontinuous function of the temperature: at finite temperatures, the entropy is given by Eq. ( 1131b . while 5(r = 0) = 0. 
Hence, the entropy has a discontinuity 6S = 5o as T — > 0. This implies that, as a result of the first-order phase transition, 
all critical fluctuations are suppressed along the quantum critical curve and the respective divergences, e.g., the divergence 
of T(T), are determined by quasiparticles and not critical fluctuations, as could be expected in the case of an ordinary 
quantum phase transition |4]. We note that according to the well-known inequality 115611 q < TSS , in our case the heat q 
of the first order transition tends to zero as its critical temperature T^l — » 0. 



10.2. The T — B phase diagram o/ YbRh2Si2, Hall coefficient and magnetization 

To study the T - B phase diagram of strongly correlated electron liquid, we examine the case where NFL behavior 
emerges when the AF phase is suppressed by an external magnetic field B, as it is in the HE metals YbRh2(Sio.95Ge() 05)2 
and YbRh2Si2 ifllfll. 

The antiferromagnetic phase is LEL with the entropy vanishing as T 0. Eor magnetic fields higher than the critical 
value Bed at which the Neel temperature Tnl{B — > B^o) — > 0, the antiferromagnetic phase transforms into a weakly 
polarized paramagnetic strongly correlated electron liquid [Tsll . As shown in Section |7] a magnetic field applied 
to the system with 7 = splits the EC state occupying the interval (pf - pi) into Landau levels and suppresses the 
superconducting order parameter K{p). The new state is specified by a multiply connected Eermi sphere, on which a 
smooth quasiparticle distribution function «o(p) in the interval {pj - pi) is replaced with a distribution v(p) as seen from 
Fig- 19] Hence, the behavior of LEL is restored and characterized by quasiparticles with the effective mass M*(B) given 
by Eq. ( |73] ). When the temperature increases so high that T > T*(B) with T*{B) given by Eq. ( l76b , the entropy of the 
electron liquid is determined by Eq. ( I131l l. The described behavior of the system is shown in the T - B diagram in Eig. 
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Figure 40: The T - B phase diagram of a strongly correlated electron liquid. The line Tr^(B/Bco) represents the dependence of the Neel temperature on 
the field sfi'ength B. The black dot at T = T^n, marks the critical temperature at which the second-order AF phase transition becomes a first-order one. 
For T < Tcrit, the heavy solid line represents the function TfjiB/Bco), when the AF phase transition becomes a first-order one. The strongly correlated 
liquid in the NFL region is characterized by the entropy S nfl given by Eq. )13U . The line separating the strongly correlated liquid (NFL) from the 
weakly polarized electron liquid, which behaves like the Landau Fermi liquid, is described by the function T'(B/BcO - 1) oc sJB/BcO - 1 [see Eq. j76H . 
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In accordance with the experimental data, we assume that at relatively high temperatures, such that T/Tno ~ 1, where 
Tjsfo is the Neel temperature in a zero magnetic field, the antiferromagnetic phase transition is a second-order one 1 15]. In 
this case, the entropy and other thermodynamic functions at the transition temperature T^i are continuous. This means 
that the entropy Saf of the antiferromagnetic phase coincides with the entropy S nfl of the strongly correlated liquid 
given by Eq. (1131b : 

Saf(T ^ TNim = Snfl(T ^ TMdB)). (136) 

Since the antiferromagnetic phase behaves like LFL, with its entropy Saf{T — > 0) — > 0, Eq. ( I136l l cannot be satisfied 
at sufficiently low temperatures T < Tdt because of the temperature-in dependen t term Sq. Hence, the second order 



antiferromagnetic phase transition becomes the first order one at T = Tcrit II183LI18411 as shown by the arrow in Fig. |40] A 
detailed consideration of this item is given in SectionfT4l 

At r = 0, the critical magnetic field Bco in which the antiferromagnetic phase becomes LFL is determined by the 
condition that the ground-state energy of the antiferromagnetic phase be equal to the ground-state energy E[no{p)] of the 
HF liquid with FC, since, as it was shown in Subsection llO.il the heat of the transition ^ = 0. This means that the ground 
state of the antiferromagnetic phase is degenerate at B = Bco- Hence, at B ^ Bco the Neel temperature T^i tends to zero 
and the behavior of the effective mass M*(B > B^o) is determined by Eq. ( f73] l. so that M*(B) diverges as B — » Z?cO from 
top. As a result, at T = 0, the phase transition separating the antiferromagnetic phase existing at B < Bcq from LFL taking 
place at B > B^o is the first order quantum phase transition. The driving parameter of this phase transition is the magnetic 
field strength B. We note that the respective quantum and thermal critical fluctuations disappear atT < Tcrit because the 
first-order antiferromagnetic phase transition occurs at such temperatures. 

We now examine the jump in the Hall coefficient detected in measurements involving YbRhaSia |185]. The Hall 
coefficient Rh(B) as a function of B experiences a jump as T — > when the applied magnetic field reaches its critical 
value B - Bco, and the n be comes even higher than the critical value at B = Bco + SB, where 6B is an infinitely small 



magnetic field strength lIlSSll . As shown in Section|7] when T -Q, the application of the critical magnetic field B^o, which 
suppresses the antiferromagnetic phase with the Fermi momentum pf restores LFL with the Fermi momentum pj > pf. 
When B < Bco, the ground-state energy of the antiferromagnetic phase is lower than that of the LFL state induced by 
the application of magnetic field, but for B > B^o we are confronted with the opposite case, where the LFL state has the 
lower energy. At B = B^o and T - 0, both phases have the same ground state energy and T^l = 0, because the phases are 
degenerate, being separated by the first order phase transition as shown in Fig. |40] 

Thus, at r = and B = B^o, an infinitely small increase 6B in the magnetic field leads to a finite discontinuity in 
the Fermi momentum. This is because the distribution function becomes multiply connected (see Fig. |9]i and the number 
of mobile electrons does not change. Thus, the antiferromagnetic ground state can be viewed as having a "small" Fermi 
surface characterized by the Fermi momentum pf, correspondingly the paramagnetic ground state at B > B^o has a 
"large" Fermi surface with pf > pf. As a result, the Hall coefficient experiences a sharp jump because Rh{B) oc 1//?^ 
in the antiferromagnetic phase and Rh{B) VPy in the paramagnetic phase. Assuming that Rh{B) is a measure of the 



Fermi momentum 111 8511 (as is the case with a simply connected Fermi volume), we obtain 



Ml^If^^l^^'^l^^l.d^, (137) 
Rh{B = B,o + 6) Pf xfc 

where Sq/xfc is the entropy per heavy electron and disa constant li ~ 5. It follows from Eq. ( I137l i that the discontinuity 
in the Hall coefficient is determined by the anomalous behavior of the entropy, which can be attributed to So- Hence, 
the discontinuity tends to zero as r — > and disappears when the system is on the disordered side of FCQPT, where the 
entropy has no temperature-independent term il83ll . 

We now turn to the magnetization which is determined by Eq. (II 10b . For T « T*{B), the effective mass is given by 
Eq. (l73T l and the static magnetization is 

M(B) ^aM^JB- B,o. (138) 

Figure|4T]shows that the function M(B) determined by Eq. ( 11381 ) is in good agreement with the data obtained in measure- 
ments on YbRh2(Sio.95Geo.o5)2 f^']. We note that B^o - in this case. 

We examine the experimental T - B diagram of the heavy-fermion metal YbRh2Si2 HE] shown in Fig. |42l In the 



LFL state, the behavior of the metal is characterized by the effective mass M*(B), which diverges as 1/ - ^cO lll5ll . 
It is quite evident that Eq. ( f73T l provides a good description of this experimental fact: M*(B) diverges as B — > B^o at 
rjv(B - Bco) = and, as Fig. |4T] shows, the calculated behavior of the magnetization agrees with the experimental data. 
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Figure 41: The values of magnetization M{B) obtained in measurements involving YbRh2(Sio.95 Geo 05)2 (black squares) [Tj- The curve represents the 
field-dependent function MiB) = am Vb given by Eq. )138t . where om is a fitting parameter. 



The magnetic-field dependence of the coefficient A{B) shown in the left panel of Fig. [TT]is also in good agreement with 
experimental facts collected on YbRhaSia iflSl . Figure |42] shows that in accordance with ( 176] ), the curve separating the 
LFL region from the NFL region can be approximated by the function c yjB — Bg o with a f itting parameter c. Bearing 



in mind that the behavior of YbRhaSia is like that of YbRh2(Sio.95Ge()()5)2 l?, 14, 170l 186 1. we also conclude that the 
thermal expansion coefficient a{T) is temperature-independent and that the Griineisen ratio diverges as a function of T in 
the NFL state fTT). We conclude that the entropy in the NFL state is determined by Eq. ( I131l l. Since the antiferromagnetic 
phase transition is the second order at relatively high temperatures [ 15J, we can predict that as the temperature decreases, 
the phase transition becomes the first o rder. The above description of the behavior of the Hall coefficient Rh{B) also 
agrees with the experimental facts 111 8511 . 

Thus, we conclude that the T — B phase diagram of the strongly correlated electron liquid shown in Fig. l40l agrees 
with the experimental T - B diagram obtained from experiments involving the heavy-fermion metals YbRh2Si2 and 
YbRh2(Si().95Geo.o5)2 and shown in Fig. |42] 

10.3. Heavy-fermion metals in the immediate vicinity of QCP 

We now consider the case where 6pfc = iPf - pd/PF ^ 1 and the electron system of HF metal is in a state close 
to QCP while remaining on the ordered side that is at the quantum critical line, see Fig. |6] It follows from Eq. ( l78l l that 
when the system is placed in a magnetic field (B - Bco)/Bco > B^r, the system passes from the ordered side of FCQPT 
to the disordered side, or the strongly correlated liquid transforms into the highly correlated one. As a result, when 
T < T*(B), the effective mass M*{B) is determined by Eqs. ( |93] | and (|99] |; thus both the Kadowaki-Woods relation and the 
Wiedemann-Franz law remain valid, and there are quasiparticles in the system. The resistivity then behaves as described 
in Subsection |9.5l 

In magnetic field with B ^ Bco and at temperatures Tf » T > T*{B), the system behaves like the strongly correlated 
Fermi liquid, the effective mass M*(T) is given by Eq. dSTT i. and the entropy is determined by Eq. ( 11311 ). The thermal 
expansion coefficient a(T) is temperature-independent [as follows from Eq. (Il33b 1, and the Griineisen ratio diverges, as 
follows from Eq. ( 11351 ). It follows from Eq. ( |3TI ) that the width y(T) oc T (see also Section l53T l. Hence, at » T » 
T*{B), the temperature-dependent part of the resistivity behaves as Ap{T) oc y{T) oc r in either case, when the electron 
system is in the highly correlated state or in the strongly correlated state. 



65 



120- 




Figure 42: T - B phase diagram for YbRh2Si2; the symbols denote the experimental data It HT^ . The line Tjv depicts the field dependence of the Neel 
temperature Tnl(B). In the NFL region, the behavior of the strongly correlated liquid is characterized by the entropy Snfl determined by Eq. )13H . 
The line sepai'ating the NFL region from the LFL region is approximated by the function T*{B - B^o) = f ~ ^cO given by Eq. )76) where c is a fitting 
parameter. 



We assume that the system becomes superconducting at a certain temperature Tc- In contrast to the jump SCiJc) of 
the heat capacity at Tc in ordinary superconductors, which is a Hnear function of Tc, the value of SCiJc) is independent 
of Tc in our case. Equations ( l62b and ( |63] ) show that both SCiJc) and the ratio dC(Tc)/C„(Tc) can be very lar ge compare d 



to the corresponding quantities in the ordinary BCS case as it was observed in the HF metal CeCoIns 11791 llOlL 118711 



Experiments show that the electron system in CeCoIns can be considered as a strongly correlated electron liqui d. Indeed, 
for T > T*(B), the linear thermal expansion coefficient a(T) oc const and the Griineisen ratio diverges fl 18911 [see Eqs. 



( 1133b and ( Il35b 1. so we may assume that the entropy is given by ( 1131b . 

A finite magnetic field takes the system to the disordered side of FCQPT; for T < T*{B), the system behaves like 
the highly correlated liquid with the effective mass given by Eq. ( 1931 ). Estimates of 5pfc based on calculations of the 
magnetic susceptibility show that 6pfc - 0.044 [79]. We conclude that Bcr ~ 0.01, as follows from Eq. ( iTSl l, and the 
electron system of the heavy-fermion metal CeCoIns passes, in relatively weak magnetic fields, to the disordered side of 
FCQPT and acquires the behavior characteristic of highly correlated liquid. We note that the estimated value of dpfc 
provi des an explanation for the relatively large jump 5C{Tc) [79] observed at Tc - 2.3 K in experiments with CeCoIns 

lisl. 

As Fig. |43] shows, the behavior A(Z?) oc Bh{B) oc M*(B) oc (B - BcqY'^^^ specified by Eq. dl 13b is in good agreement 
with the experimental results [ 30, 3 1 ] . The coefficient Bh(B) determines the T^-dependence of the thermal resistance, and 



the ratio A{B)IBh{B) is field-independent, with A/B^ ^ 0.70 (lHIMl- In the LFL state, the Kadowaki- Woods relation 



and the Wiedemann-Franz law hold, and the behavior of the system is determined by quasiparticles SmEIi. Thus, 
we conclude that our description is in good agreement with the experimental facts. 

At low temperatures and in magnetic fields Bycd ~ Be, [see Eq. (fTSTll. the electron system is in its LFL state. As the 
temperature increases, the behavior of the strongly correlated liquid determined by the entropy S[) is restored at T*{B), 
and the effective mass becomes temperature-dependent, according to Eq. (ISTb . To calculate T*{B), we use the fact that 
the behavior of the effective mass is given by Eq. ( |93l ) for T < T*(B) and by Eq. (|3TI ) for T > T*(B). Since the effective 
mass cannot change atT - T*{B), we can estimate T*(B) by equating these two values of the effective mass. As a result, 
we obtain 

T*{B)^{B-Bcafl\ (139) 
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The function T*(B) ( |139l l is shown by the dotted Hne in Fig. |43] As the magnetic field becomes stronger, B » Bcr, the 
system becomes the highly correlated liquid in which the behavior of M*{T) is given by Eq. ( |99l ) and that of M*{E) by 
Eq. (|93T l. Comparison of these two types of behavior yields Eq. (llOOI i. The function T*(B) given be Eq. ( llOOl i is depicted 
by the light solid line in Fig. |43] Clearly, both lines match the experimental results. 




Figure 43: T ~ B phase diagram for CeCoIns . In the left panel, shown are A(B) and Bh(B) that determine the T^-dependence of the resistance and heat 
transfer in the LFL state induced by the magnetic field; the symbols mark the experimental data. The right panel depicts the curves of phase transitions 
in a magnetic field; the line separates the normal (NFL) state from the superconducting (SC) state 1 188]; the solid curve coiresponds to the second-order 
phase transitions, the dashed curve coiresponds to the first-order phase transitions, the black square (at Tq) is the point where second-order transitions 
become first-order transitions. The dotted line represents the function T'(B) calculated in accordance with U39t for the transition region between the 
LFL and NFL states. The light solid line represents the function T'{B) calculated according to Eq. jlOOt for the transition region (when B > B„) 
between the highly coiTelated and strongly coiTelated liquids; the black squares mark the experimental results obtained from resistivity measurements 

da El. 



Using Eq. (I136I I to study the superconducting phase transition, we can explain the main universal properties of the 
T - B phase diagram of the HF metal CeCoIns shown in Fig. |43] The latter substance is a cZ-wave superconductor 
with Tc - 2.3 K, while field tuned QCP with a critical field of Bco - 5.1 T coincides with Bc2, the upper critical field 
where superco nduct ivity vanishes Under the application of magnetic fields Bco, CeCoIns demonstrates the 



NFL behavior 1118911 . It also follows from the above consideration given in Subsection 19.91 that Bc2 > Bco- Therefore, 
the approximate equality Bc2 - SeO observed in CeCoIns is an accidental coincidence that has to disappear under the 
application of external factors. Indeed, Bc2 is determined by Aq which in turn is given by the coupling of electrons with 
magnetic, phonon, etc excitations rather than by B^o- As a result, under the application of pressure influencing diff'erently 
the couphng constant Aq and Bco, the above coincidence is lifted in complete agreement with experimental facts, so that 
Bc2 > Bco 03311 as has been shown in Subsection l9.9l At relatively high temperatures, the supercond ucting-nor mal phase 
transition in CeCoIns shown by the solid line in the right panel of Fig. |43]is of the second order 1 171 , 190ll so that S 



and the other thermodynamic quantities are continuous at the transition temperature Tc(B). Since Bc2 — Bco, upon the 
application of magnetic field, the HF metal transits to its NFL state down to lowest temperatures as it is seen from Fig. 
|43] As long as the phase transition is of the second order, the entropy of the superconducting phase SsciT) coincides with 
the entropy S nfl{T) of the NFL state and Eq. ( 1136b becomes 

SsciT ^ T,{B)) = Smfl(T ^ T,{B)). (140) 

Since Ssc(T ^ 0) ^ 0, Eq. (1140b cannot be satisfied at sufficiently low temp eratu r es du e to the presence of the 
temperature independent term So- Thus, in accordance with experimental results 1 171 , 190ll , the second order phase 



transition converts to the first order one below some temperature To(B) Il97ll . To estimate To(B), we use the scaling idea 
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of Volovik (see Ref. lll9lll for details), who derived the interpolation formula for the entropy of a li-wave superconductor 
in a magnetic field B, while S nfl has been estimated in |79']. As a result, we obtain To(B)/Tc ^ 0.3. This point coincides 
well with the experimental value, shown on the Fig. |43] Note that the prediction that the superconducting phase transition 
may change its order had been made in the early 1960-s [ 192]. Since our consideration is based on purely thermodynamic 
reasoning, it is robust and can be generalized to the cases when the superconducting phase is replaced by another ordered 
state, e.g. ferromagnetic state or antiferromagnetic one. 

Under constant entropy (adiabatic) conditions, there should be a temperature step as a magnetic field crosses the 
phase boundary due to the above thermodynamic inequality. Indeed, the entropy jump would release the heat, but since 
S - const the heat q is absorbed, causing the temperature to decrease in order to keep the constant entropy of the NFL 
state. Note that the minimal jump is given by the temperature-independent term So, and q can be quite large so that the 
corresponding HF metal can be used as an effective cooler at low temperatures. 



11. Scaling behavior of heavy fermion systems 

As we have seen in Section |2] and Subsection 19.11 the core of Landau Fermi liquid theory, the effective mass 
practically does not depend on temperature T, magnetic field B etc, M*^{T,E) - M*^ - const [19J. The thermodynamic 
functions such as the entropy S , heat capacity C, magnetic susceptibility x behave as in the case of noninteracting Fermi 
gas, namely low temperatures S jT oc CfT cc ^ cc M*^. In other words, when the inter-particle interaction is switching 
on and its strength A is increasing, a noninteracting Fermi gas continuously transforms into LFL with S(A), M*(/l) etc. 
becoming functions of A, while the main scaUng behavior of LFL, S oc M*^{A)T, remains untouched. This fact imposes 
strict conditions on the low temperature thermodynamic properties causing LFL exhibit the scaling behavior, which could 
be represented by some reference LFL with a normalized effective mass M]^^ = M*^{T, B)/M*^ ^ 1. As seen from Fig. |2] 
in the case of HF metals the scaling behavior of M]^ is different from that of M^. 

Here we show that despite of the very different microscopic nature of 2D "^He and HF metals with various ground 
state ma g neti c proper ties their NFL behavior is universal and can be captured well within the framework of FCQPT 



S 41 , 42 , 47 , 22 i 164ll that supports the extended quasiparticles paradigm. We concentrate on the NFL behavior observed 



when heavy fermion systems transit from their LFL to NFL states. This area is mostly puzzling and important because 
the behavior of the system in its transition state strongly depends on the scenario shaping the corresponding QCP. For 
example, if the transition region is described by theories based on quantum and thermal critical fluctuations there are no 
theoretical grounds to expect that these systems with different magnetic ground states could exhibit a universal scaling 
behavior 1 1, 3,15,11, 25J. 

There are many measurements of the heat capacity C{T, B), thermal expansion coefficient a{T, B) and the magnetic 
AC susceptibility B) on strongly correlated Fermi systems such as HF metals, high-Tc superconductors and 2D ^He 
carried out at different temperatures T, fixed magnetic fields B and the number density (or doping) x. Many of these 
measurements allow to explore the systems at their transition from the LFL state to the NFL one. Due to the equation 

C/r cc SIT oc ^ oc X oc alT ocM", (141) 

relating all the above quantities to the effective mass, these can be regarded as the effective mass M*(T, B, x) measurements 
producing information about the scaling behavior of the normalized effective mass Mj^. 

Experimental facts show that the effective mass extracted from numerous measurements on different strongly corre- 
lated Fermi systems upon using Eq. ( 11411 ) depends on magnetic field, temperature, number density and composition. As 
we have seen and shall see, a 4D function describing the normalized effective mass is reduced to a function of a single 
variable. Indeed, the normalized effective mass depends (as the effective mass does) on magnetic field, temperature, num- 
ber density and the composition of a strongly correlated Fermi system such as HF metals and 2D Fermi systems, and all 



these parameters can be merged into the single variable y by means of Eq. ( 1106b . 

77.7. Quantum criticality in 2D ^He 

We now discuss how the scaling be havior of the normalized effective mass MJ^ given by Eq. ( 11061 ) describes the 



quantum criticality observed in 2D ^He Il72lll52lll53l] . This quantum criticality is extremely significant as it allows us to 
check the possibility of the scaling behavior in the 2D system formed by ^He atoms which are essentially different from 
electrons. Namely, the neutral atoms of ^He are fermions interacting with each other by Van der Waals forces with strong 
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hardcore repulsion and a weakly attractive tail. The different character of the inter-particle interaction along with the fact 
that the mass of the ^He atom is 3 orders of magnitude larger than that of an electron, makes ^He systems have drastically 
different properties than those of HF metals. Because of this difference nobody can be sure that the macroscopic physical 
properties of these systems will be more or less similar to each other at their QCP. The 2D ^He has a very important 
feature: a ch ange in th e total density of -^He film drives it towards QCP at which the quasiparticle effective mass M* 



diverges ll72lll52lll53ll as seen from Figs. |4]and[T3] This peculiarity permits to plot the experimental dependence of 
the normalized effective mass versus temperature as a function of the number density x, which can be directly compared 
with M*i^ given by Eq. ( 1106b . As a result, 2D ''He, being an intrinsically isotropic Fermi-liquid with negligible spin- 
orbit interaction becomes an ideal system to test a theory describing the NFL behavior Note that the bulk liquid ^He is 
historically the first object to which the LFL theory had been applied [19]. One may speculate that at a sufficiently high 
pressure the liquid ^He would exhibits the NFL behavior. Unfortunately, the application of pressure causes 3D ^He to 
solidify. 

Let us consider HF liquid at T = characterized by the effective mass M* . As it was shown in Section [8] at QCP 
X - Xfc the effective mass diverges at T = and the system undergoes FCQPT. The leading term of this divergence reads 

M'jx) 02 X 

= ai -H , z = , (142) 

m I - z Xfc 

where m is the bare mass. Equation ( 11421 ) is valid in both 3D and 2D cases, while the values of the factors oi and 02 
depend on dimensionality and inter-particle interaction At jc > Xfc (or z > 1) FCQPT takes place. Here we confine 
ourselves to the case x < Xfc- It is seen from Eq. (11421) that FCQPT takes place in 2D ^He at elevated densities due to 
Van der Waals forces with strong hardcore repulsion. This strong hardcore repulsion makes the potential energy produce 
the main contribution to the ground state energy resulting in strong rearrangement of the single-particle spectrum and 
FCQPT. We recall that in the heavy electron liquid FCQPT occurs at diminishing densities due to Coulomb interaction. 

When the system approaches QCP, the dependence of quasiparticle effective mass on temperature and number density 
X is governed by Eq. ( l89b . It follows from Fig. [15] that the effective mass M*(T) as a function of T at fixed x reveals 
three different regimes at growing temperature. At the lowest temperatures we have the LFL state. The effective mass 
grows, reaching its maximum M*^(T, x) at some temperature Tm{x) and subsequently diminishing as T^^^^ as seen from 
Eq. ( |99] l. Moreover, the closer is the number density x to its threshold value Xc, the higher is the rate of the growth. The 
peak value M^^ grows also, but the maximum temperature Tm diminishes. Near the Tm temperature the last "traces" of 
the LFL state disappear, manifesting themselves in substantial growth of M*(x). The temperature region beginning near 
Tm{x) signifies the crossover between the LFL state with almost constant effective mass and the NFL behavior with the 
r^^^^ dependence. Thus the Tm point can be regarded as crossover between the LFL and NFL states or regimes. 

As we have seen, M*{T, x) in the T and x range can be well approximated by a simple universal interpolating function. 
The interpolation occurs between the LFL (M* cc T~) and NFL (M* oc T^^^^) states, thus describing the above crossover 
Substituting T by the dimensionless variable y = T/Tm, we obtain the desired expression ( 1106b . It is possible to calculate 
Tm as a function of z- Equation (1142b shows that M^j cc 1/(1 - z) and it follows from (|99] | that oc T^^^^. As a result, 
we obtain 140,] 

Tm^(1-z)"\ (143) 



Equation (1141b demonstrates that M*(T) can be measured in experiments on strongly correlated Fermi systems. Upon 
normalizing both M*(T) by its peak value at each x and the temperature by Tm, we see from Eq. ( 1106b that all the curves 
merge into a single one demonstrating a scaling behavior. 

In Fig. |44] we show the phase diagram of 2D ^He in the variables T and z (see Eq. ( 1142b ). For the sake of comparison, 
the plot of the effective mass versus z is shown by the dashed line. The part of the diagram where z < 1 corresponds to 
HF behavior and consists of LFL and NFL parts, separated by the line Tm(z) oc (1 - z)^^^- We note here that our exponent 



3/2 - 1.5 is exact as compared to that from Ref. 1115311 1.7 ± 0.1. The good agreement between the theoretical and 
experimental exponents supports our FCQPT description of the NFL behavior of both 2D ^He and HF metals; the former 
system is in great detail similar to the latter Th e regi me for z > 1 consists of a low-temperature LFL piece (shadowed 
region, beginning in the intervening phase z < 1 115311 ) and the NFL state at higher temperatures. The former LFL piece 
is related to the peculiarities of the substrate on which the 2D ''He film is placed. Namely, it is related to weak substrate 
heterogen eity (steps and edges on its surface) so that quasiparticles, being localized (pinned) on it, give rise to the LFL 



behavior 015311 . The competition between thermal and pinning energies returns the system back to NFL state and hence 

69 



3 
*^ 

m 

0) 
Q. 

E 

0) 



\ NFL ; 


NFL 


\ / 

lfl\ 





10 
M 

n 
E 

> 

u 

0) 

!t 

0) 

o 



n 

Q. 
3 

o 



z=x/x. 



Figure 44: The phase diagram of the 2D 'He system. The pail for z < 1 coiresponds to HF behavior divided into the LFL and NFL pails by the line 
Tm{z) (1 - z)''^, where Tm is the temperature at which the effective mass reaches its maximum. The exponent 3/2 =1.5 coming from Eq. )143t is 
in good agreement with the experimental value 1.7 ±0.1 1153]. The dependence M'{z) oc (1 - z)"' is shown by the dashed line. The regime for z > 1 
consists of the LFL piece (the shadowed region, beginning in the intervening phase z < 1 |l53!|, which is due to the substrate inhomogeneities, see text) 
and NFL regime at higher temperatures. 



restores the NFL behavior. Note, that the presence of the substrate can be considered as the main difference between 2D 
^He and HF metals. Namely, the latter metals do not have a substrate, the above LFL piece would be absent or very thin 
if some 3D disorder (like point defects, dislocations etc) is present in a HF metals. 
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Figure 45: The normalized effective mass M*^ as a function of the normalized temperature T ITm at densities shown in the left lower comer. The 
behavior M'^ is extracted from exper imental data for the entropy in 2D 'He )153||| and 3D HF compounds with different magnetic ground states such as 
CeRu2Si2 and CePdi_xRlix (Ulnil, fitted by the universal function l fT06l . 



In Fig. 03 we report the experimental values of the effective mass M* (z) o btained by the measurements on ^He 
monolayer 17211 . These measurements, in coincidence with those from Ref. lll53ll . show the divergence of the effective 
mass at jc = x^.. To show, that our FCQPT approach is able to describe the above data, we represent the fit of M*{z) 
by the fractional expression coming from Eq. (1142b and the reciprocal effective mass by the linear fit M/M*{z) k aiz. 
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We apply the universal dependence ( 11061 ) to fit the experimental data not o nly in 2D He but in 3D HF metals as well 



M*^(y) extracted from the entropy measurements on the He film 11531115411 at different densities x < Xc smaller then the 



critical point Xc = 9.9 + 0.1 nm 
the ferromagnet CePdo.iRho 
different magnetic fields 



is reported in Fig. 



In the same figure, the data extracted from the heat capacity of 
II174II and the AC magnetic susceptibility of the paramagnet CeRu2Si2 ll92ll are plotted for 
It is seen that the scaling behavior of the normalized effective mass given by Eq. ( 1106b is in 
accord with the experimental facts. All substances are located at QCP, where the system progressively disrupts its LFL 
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Figure 46: The dependence of M^{T/Tm) on T/Tm at densities shown in the left lower corner. The behavior M'^ is extracted from experimental data 
for C(T)/T in 2D ^He (tI and for the magnetization Mq in 2D ^He 1 153]. The solid curve shows the universal function, see the caption to Fig. [45] 

behavior at elevated temperatures. In that case the control parameter, which drives the system towards its QCP Xfc is 
represented merely by a number density x. It is seen that the be havio r of the effective mass M^(y), extracted from S{T)/T 
in 2D ^He (the entropy S(T) is reported in Fig. S8 A of Ref. lll54ll ) looks very much like that in 3D HF compounds as 
was shown in Sections |9] 

The attempt to fit the available experimental data for C(T)/T in ^He 17211 by the universal function M^(y) is reported 
in Fig. |46l Here, the data extracted from heat capacity C(T)/T for ^He monolayer and magnetization Mq for bilayer 
JJ4], are reported. It is seen that the normaUzed effective mass extracted from these th ermo dynamic quantities can 
be well described by Eq. ( 1106b . We note the qualitative similarity between the double layer il53[l and monolayer ll72ll of 
■'He seen from Fig. |46] 

On the left panel of Fig. |47l we show th e density dependence of Tmax, extracted from measurements of the mag- 
netization M(){T) on ^He bilayer 111531 Il5411 . The peak temperature is fitted by Eq. ( 1143b . On the same figure, 
we have also reported the maximal magnetization M^ax- It is seen that M^ax is well described by the expression 
A/max 0^ (5'/r)n,ax oc (1 - 2)"', see Eq. (1142b . The right panel of Fig. |47]reports the peak temperature T^ax and the 
maxima l ent r opy (5/ r)max versus the number density x. They are extracted from the measurements of S{T)/T on ^He 



bilayer 111531 Il54tl . The fact that both the left and right panels extracted from Mo(T) and S/T demonstrate the same 



behavior shows once more that th ere a re indeed the quasiparticles, which determine the thermodynamic behavior of 2D 
^He (and also 3D HF compounds lll6411 ') near the point of their effective mass divergence. 

As seen from Fig. |47] the amplitude and positions of the maxima of magnetization M()(r) and S{T)/T in 2D ^He 
follow well Eqs. ( 1142b and ( 11431 ). while Eq. ( 1106b describes the scaling behavior of the normalized thermodynamic 
functions. We recall that we can calculate only relative values of the effective mass, that is the normalized effective mass, 
since the real values of Tm and M*^ are determined by the specific properties of the system in question. Thus, with only 
two values defining both the real value, for example, of the entropy and the corresponding temperature, it is possible to 
calculate the thermodynamic or transport properties of HF metals or 2D ^He. We conclude that Eq. ( 11061 ) allows us to 
reduce a 4D function describing the normalized effective mass to a function of a single variable. Indeed, the normalized 
effective mass depends on magnetic field, temperature, number density and the composition of a strongly correlated Fermi 
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system such as HF metals and 2D Fermi systems, and as we have seen above, all these parameters can be merged into the 
single variable by means of Eq. ( IIO6I 1 [M)]. We note that the validity of Eq. (IIO6I 1 is confirmed by numerical calculations 
as described in Subsection l9.3.1l 

In conclusion of this Subsection, we have described the diverse experimental facts related to temperature and number 
density (2D number density) dependencies of different thermodynamic characteristics of 2D ''He by the single universal 
function of one argument. The above universal behavior is also inherent to HF metals with different magnetic ground state 
properties. The amplitude and positions of the maxima of the magnetization Mo{T) and S{T)/T in 2D ''He are also well 
described. We have shown that bringing the different experimental data collected on strongly correlated Fermi systems to 
the above form immediately reveals their universal scaling behavior. 



77.2. Kinks in the thermodynamic functions 

To illuminate kinks or energy scales observed in the thermodynamic functions measured on HF metals lll7ll and 
2D ^He, we present in Fig. |48] the normalized effective mass M*j^ extracted from the thermodynamic functions versus 
normalized temperature (the left panel) and the norma lized thermodynamic functions proportional to T^M'^^ (the right 
panel) as a function of the normalized temperature [.I66I1 . 

Mt(y) extracted from the entropy S(T)/T and magnetization M measurements on the -^He fikn lll53ll at different 



densities x is reported in the left panel of Fig. 



48] In the same panel, the data extracted from the heat capacity of the 



ferromagnet CePdo.2Rho.8 1 174], CeCoIns |193|] and the AC magnetic susceptibility of the paramagnet CeRu2Si2 1921 
are plotted for different magnetic fields. It is seen that the universal behavior of the normalized effective mass given by 
Eq. (11061 ) and shown by the solid curve is in accord with the experimental facts. It is seen that the behavior of M*^(y), 
extracted from S {T)/T and magnetization M of 2D ^He looks very much like that of 3D HF compounds. In the right panel 
of Fi g. |48] t he no rmaUzed data on C(y), S(y),w(y) and M - M(y)+_yx(y) extracted from data collected on CePdi_xRhx 
ifnl , ^He jTsi], CeRu2Si2 CeCoIns f]3J\ and YbRu2Si2 (l7\ respectively are presented. Note that in the case 
of YbRu2Si2, the variable y - (B - Bc())iibITm can be viewed as effective normalized temperature. We remark that in 
Subsection l9.4.2l we calculate M as a function of magnetic field. 

It is seen from the right panel of Fig. |48]that all the data exhibit the kink (shown by arrow) at y > 1 taking place 
as soon as the system enters the transition region from the LFL state to the NFL one. This region corresponds to the 
temperatures where the vertical arrow in Fig. [T5]a crosses the hatched area separating the LFL from NFL behavior. It is 
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Figure 48: Energy scales in HF metals and 2D ^He. The left panel. The normaUzed effective mass versus the normalized temperature y = T/Tm- 
The dependence M^(>') is extracted from measurements of S{T) /T an d magnetizatio n M on 2D ^He |153|| ). from AC susceptibility xC^) collected on 
CeRu2Si2 i93l and from C(T)/T collected on both CePdi_xRhx Il74il and CeCoIns Il93ll . The data are collected for different densities and magnetic 
fields shown in the left bottom comer The solid curve traces the universal behavior of the normalized effective mass determined by Eq. U06t . 
Parameters C[ and q are adjusted (or xn{Tn, B) aX B = 0.94 mT. The right panel. The normalized specific heat C(y) of CePd|_xRhx and CeColns at 
different magnetic fields B, normalized entropy 5 (v) of 'He at different number densities x, and the normalized yxiy) at B = 0.94 mT versus normaUzed 
temperature y are shown. The upright triangles depict the normalized "average" magnetization M = M + Bx collected on YbRu2Si2 I I'^il. The kink 
(shown by the arrow) in all the data is clearly seen in the transition region y > 1. The soHd curve represents yM^(y) with parameters C[ and C2 adjusted 
for the magnetic susceptibility of CeRu2Si2 at B = 0.94 mT. 



also seen that the low temperature LFL scale of the thermodynamic functions (as a function of y) is characterized by the 
fast growth, and the high temperature scale related to the NFL behavior is characterized by the slow growth. As a result. 



we can identify the energy scales near QCP, discovered in Ref. II17I1 : the thermodynamic characteristics exhibit the kinks 
(crossover points from the fast to slow growth at elevated temperatures) which separate the low temperature LFL scale 
and high temperature one related to the NFL state. 

11.3. Heavy-fermion metals at metamagnetic phase transitions 

A Fermi system can be driven to FCQPT when narrow bands situated close to the Fermi surface are formed by the 
application of a high critical magnetic field Th e emergence of such state is known as metamagnetism that occurs 



when this transformation comes abruptly at 11194 1 



Let us assume that the magnetic field B„, is similar to that of Bcq driving a HF metal to its magnetic field tuned QCP. 



In our simple model both BcO and B,„ are taken as parameters. To apply equation ( 1106b when the critical magnetic field 
is not zero, we have to replace Bhy {B - B„,). Acting as above, we can extract the normalized effective mass M*^(Tn) 
from data collected on HF metals at their metamagnetic QCP. In Fig. |49]the extracted normalized mass is displayed. 
M*^{Tn) is extracted from measurements of C/T collected on URui ^aRho. osSia, Ce Ru2Si2 and CeRuaSii.gGeo.i at their 
metamagnetic QCP with B„ ^ 35 T, B,,, ^ 7 T and B,„ ^ 1 .2 T respectively 1 195 , 196ll . As seen from Fig. |49] the effective 



mass M*^{Tyi) in different HF metals reveals the same form both in the high magnetic field and in low ones as soon as 
the corresponding bands become flat, that is the electronic system of HF metals is driven to FCQPT. This observation is 
extremely significant as it allows us to check the universal behavior in HF metals when these are under the application of 
essentially different magnetic fields. Namely, the magnitude of the applied field (B ~ 10 T) at the metamagnetic point is 
four orders of magnitude larger than that of the field applied to tune CeRu2Si2 to the LFL behavior (B ~ \ mT). Relatively 
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Figure 49: The normalized eifective mass as a function of magnetic field versus the normalized te mperature. M'^(Tn) is extracted from measurements 
of CjT collected on URui i)2Rho.08Si2, CeRu2Si2 and CeRu2Sii 8Geo.2 at different magnetic fields Il95lll9al shown in the right panel. The solid curve 
gives the universal behavior of given by Eq. UQ6t . see also the caption to Fig. 1321 



small values of M^(Tm) observed in URui 92Rho.o8Si2 and CeRu2Si2 at the hig h fie lds and small temperatures can be 
explained by taking into account that the narrow band is completely polarized 119511 . As a result, at low temperatures 
the summation over the spins "up" and "down" reduces to a single direction producing the coefficient 1 /2 in front of the 
normalized eff'ective mass. At high temperatures the summation is restored. As seen from Fig. |49] these observations are 
in accord with the experimental facts. 



12. Asymmetric conductivity in HF metals and high-r^ superconductors 

The main subjects of investigation in experiments on HF metals are the thermodynamic properties. Therefore, it seems 
reasonable to study the properties of HF liquids that are determined by the quasi narticle distribution function n(p, T) and 
not only by the density of states or by the behavior of the effective mass M* llil Il72 , 197 , 198 , 199ll . As we shall see in 
this Section, the FC solutions no(p) leads to the NFL behavior and violate the particle-hole symmetry inherent in LFL 
and generate dramatic changes in transport properties of HF metals, particularly, the differential conductivity becomes 
asymmetric. As was shown in Section |7] the LFL behavior is restored under the application of magnetic field. Thus, 
we expect that in magnetic fields the asymmetric part of the diff'erential conductivity is suppressed. Scanning tunnel 
microscopy and point-contact spectroscopy closely related to the Andreev reflection are sensitive to b oth the de nsity of 
states and the probability of the population of quasiparticle states determined by the function n(p, T) ll200i 201 1. Thus, 



scanning tunnel microscopy and point-contact spectroscopy are ideal tools for studying specific features of the NFL 
behavior of HF metals and high-Tc superconductors. 



72.7. Normal state 

The tunnel current / running through a point contact of two ordinary metals is proportional to the applied volt- 
age V and to the square of the absolute value of the quantum mechanical transition amplitude t times the difference 
A^i(0)A^2(0)(ni(p, T) - n2(p, T)) 1 199], where A^i(O) A^2(0) are the density of states of the respective metals and n2(p, T)) 
and n2ip, T) are respectively the distribution functions of the respective metals. On the other hand, in the semiclassi- 
cal approximation, the wave function that determines the amplitude t is proportional to (A^i(0)A^2(0))"'^^. Therefore, the 
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density of states drops out from the final result and the tunnel current becomes independent of Ni{0)N2(0). Because the 
distribution n(p, T — » 0) — » 6(pf - p) as T — > 0, where 6(pf - p) is the step function, it can be verified that the differential 
tunnel conductivity cr^(y) - dl/dV is a symmetric or even function of V in the Landau Fermi-liquid theory. Actually, the 
symmetry of cr^(y) is obeyed if there is the hole-quasiparticle symmetry (which is present in the LFL theory). Hence, the 
fact that (Td{V) is symmetric is obvious and is natural in the case of metal-metal contacts for ordinary metals that are in 
the normal or superconducting state. 




Figure 50: The single particle energy £(k, T) (a) and the distribution function «(k, T) (b) at finite temperatures as functions of the dimensionless variable 
k = p/pf. The arrows show temperature measured in T/Ef. At T = 0.0001 the vertical line shows the position of the Fermi level Ef at which 
n(k, T) = 0.5 as depicted by the horizontal line. At diminishing temperatures 7 — > 0, the single particle energy £(k, T) becomes more flat in the region 
iPf ~ Pi) the distribution function n(k, T) in this region becomes more asymmetrical with respect to the Fermi level Ef producing the particle-hole 
asymmetry related to the NFL behavior 

We study the tunnel current at low temperatures, which for ordinary metals is given by the expression lfl99ll20(il] 

liV) = 2|f|2 J [n(£ - y) - n(e)] de. (144) 

where we use the atomic system of units e - m - h - \ and normalize the transition amplitude to unity, |fp = 1. Since 
the temperatures are low, we can approximate the distribution function n(e) by the step function 9{p - s); Eq. (I144l i then 
yields I{V) - a\V, and hence the differential conductivity cr^(y) = dlfdV - ai - const is a symmetric function of the 
applied voltage V. 

To quantitatively examine the behavior of the asymmetric part of the conductivity a-^iy), we find the derivatives of 



both sides of Eq. (1144b with respect to V. The result is the following equation for cr^(y): 



r,, = 1 J nieiz) - y r)(l - n{s{z) - V, T))^dz, (145) 



dz 



In the integrand in Eq. (I145I I. we used the dimensionless momentum z = p/pf instead of s for the variable, because n is no 
longer a function of s in the case of a strongly correlated electron liquid; it depends on the momentum as shown in Figs. 
|5]and|50] Indeed, the variable s in the interval (pf - pi) is equal to p, and the quasiparticle distribution function varies 
within this interval. It is seen from Eg. (11451) that the violation of the particle-hole symmetry makes (TdiV) asymmetric as 
a function of the applied voltage V lf6l [l72lll97l[l98ll . 



The single particle energy e(k, T) shown in Fig. |50] (a) and the corresponding n(k, T) shown in the panel (b) evolve 
from the FC state characterized by «o(k, T - Q) determined by Eq. (ITTT i. It is seen from Fig. |50](a), that at elevated 
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temperatures the dispersion e(k, T) becomes more inclined since the effective mass M*{T) diminishes as seen from Eq. 
( ISTT i. At the Fermi level eip, T) = yu, then from Eq. ^ the distribution function n(p, T) - 1/2. The vertical line in 
Fig. |50] crossing the distribution function at the Fermi level illustrates the asymmetry of the distribution function with 
respect to the Fermi level at T = 0.0001. It is clearly seen that the EC state strongly violates the particle-hole symmetry at 
diminishing temperatures. As a result, at low temperatures the asymmetric part of the differential conductivity becomes 
larger Under the application of magnetic fields the system transits to the LFL state that strongly supports the particle hole 
symmetry. Therefore, the application of magnetic fields restoring the symmetry suppresses the asymmetric part of the 
differential conductivity. 




Voltage (mV) 



Figure 51: Differential conductivity <tj{V) measured in the case of point contacts Au/CeCoIn5. The curves (Tj{V) are displaced along the vertical axis 
by 0.05. The co nductivity is normalized to its value atV = -2 mV. The asymmetry becomes noticeable at T < 45 K and increases as the temperature 
decreases (2031. 



After performing fairly simple transformations in Eq. (I145l l, we find that the asymmetric part 

AcrAV) = icTdiV) - (rA-V))/2 
of the differential conductivity can be expressed as 

2 J [« 



(z, T) + a[l-n(z, T)]^ 
dn(z, T) 1 - 2n(z, T) 



dz [a«(z,r) + [l -«(z,r)]]2 



dz, (146) 



where a - exp(-V7r). 

Asymmetric tunnel conductivity can be observed in measurements involving metals whose electron system is lo- 
cated near FCQPT or behind it. Among such metals are \\\g\\-Tc superconductors and heavy-fermion metals, e.g., 
YbRh2(Sio.95Geo.o5)2, CeCoIns, YbCu5_xAlx or YbRh2Si2. The measurements must be conducted when the heavy- 
fermion metal is in the superconducting or normal state. If the metal is in its normal state, measurements of AcrdiV) 
can be done in a magnetic field B > Bco at temperatures T*(B) < T < Ty or in a zero magnetic field at temperatures 
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Figure 52: The asymmetric conductivity Actj(V) as a function of V/fi for three values of the temperature T/fi (normalized to ji). The inset shows the 
behavior of the asymmetric conductivity extracted from the data in Fig. \5l\ 

higher than the corresponding critical temperature when the electron system is in the paramagnetic state and its behavior 
is determined by the entropy Sq. 

Recent mea surements of the differential conductivity in CeCoIns carried out using by the point-contact spectroscopy 
technique i202ll have vividly revealed the asymmetry in the differential conductivity in the superconducting (Tc - 2.3 
K) and normal states. Figure |5T] shows the results of these measurements. Clearly, Acr,/(y) is nearly constant when the 
heavy-fermion metal is in the superconducting state, experiencing no substantial variation near Tc, see also Fig. l56lbelow. 
Then it monotonically decreases as the temperature increases a2Q2ll . 

Figure [52l s hows the results of calculations of the asymmetric part Acrrf(y) of the conductivity cTdiV) obtained from 
Eq. (I146I I 1I172I1 . In calculating the distribution function n(z, T), we used the functional ( I102l i (with the parameters P - 3 
and g - 8). In this case, (pf - pd/ Pf - O.l. Figure |52] also shows that the asymmetric part Acr^CV) of the conductivity is 
a linear function of V for small voltages. Consistent with the Fig. |50]showing that the asymmetry of n(k, T) diminishes at 
elevated temperatures, the asymmetric part decreases with increasing temperature, which agrees with the behavior of the 
experimental curves in the inset in Fig. |52] 

We now derive an estimate formula for analyzing the asymmetric part of the differential conductivity. It follows from 
Eq. ( 1146b that for small values of V, the asymmetric part behaves as Acr^(y) oc V. Here, it is appropriate to note that 
the asymmetric part of the tunnel conductivity is an odd function of V, and therefore Acrj(y) must change sign when V 
changes sign. The natural unit for measuring voltage is 2T, because this quantity determines the characteristic energy 
for FC, as shown by Eq. (l33T l. Actually, the asymmetric part must be proportional to the size (pf - pd/pp of the region 
occupied by FC: 

y Pf-Pi V So 

2T Pf 2T Xfc 

where So/xfc ~ ipf - pd/ Pf is the temperature-independent part of the entropy [see Eq. (Il30b 1 and c is a constant of the 
order of unity. For instance, calculations of c using the distribution function displayed in Fig. |50]yield c ~ I. From Eq. 
( 1147b we see that when y ^ 2r and FC occupies a sizable part of the Fermi volume, {pf - pi)! pf ^ 1, the asymmetric 
part becomes comparable to the differential tunnel conductivity Acrj{V) ~ Vd{V). 
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Figure 53: Characteristic temperature behavior of (a) symmetiic dV/dl'iV) and (b) asymmetric dV I dl'"^ (V) parts of dV/dI(V) for heterocontact 
YbCu3 5 All 5 - Cu at B = T and different temperatures shown by the arrows. The inset shows the bulk resistivity p(T) of YbCu3 5AI1 5 L2Q3.1 . 



12.1.1. Suppression of the asymmetrical differential resistance in YbCus-^Alx in magnetic fields 

Now consider the behavior of the asymmetric part of the differential conductivity Acrj(y) under the application of a 
magnetic field B. Obviously, the differential conductivity being a scalar should not to depend on the direction of current /. 
Thus, the non-zero value of Acrj(V) manifests the violation of the particle-hole symmetry on a macroscopic scale. As we 
have seen in SectionQand Subsection |9.1| at sufficiently low temperatures T < T*(B), the application of a magnetic field 
B > BcQ leads to restoration of the LFL beha vior eUmin ating the particle-hole asymmetry, and therefore the asymmetric 
part of the differential conductivity disappears lll72[ll97ll . This prediction is in accord with the experimental facts col lecte d 
in measurements on YbCu5_xAlx of the differential resistance dV/dI(V) under the application of magnetic fields ll203ll . 
Representing the differential resistance as the sum of its symmetrical dVldlHV) and the asymmetrical part dV/dF^V), 



dVldliV) = dV/dPiV) + dV/dr{V), 
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we obtain the equation 



dV/dl^'iV) 



(148) 



[dV/dP(V)]^' 

Deriving Eq. (1148b . we assume that dV/dr\V) » dV/dF^V). Figure |53] [203,1 shows the temperature evolution of (a) 
the symmetric dVldlHV) and (b) the asymmetric dV/dI"''(V) parts at zero applied magnetic field. Also for the case of a 
heterocontact, the behavior of the symmetric part does not show a decrease in piT), while the asymmetric part decreases 
at elevated temperatures |203]. It seen from Fig. |53]that the behavior of the asymmetric part of the differential resistance 
given by Eqs. (11471 ) and ( I148I I is in accord with the experimental facts. 
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Figure 54: Characteristic magnetic-field beliavior of the asymmetiic part dV/dI'"'(V) of the differential conductivity is shown versus magnetic fields 
displayed in the legends for heterocontacts with different ,v = 1.3, 1.5, and 1.75 at 1.5 K t203ll . 

It is seen from Fig. |54]|203] that increasing magnetic fields suppress the asymmetric part. Thus, the application of 
magnetic fields destroys the NFL behavior and recovers both the LFL state and the particle-hole symmetry. Correspond- 
ingly, we conclude that the particle-hole symmetry is macroscopically broken in the absence of applied magnetic fields, 
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while the apphcation of magnetic fields restores both the particle-hole symmetry and the LFL state. It is see n from Figs. 
I53]m dl54lthat the asymmetric part shows a linear behavior as function of the voltage below about 1 mV 120311 as predicted 



72.2. Superconducting state 

Tunnel conductivity may remain asymmetric as a high-Tc superconductor or a HP metal pass into the superconducting 
state from the normal state. The reason is that the function no(p) again determines the differential conductivity. As 
we saw in Section |5] «o(p) is not noticeably distorted by the pairing interaction, which is relatively weak compared to 
the Landau interaction, which forms the distribution function no(p). Hence, the asymmetric part of the conductivity 
remains practically unchanged for T < T^, which agrees with the results of experiments (see Fig. ISTT ). In calculating the 
conductivity using the results of measurements with a tunneling microscope, we must bear in mind that the density of 
states in the superconducting state 

NsiE) = N(s-^) ^ , (149) 

determines the conductivity, which is zero for E < |A|. Here, E is the quasiparticle energy given by Eq. (l38T l. and 
e - fi = Ve^-A^. Equation ( |149t implies that the tunnel conduct ivity may be asymmetric if the density of states in 
the normal state N(s) is asymmetric with respect to the Fermi level 120411 . as is the case with strongly correlated Fermi 



systems with EC. Our calculations of the density of states based on model functional ( 1102b with the same parameters as 
those used in calculating Acrj(y) shown in Fig. l52l corroborate this conclusion. 




Figure 55: Density of states N{^, T) as a function of f = (£ - fi)/fi, calculated for three values of the temperature T (normalized to //). 

Figure |55] shows the results of calculations of the density of states A^(^, T). Clearly, A^(^, T) is strongly asymmetric 
with respect to the Fermi level. If the system is in the superconducting state, the values of the normalized temperature 
given in the upper right corner of the diagram can be related to Ai. With Ai ^ 2Tc, we find that 2T/jj ^ Ai/jj. Because 
A^(^, T) is asymmetric, the first derivative dN(^, T)/d^ is finite at the Fermi level, and the function N(^, T) can be written 
as N{^, T) ^ qq + ai^ for small values of ^. The coefficient ao contributes nothing to the asymmetric part. Obviously, the 
value of Acrj(y) is determined by the coefficient ai cc M*{^ - 0). In turn, M*(^ - 0) is determined by Eq. dTO . As a 
result, Eq. (11491 ) yields 

V So 

AcraiV)-ci — ^, (150) 
|A| xfc 
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Figure 56: Temperature dependence of the asymmetric parts A(Trf(V) of tlie conductance spectra extracted from measurements on CeCoInj l202ll . The 
temperatures are boxed and shown by the arrow for T < 2.60 K, otherwise by numbers near the curves. 

because {pf- pdlpF - So/xpc, the energy E is replaced by the voltage V, and ^ - VV^ - A^. The entropy So here refers 
to the normal state of a heavy-fermion metal. 
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Figure 57: Spatial variation of the spectra of the differential tunnel conductivity measured in Bi2Sr2CaCu208+x- Lines 1 and 2 belong to regions in 
which the integrated local density of states is very low. Low differential conductivity and the absence of a gap are indications that we are dealing with 
an insulator. Line 3 corresponds to a large gap (65 meV) with mildly pronounced peaks. The integrated value of the local density of states for curve 3 is 
small, but is larger than that for Hnes 1 and 2. Line 4 con'esponds to a gap of about 40 meV, which is close to the avera ge va lue. Line 5 con'esponds to 
the maximum integrated local density of states and the smallest gap about of 25 meV, and has two sharp coherent peaks l205ll . 
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Actually, Eq. ( llSOI l coincides with Eq. (11471 ) if we use the fact that the characteristic energy of the superconducting 
state is determined by Eq. ( l42l i and is temperature-independent. In studies of the universal behavior of the asymmetric 
conductivity, Eq. ( 11501 1 has proved to be more convenient than (|149t . It follows from Eqs. ( I147l i and ( 11501 ) that measure- 
ments of the transport properties (the asymmetric part of the conductivity) allow the determination of the thermodynamic 
properties of the normal phase that are related to the entropy Sq. Equation ( 1150b clearly shows that the asymmetric part of 
the differential tunnel conductivity becomes comparable to the differential tunnel conductivity at V ~ 2|A| if EC occupies 
a substantial part of the Fermi volume, (pj - pi)/pF - 1- In the case of the cZ-wave symmetry of the gap, the right-hand 
side of Eq. ( 11501 ) must be averaged over the gap distribution A(0), where (f> is the angle. This simple procedure amounts 
to redefining the gap size or the constant ci. As a result, Eq. ( 1150b can also be applied when V < Ai, where Ai is the 
maximum size of the c/-wave gap lll97ll . For the Andreev reflection, where the current is finite for any small value of V, 
Eq. ( 11501 ) also holds for V < Ai in the case of the i-wave gap. 
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Figure 58: The asymmetric part Ao"rf(y) of the differential tunnel conductivity in the high-T^ superconductor Bi2Sr2CaCu208+x, extracted from the data 
in Fig. 1571 as a function of the voltage V (mV). The lines are numbered consistent with the numbers of the lines in Fig. 1571 

It is seen from Fig. |56]that the asymmetrical part Acrj(y) of the conductivity remains constant up to temperatures of 
about Tc and persists up to temperatures well above Tc- At small voltages the asymmetric part is a linear function of V 
and starts to diminish at T > Tc- It follows from Fig. |56]that the description of the asymmetric part given by Eqs. ( 1147b 
and ( 11501 ) coincides with the facts obtained in measurements on CeCoIns. 



Low-temperature measurements with tunneling microscopy and spectroscopy techniques were used in 1120511 to detect 
an inhomogeneity in the electron density distribution in Bi2Sr2CaCu208+x- This inhomogeneity manifests itself as spatial 
variations in the local density of states in the low-energy part of the spectrum and in the size of the superconducting gap. 
The inhomogeneity observed in the integrated local density of states is not caused by impurities but is inherent in the 
system. Observation facilitated relating the value of the integrated local density of states to the concentration x of local 
oxygen impurities. 

Spatial variations in the differential tunnel conductivity spectrum are shown in Fig. [ST] Clearly, the differential tunnel 
conductivity is highly asymmetric in the superconducting state of Bi2Sr2CaCu208+x- The differential tunnel conductivity 
shown in Fig. |52]may be interpreted as measured at different values of Ai(x) but at the same temperature, which allows 
studying the Acr,/(y) dependence on Ai(ji:). Figure |58] shows the asymmetric conductivity diagrams obtained from the 
data in Fig. |57] Clearly, for small values of V, Acrj(y) is a linear function of voltage consistent with ( 11501 ) and the slope 
of the respective straight lines Acr^(y) is inversely proportional to the gap size Ai . 

Figure |59] shows the variation in the asymmetric part of the conductivity Acrj(y) as the temperature increases. The 
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Figure 59: Temperature dependence of the asymmetric part Acrj(V) of the conductivity sp ectra obtained in measurements for 
YBa2Cu307_x/Lao.7Cao.3Mn03 by the contact spectroscopy method; the critical temperature r^. =! 30 K |206|1 . 



measurements were done on YBa2Cu307_x/Lao.7Cao.3Mn03 with Tc - 30 K f206'] . Clearly, at T < in the region of the 
linear dependence on V, the asymmetric part Acr^(y) of the conductivity depends only weakly on the temperature; such 
behavior agrees with (1150b . When T > Tc, the slope of the straight Une sections of the AcrdiV) diagrams decreases as the 
temperature increases; this behavior is described by Eq. (1147b . We conclude that the description of the universal behavior 
of AcTdiV) based on the FCQPT is in good agreement with the results of the experiments presented in Figs. I52ll53ll54l 
|56]|58] and|59]and is valid for both high-Tf superconductors and heavy-fermion metals. 



13. Violation of the Wiedemann-Franz law in HF metals 

As early as in 1853, German physicists Gustav Wiedemann and Rudolph Franz discovered the empirical law 
stating that for a metal at a constant temperature the ratio of its thermal conductivity k{T) to its electrical conductivity 
(t{T) is a constant, k(T)/o-{T) =const. Later on, the Danish physicist Ludvig Valentin Lorenz showed that the above ratio 
is proportional to the temperature T, K(T)/cr(T) = LT, the proportionality constant L is known as the Lorenz number. 
What is called Wiedemann-Franz (WF) law is indeed an independence of the Lorenz number L on temperature. However, 
it was firmly established that the WF law is obeyed both at room temperatures and for low ones (several Kelvins); at the 
intermediate temperatures L - L{T). 

Strictly speaking, the Lorenz number is temperature-independent only at low temperatures; its theoretical value 

L() = lim = = "T^ (151) 

r^o To-(T) 3 



(ks and e are Boltzmann constant and electron charge, respectively) had been calculated by Sommerfeld in 1927 11208 1 
in the model of noninteracting electrons, obeying Fermi-Dirac statistics. The same result is obtained in LFL theory and 
reflects merely the fact that both thermal and electrical conductivities of a metal are determined by Landau quasiparticles. 
Due to this fact, possible deviations from the WF law can be regarded as a signature of NFL behavior in a sample. 

Actually, Eq. ( 1151b is usually referred to a s the Wiedemann-Franz (WF) law. It was shown that at T = Eq. ( 1151b 
remains vaUdfor arbitrarily strong scattering [209], disorder l210ll and interactions 121 1(1 . This law holds for ordinary 



metals ilHEll l214l ElSH and does not hold for HF metals ll3lll216l. 1217,] CeNiSn and CeCoIns, the electron-doped 



material i218ll Pr^-rCerCuOd-v. and the underdoped compound i219tl YbBaoCu^O.,. In CeNiSn, the experimental value 
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of the reduced Lorenz number L{T)/Lo ~ 1.5 changes Httle at T < 1 K. This rules out the phonon contribution to the 
violation of the WF law. In the electron-doped compound Pr2-.iCe vCu04_y the departure of L(T) from Lq at T > 0.3 
K is also by more then unity and even larger |218] than that in CeNiSn. Other experimental tests of the WF law have 
been undertaken in the normal state of cuprate superconductors. The phase diagram of these compounds shows evolution 
from Mott insulator for undoped materials towards metallic Fermi liquid behav ior for overdo ped cases. Upward shift 
L/Lq ^ 2-3 was measured in underdoped cuprates at the lowest temperatures 1 218 , 219l 22011 . In strongly overdoped 
cuprates, the WF law was found to be obeyed perfectly |26]. 

The physical mechanism for the WF law violation is usually attributed to the NFL behavior like in Luttinger and 
Laughlin liquids il44 , 145 . 221, 222ll or in the case of a marginal Fermi liquid 1 761 . Yet another possibility for the LFL 
theory and the WF law (11511) violation occurs near QCPs where the effective mass M* of a quasiparticle diverges. This 
is because at the QCP the Fermi liquid spectrum with finite Fermi velocity Vf - pf/M* becomes meaningless as in 
this case vf — » 0. In a standard scenario of the QCP |il[il|2231 the divergence of the effective mass is attributed to 
the vanishing of the quasiparticle weight z in the single-particle states close to second-order phase transitions, implying 
that the quasiparticles disappear in this region. A conventional scenario of the WF law violation, associated with critical 
fluctuations in the vicinity of the second order phase transition, has recently been suggested in i224ll . However, it has 



been shown in several works Il6|,|23l|24|], that the standard scenario of the QCP is flawed so that to describe the deviations 



of L from the WF value Lq, we apply a scenario of the QCP, where the NFL peculiarities are due to FCQPT. That is 
related to the rearrangement of the single-particle spectrum of strongly correlated electron liquid with the conservation of 
quasiparticle picture within the extended paradigm. 

Therefore, to describe theoretically the violation of the WF law within the FCQPT formalism, it is sufficient to use the 
well-known LFL formulas for thermal and electrical conducti vitie s with the substitution of the modified single particle 
spectrum into them. Such theory has been advanced in Refs. 1 225 , 226ll . The authors showed that close to the QCP the 
Lorenz number Lqcp{T = 0) = 1.81 Lq, i.e. almos t two times larger than that from the LFL theory ( 11511 ). This result 
agrees well with the experimental values (2161,1218]. Furthermore, the depend ence L(T)/Lo has been calculated for two 
topologically distinct phases (see Section [TSl) - "iceberg" phase and FC phase ll225|. |2261. Theoretical calculations have 
shown that in both phases the largest departure from the WF law occurs near QCP 022511226 1. Deep in the "iceberg" phase 
we have the reentrance of the "classical" WF law in a sense that L - Lq while in the deep FC phase the Lorenz number 
is temperature independent at lo w temperatures , but its value is slightly larger than Lq. This is due to the particle-hole 
symmetry violation in FC phase 1 17211 1971 122711 . 

Rece ntly, t he anisotropy of the WF law violation near the QCP has been experimentally observed in the HF metal 
CeCoInj ll217ll . In that paper, the above HF compound has been studied experimentally in external magnetic fields, close 
to the critical value Hc2, suppressing the superconductivity. Under these conditions, the WF law was found to be violated. 
The violation is anisotropic and cannot be attributed to the standard scenario of quasiparticle collapse. At the same 
time, close to the QCP, sufficiently large external magnetic fields reveal the anisotropy of the electrical conductivities 
aik oc< ViVk > (v, are the components of the group velocity vector) and thermal conductivities Kik oc< £(p)v,v^ > of a 
substance. This is because the magnetic field does not affect the z-components of the group velocity v so that the QCP 
T-dependence of the transport coefficients holds, triggering the violation of the WF relation L,, - ct^^ITk-,, - n^kB/3e^. 
On the other hand, the magnetic field B alters substantially the electron motion in the perpendicular direction, yielding 
considerable increase of the x and y components of the group velocity so that the corresponding components L,t do not 
depart from their WF value. 

Therefore, the flattening of the single particle spectrum e(p) of strongly correlated electron systems considerably 
changes their transport properties, especially beyond the point of FCQPT due to breaking of the particle-hole sym metry . 
Also, in topologically different "iceberg" phases the WF law is also violated near its QCP. The results of theoretical 11225 1 
and experimental investigations demonstrate that the FCQPT scenario with further occurrence of both "iceberg" and FC 
phases give natural and universal explanation of the NFL changes of the transport properties of HF compounds and the 
WF law violation in particular. 



14. The impact of FCQPT on ordinary continuous phase transitions in HF metals 

The microscopic nature of quantum criticality determining the NFL behavior in strongly correlated fermion systems of 
different types is still unclear. Many puzzling and common experimental features of such seemingly different systems as 
two-dimensional (2D) electron systems and liquid ^He as well as 3D heavy-fermion metals and high-Tc superconductors 
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suggest that there is a hidden fundamental law of nature, which remains to be recognized. To reveal this hidden law "the 
projection" of microscopic properties of the above materials on their observable, macroscopic characteristics is needed. 
One such "projections" is the impact of the FCQPT phenomenon on the ordinary phase transitions in HF metals. As we 
have seen in Subsection [103] the main peculiarity here is the continuous magnetic field evolution of the superconductive 
phase transition from the second order to the first one Il7llfl9(il I192I1 . The same changing of the order is valid for 



magnetic phase transitions. 

Exciting measurements on YbRh2Si2 at antiferromagnetic (AF) phase transition revealed a sharp peak in low- 
temperature specific heat, which is characterized by the critical exponent a - 0.38 and therefore differs drastically from 



those of the conventional fluctuation theory of second order phase transitions 1122811 . where a ^ Q.I 1156] . The obtained 
large value of a casts doubts on the applicability of the conventional theory and sends a real challenge for theories de- 
scribing the second order phase transitions in HF metals |228j]jjgniting_sttong theoretical effort to explain the violation of 



the critical universality in terms of the tri critical point 1 229i, i230l 123 IL 123211 . 



The striking feature of FCQPT is that it has profound influence on thermodynamically driven second order phase 
transitions provided that these take place in the NFL region formed by FCQPT. As a result, the curve of second order 
phase transitions passes into a curve of the first order ones at the tricritical point leading to a violation of the critical 
universality of the fluctuation theory. For example, as we have seen in Subsection l9.9l the second order superconducting 
phase transition in CeCoIng changes to the first one in the NFL region. It is this feature that provides the key to resolve 
the challenge. 

14.1. T — B phase diagram for YbRh2Si2 versus one for CeCoInj 

In Fig. |60] the Tffi line represents temperature TML{B)ITm) versus field B/Bco in the schematic phase diagram for 
YbRh2Si2, with Tffo - Tnl{B = 0). There Tni{B) is the Neel temperature as a function of the magnetic field B. The 
solid and dashed curves indicate the boundary of the AF phase at B/B^o < 1 1 15]. For B/Bco > 1, the dash-dot Une marks 
the upper limit of the observed LFL behavior This dash-dot line coming from Eq. (|76] | separates the NFL state and the 
weakly polarized LFL state, and in that case is represented by 



T* / B 



where ai is a parameter. We note that Eq. (1152b is in good agreement with experimental facts Ill5ll . Thus, YbRh2Si2 
demonstrates two different LFL states, where the temperature-dependent electrical resistivity Ap follows the LFL behavior 
Ap oc r^, one being weakly AF ordered {B < B^q and T < Tnl{B)) and the other being the weakly polarized (B > Bco and 



T < T*{B)) II15I1 . At elevated temperatures and fixed magnetic field, during which the system moves along the vertical 
arrow shown in Fig. |60l the NFL state occurs which is separated from the AF phase by the curve T^^ of the phase 
transitions. Consistent with the experiment al fac ts we assume that at relatively high temperatures T/Tnl{B) ^ 1 the AF 



phase transition is of the second order L15L 122811 . In that case, the entropy and the other thermodynamic functions are 
continuous functions at the line of the phase transitions T^i shown in Fig. |60l This means that the entropy of the AF 
phase Saf(T) coincides with the entropy Snfl{T) of the NFL state. Since the AF phase demonstrates the LFL behavior, 
that is Saf{T ^ 0) — > 0, while SNfiiT) contains the temperature-independent term given by Eq. ( l27T i. Thus, in the 
NFL region formed by FCQPT Eq. ( I136l l cannot be satisfied at diminishing temperatures and the second order AF phase 
transition inevitably becomes the first order one at the tricritical point with T = as shown in Fig. |60l At T = 0, the 
heat of the transition ^ = as was shown in Subsection 110.11 thus the critical field Bco is determined by the condition 
that the ground state energy of the AF phase coincides with the ground state energy of the weakly polarized LFL, and the 
ground state of YbRh2Si2 becomes degenerate at B = Bco. Therefore, the Neel temperature Tml{B Beo) — > 0. That 
means that at 7" = the system moving along the horizontal arrow shown in Fig. |60] transits to its paramagnetic state 
when the applied magnetic field reaches its critical value B = Bcq, and becomes even higher B = Bco + 6B, whe re dB is an 



infinitesimal magnetic field increment, while the Hall coefficient experiences the jump as seen from Eq. ( 11361) 123311 . 

Upon comparing the phase diagram of YbRh2Si2 depicted in Fig. |60]with that of CeCoIng shown in Fig. [3T] it is 
possible to conclude that they are similar in many respects. Indeed, the line of the second order superconducting phase 
transitions changes to the line of the first ones at the tricritical point shown by the square in Fig. [31] This transition takes 
place under the application of magnetic fields B > Bc2 > Bco (see Subsections 19. 9l and l 10. 31 ). where Bc2 is the critical field 
destroying the superconducting state, and Bcq is the critical field at which the magnetic field induced QCP takes place 
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Figure 60: Schematic T - B phase diagram for YbRh2Si2. The solid and dashed T^l curves separate the AF and NFL states representing the field 
dependence of the Neel temperature. The black dot at T = T^- shown by the arrow in the dashed curve is the tiicritical point, at which the curve 
of second order AF phase transitions shown by the solid line passes into the curve of the first ones. At T < T^-, the dashed line represents the field 
dependence of the Neel temperature when the AF phase transition is of the first order. The NFL state is characterized by the entropy S o given by Eq. 
1130) . The dash-dot line separating the NFL state and the weakly polarized LFL state is represented by T* given by Eq. 1152) . The horizontal solid 
arrow represents the direction along which the system transits from the NFL behavior to the LFL one at elevated magnetic field and fixed temperature. 
The vertical solid arrow represents the direction along which the system transits from the LFL behavior to the NFL one at elevated temperature and 
fixed magnetic field. The hatched circles depict the transition temperature T' from the NFL to LFL behavior 



We note that the superconducting boundary hne BdiT) at lower temperatures acquires the tricritical point due 



to Eq. (I136|l tha t cannot be satisfied at diminishing temperatures T < T^r, i e. the corresponding phase transition becomes 



first order 1117111 . This permits us to conclude that at lower temperatures, in the NFL region formed by FCQPT the curve 



of any second order phase transition passes into the curve of the first order one at the tricritical point. 

74.2. The tricritical point in the B — T phase diagram o/YbRh2Si2 

The Landau theory of the second order phase transitions is applicable as the tricritical point is approached, T ^ Tcr, 
since the fluctuation theory can lead only to further logarithmic corrections to the values of the critical indices. Moreover, 
near the tricritical point, the difference TffiiB) - Tc- is a second order small quantity when entering the term defining the 
divergence of the specific heat [15^. As a result, upon using the Landau theory we obtain that the Sommerfeld coefficient 
70 - C/T varies as yo " where f = T/Tni{B) with the exponent being a ^ 0.5 as the tricritical point is approached 

at fixed magnetic field lll56ll . We will see that a - 0.5 gives a good description of the facts collected in measurements of 
the specific heat on Yb Rh2Si 2. Taking into account that the specific heat increases in going from the symmetrical to the 



asymmetrical AF phase 015611 . we obtain 

ro(f)=Ai + -^. (153) 
V|f- i| 

Here, Bi = B-t are the proportionality factors which are different for the two sides of the phase transition. The parameters 
Ai = A±, related to the corresponding specific heat {C/T)±, are also different for the two sides, and "+" stands for f > 1, 
while "-" stands for f < 1 . 

The attempt to fi t the available experimental data for yo = C{T)/T in YbRh2Si2 at the AF phase transition in zero 



magnetic fields 1122 811 by the function (I153I I is reported in Fig. |6T| We show there the normalized Sommerfeld coefficient 



7oM+ as a function of the normalized temperature t = T /Tnq. It is seen that the normalized Sommerfeld coefficient jq/A-^ 



extracted from C/T measurements on YbRh2Si2 1122811 is well described in the entire temperature range around the AF 



phase transition by the formula ( I153l l with A+ - 1. 
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Figure 61: The normalized Sommerfeld coefficient yolA-^ as a function of the normalized temperature t = T/Tm given by the formula 1153) is shown 
by th e soHd curve. The normalized Sommerfeld coefficient is extracted from the facts obtained in measurements on YbRh2Si2 at the AF phase transition 
|228|1 and shown by the triangles. 



Now transform Eq. (11531) to the form 



It follows from Eq. ( 1154b that the ratios (yo 



1 



yo(0-Ai ^ 

Bi ^/\^^\^ 
Ai)IBi for f < 1 and f > 1 versus 



(154) 



1 - f| collapse into a single line on 
logarithmicxlogarithmic plot. The extracted from experimental facts 112281] ratios are depicted in Fig. |62] the coefficients 
A\ and B\ are taken from fitting yo shown in Fig. [61] It is seen from Fig. |62]that the ratios (yo - A.\)IBi shown by the 
upward and downwards triangles for f < 1 and f > 1, respectively, do collapse into the single line shown by the solid 
straight line. 

A few remarks are in order here. The good fitting shown in Figs. |6T]and|62]of the experimental data by the functions 
153 1 and ( 11541 ) with the critical exponent a = 1 /2 allows us to conclude that the specific-heat measurements on YbRh2Si2 



1 22811 are taken near the tricritical point and to predict that the second order AF phase transition in YbRh2Si2 changes to 
the first order under the application of magnetic field as it is shown by the arrow in Fig. |60] 1 184]. It is seen from Fig. 
l6n that at f ^ 1 the peak is s harp, while one would expect that anomalies in the specific heat associated with the onset 
of magnetic order are broad 112281 l234l 123511 . Such a behavior represents fingerprints that the phase transition is to be 
changed to the first order one at the tricritical point, as it is shown in Fig. |60] As seen from Fig. [61] the Sommerfeld 
coefficient is larger below the phase transition than above it. This fact is in accord with the Landau theory that states that 
the specific heat increases when passing from f > 1 to f < 1 il56ll . 



14.3. Entropy in YbRh2Si2 at low temperatures 

It is instructive to analyze the evolution of magnetic entropy in YbRh2Si2 at low temperatures. We start with consid- 
ering the derivative of magnetic entropy dS(B, T)/dB as a function of magnetic field B at fixed temperature Tf when the 
system transits from the NFL behavior to the LFL one as shown by the horizontal solid arrow in Fig. |60] Such a behavior 
is of great importance since exciting experimental facts [ 161] on measurements of the magnetic entropy in YbRh2Si2 al- 
low us to analyze the reliability of the theory employed and to study the scaling behavior of the entropy when the system 
is in its NFL, transition and LFL states, correspondingly. 

The quantitative analysis of the scaling behavior of dS{B, T)/dB is given in Subsection 19.4.41 Fig. |2T]reports the 
normalized {dS /dB)^ as a function of the normalized magnetic field. It is seen from Fig. |2T|that our calculations shown 
by the solid line are in good agreement with the measurements and the scaled functions (AM/AT)^ extracted from the 
experimental facts show the scaling behavior in a wide range variation of the normalized magnetic field B/Bm- The other 
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Figure 62: The temperature dependence of the ratios (yo - A[)/Bi for ; < 1 and r > 1 versus |1 - ^| given by the formula )154t is shown by the solid 
line. The ratios are extracted from the facts obtained in measurements of yo on YbRh2Si2 at the AF phase transition I228fl and depicted by the triangles 
as shown in the legend. 



thermodynamic and transport properties of YbRh2Si2 analyzed in Subsection 19.41 are also in good agreement with the 
measurements. These developments make our analysis of the AF phase transition quite substantial. 

Now we are in a position to evaluate the entropy S at temperatures T < T* in YbRh2Si2. AtT < T* the system in its 
LFL state, the effective mass is independent of T, is a function of the magnetic field B. As a result, Eq. (l73T i reads 



= fl2A/— -1' (155) 

where 02 is a parameter In the LFL state at T < T* when the system mo ves a long the vertical arrow shown in Fig. 
l60l the entropy is given by the well-known relation, S - M*Tn^ I p^^ - jqT il56ll . Taking into account Eqs. (1152b and 
( II55I 1 we obtain that atT^T* the entropy is independent of both the magnetic field and temperature, S(T*) ^ joT* - 
So - aimTNiTT^ /(a2p}-). Upon using the data |15], we obtain that for fields applied along the hard magnetic direction 
SoiBcoWc) ~ 0.037? In2, and for fields applied along the easy magnetic direction 5 o(ficO -L c) ~ 0.0057? In2. Thus, as it 
follows from Fig. |42]and in accordance with the data collected on YbRh2Si2 [15] we conclude that the entropy contains 
the temperature-independent part So [6^29^ which gives rise to the tricritical point. 

To summarize this Section, we remark that a theory is an important tool in understanding what we observe; we have 
demonstrated that the obtained value of a is in good agreement with the specific-heat measurements on YbRh2Si2 and 
conclude that the critical universality of the fluctuation theory is violated at the AF phase transition since the second 
order phase transition is about to change to the first order one, making a — > 1/2. We have also shown that in the NFL 
region formed by FCQPT the curve of any second order phase transitions passes into a curve of the first order ones at the 
tricritical point leading to the violation of the critical universality of the fluctuation theory. This change is generated by 
the temperature-independent entropy formed behind FCQPT. 



15. Topological phase transitions related to FCQPT 

We have now investigated the structure of the Fermi surface beyond QCP within the extended quasiparticle paradigm. 
We have shown that at T = there is a scenario that entails the formation of FC, manifested by the emergence of a 
completely flat portion of the single-particle spectrum. 
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Figure 63: Schematic plot of the single particle spectrum e(p) (a) and occupation numbers n(p) (b), con'esponding to LFL (curves 1), FC (curves 2) 
and iceberg (curves 3) phases at T = 0. For LFL the equation, e(p) = /j, has a single root equal to Fermi momentum pf. For iceberg phase, the above 
equation has countable set of the roots pi ...p^..., for FC phase the roots occupy the whole segment (pf - p, ). We note that p, < pf < pf and the states, 
where e(p) < /i are occupied (n=l), while those with e(p) > /i are empty (n=0). 



In this Section we consider different kinds of instabilities of normal Fermi liquids relative to several perturbations 
of initial quasiparticle spectrum e( p) and occupat ion numbers n(p) associated with the emergence of a multi-connected 
Fermi surface, see e.g. 1 23, 94, 95, 135 . 137 . 23 d]. Depending on the parameters and analytical properties of the Landau 
amplitude, such instabilities lead to several possible types of restructuring of initial Fermi liquid ground state. This 
restructuring generates topologically distinct phases. One of them is the FC discussed above, another one belongs to a 
class of topological transitions (TT) and will be called "iceberg" phase, where the sequence of rectangles ("icebergs") 
n(p) = and n{p) = 1 is realized at T = 0. 

In such considerations, we analyze stability of a fermion system with model repulsive Landau amplitude allowing 
us to carry out an analytical consideration of the emergence of a multi-connected Fermi surface |94n95il. We show, in 
particular, that the Landau amplitude given by the screened Coulomb law does not generate FC phase, but rather iceberg 
TT phase. For this model, we plot a phase diagram in the variables "screening parameter - coupling constant" displaying 
two kinds of TT: a 5/2-kind similar to the known Lifshitz transitions in metals, and a 2-kind characteristic for a uniform 
strongly interacting system. 

The common ground state of isotropic LFL with density is described at zero temperature by the stepwise Fermi 
function nf{p) - 9{pf - p), dropping discontinuously from 1 to at the Fermi momentum pf. The LFL theory states 
that the quasiparticle distribution function nip) and its single particle spectrum eip) are in all but name similar to those 
of an ideal Fermi gas with the substitution of real fermion mass m by the effective one M* list]. These ripip) and 
e(p) can become unstable under several circumstances. The best known example is Cooper pairing at arbitrarily weak 
attractive interaction with subsequent formation of the pair condensate and gapped quasiparticle spectrum fj^. However, 
a sufficiently strong repulsive Landau amplitude can also generate non-trivial ground states. The first example of such 
restructuring for a Fermi system with model repulsive interaction is FC t41i] . It reveals the existence of a critical value a„ 
of the interaction constant a such that at a = the stability criterion s{p) - (s(p) - Ef)/(p^ - p^) > fails at the Fermi 
surface s{pf) - (pf -instability). We recall that in the case of this instability the single particle spectrum s(p) possesses 
the inflection point at the Fermi surface, see Subsection l9.3.1l Then at a > an exact solution of a variational equation 
for n{p) (following from the Landau functional E(n{p))) exists, exhibiting some finite interval {pj - pd around pf where 
the distribution function n{p) varies continuously taking intermediate values between 1 and 0, while the single-particle 
excitation spectrum s{p) has a flat plateau. Equation (1211 1 means actually that the roots of the equation s(p) = p form an 
uncountable set in the range < p < pj, see Fig. |63] It is seen from Eq. (ISTT l that the occupation numbers n{p) become 
variational parameters, deviating from the Fermi step function to minimize the energy E. 
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The other type of phase transition the so-called iceberg phase occurs when the equation s{p) = ij has discrete countable 
number of roots, either finite or infinite. This is reported on Fig. |63]and related to the situation when the Fermi surface 
becomes multi-connected. Note that t he idea of mul ti-connected Fermi surface, with the production of new, interior 
segments, had already been considered lfl3lfl36[fl37ll . 

Let us take the Landau functional E{n{p)) of the form 



(2nf 



+ ^jjn(p)Ui\p-T?'Mp')^^, (156) 
which, by virtue of Eq. (O, leads to derivation of the quasiparticle dispersion law: 



The angular integration with subsequent change to the dimensionless variables x = p/pp, y - y{x) = 2Ms{p)/ pj^, 
z - 27T^ME/p^p, leads to simplification of the equations ( 11561 ) and ( 1157b 

z[v(x)] = J^^*^ ^x^V(x)]v(x)dx, (158) 
y(x) = x^ + V(x), (159) 

where 

1 r 

y(x) = — I x'v(x')u(x,x')dx' , 
X J 

x+.x' 

M r 

u{x,x)-—^ — u{t)tdt. (160) 
J^-^Pf J 

\x-x'\ 

Here u(x) = U(pfx) and the distribution function v(x) = n(pfx) is positive, obeys the normaUzation condition 



xMx)dx=l/3, (161) 



and the Pauli principle limitation v(x) < 1. The latter can be lifted using, e.g., the ansatz: v(x) -[1 + tanh77(jic)]/2. In the 
latter case the system ground state gives a minimum to the functional 

/Wx)] = J[l.tanh,(.)]{x^-,x^ 

+ x'[l + tanhri(x')]u(x, x')dx'}dx, (162) 

containing a Lagrange multiplier ju, with respect to an arbitrary variation of the auxiliary function ri(x). This allows us to 
represent the necessary condition of extremum 6f -0 in the form 

xMx)[l-y(x)][y(x)-fi]^0. (163) 



This means that either v(x) takes only the values and 1 or the dispersion law is flat: y(x) - 114 111 , in accordance with 
Eq. (|2TJ. The former possibility corresponds to iceberg phase, while the latter to EC. As it is seen from Eq. ( 1211 1. the 
spectrum E{p) in this case cannot be an analytic function of complex p in any open domain, containing the FC interval 
{pf — Pi). In fact, all the derivatives of e{p) with respect to p along the strip (pf - pi) should be zero, while this is not the 
case outside {pf - pi). For instance, in the FC model with U{p) - Uo/p 14111 the kernel, Eq. ( 1160b . is non-analytic 

u(x, x ) — — — -{x + x +\x — x'l), (164) 
n-^Pf 
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Figure 64: Occupation function for a multiconnected distribution. 



which eventually causes non-analyticity of the potential V{x). It follows from Eq. (|159l l, that the single particle spectrum 
is an analytic function on the whole real axis if V{x) is such a function. In this case FC is forbidden and the only 
alternative to the Fermi ground state (if the stability criterion gets broken) is iceberg phase corresponding to TT between 
the topologically unequal states with v{x) = 0, 1 [73]. 

On the other hand, applying the technique of Poincare mapping, it is possible to analyze the sequence of iterative maps 
generated by Eq. ( fTST l for the single-particle spectrum at zero temperature [23]. If the sequence of maps converges, the 
multi-connected Fermi surface is formed. If it fails to converge, the Fermi surface swells into a volume that provides a 
measure of entropy associated with the formation of an exceptional state of the system characterized by partial occupation 
of single-particle states and dispersion of their spectrum proportional to temperature as seen from Eq. (ISTT i. 

Generally, all such states related to the for mation of iceberg phases are classified by the indices of connectedness 



(known as Betti numbers in algebraic topology Il237ll238ll ') for the support of v(x). In fact, for an isotropic system, these 
numbers simply count the separate (concentric) segments of the Fermi surface. Then the system ground state corresponds 
to the following multi-connected distribution shown in Fig. l64l 

n 

v(x) = J] 0(^ - X2i-i)eix2i - x), (165) 

;=1 

where the parameters < xi < X2 < . . . < JCa,, obey the normalization condition 

n 

X^4-4-i) = l- (166) 
1=1 



The function z, Eq. ( 1158b . 

^"9^ x^lx^+y(x)]dx, (167) 

has the absolute minimum with respect to xi, . . . , X2„-i and to n > 1. To obtain the necessary condition of extremum, we 
use the relations 

l<;t<2«-l, (168) 

OXk \X2nj 
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following from Eq. ( 1166b and the dependence of the potential V(x) in the dispersion law y(x) on the parameters 

X], . . .,X2„-l 



V(x) 



1 " r 

~Z1 J ^'^(X'X'^dx'. (169) 



Subsequent differentiation of Eq. ( I167l i with respect to the parameters xi, . . . , X2n-i and the use of Eqs. ( II68I 1 and ( I169l l 
yield the necessary conditions of extremum in the following form 



(-I)\^[yfe)-3;(X2„)] = 0, 1<A:<2«-1. (170) 



This means that a multi-connected ground state is controlled by the evident rule of unique Fermi level yix^) = y(x2n) for 
all 1 < k < 2n - I (except for xi - 0). In principle, given the dispersion law y(x) all the 2« - 1 unknown parameters Xk 
can be found from Eq. ( 11701 ). Then, the sufficient stability conditions d^z/dxidxj = jiSij, y, > generate the generalized 
stability criterion. Namely, the dimensionless function 

a-{x) = 2Ms{p) = YZ~~2 ' 

X x^^ 

should be positive within filled and negative within empty intervals, turning to zero at their boundaries in accordance 
with Eq. ( 1170b . It can be proved rigorously that, for given analytic kernel u{x, x'), Eq. (1171b uniquely defines the system 
ground state. 

In what follows we shall label each multi-connected state, Eq. (1165b . by an entire number related to the binary 
sequence of empty and filled intervals read from X2„ to 0. Thus, the Fermi state with a single filled interval {x2 - 1 , xi =0) 
reads as unity, the state with a void at the origin (filled [x2, x\\ and empty [xi, 0]) reads as (10) = 2, the state with a single 
gap: (101) - 3, etc. Note that all even phases have a void at the origin and odd phases have not. 

For free fermions V{x) - 0, ^(x) - x^, Eq. (1170b only yields the trivial solution corresponding to the Fermi state 1. To 
obtain non-trivial reaUzations of TT, we choose U{p) to correspond to the common screened Coulomb potential: 

p^ + pI 

The related explicit form of the kernel, 

(x + x')^ + Xq 

M(x,x) = aln- — -, (173) 

(X - X'Y + Xq 

with the dimensionless screening parameter xq - polpp and the coupling constant a - 2Me^ jnpf, evidently displays the 
necessary analytical properties for existence of iceberg phase. Equations ( 1169b and ( 1173b permit to express the potential 
y(x) in elementary functions |94]. Then, the straightforward analysis of Eq. (11701 ) shows that their nontrivial solutions 
appear only when the coupling parameter a exceeds a certain critical value a* . This corresponds to the situation when the 
stability criterion |41] cr(x) - (yri^) — 3'f(-^))/(1 - x^) > calculated with the Fermi distribution, yrix) - x^ + Y{x; 0, 1), 
fails in a certain point < x, < 1 within the Fermi sphere: cr(x,) — > 0. There are two different types of such instabilities 
depending on the screening parameter xq (Fig. |65]) . For xq below a certain threshold value Xth ~ 0.32365 (weak screening 
regime, WSR) the instability point x, sets rather close to the Fermi surface: 1 - x, <K 1, while it drops abruptly to 
zero at xq — > Xth and equals zero for all xq > Xth (strong screening regime, SSR). The critical coupling (^'(xo) results 
in a monotonously growing function of xq with the asymptotic a* (ln2/xo - 1)"' at xq — > and staying analytic at 
ttth = Q''(^th) ~ 0.91535, where it only exhibits an inflection point. 

These two types of instabilities give rise to different types of TT from the state 1 at a > a*: at SSR a void appears 
around x = (1 ^ 2 transition), and at WSR a gap opens around x, (1 — » 3 transition). Further analysis of Eq. ( 1170b 
shows that the point Xth, ffth represents a triple point in the phase diagram in the variables xq, a (Fig. l65T l where the 
phases 1, 2, and 3 meet one another Similarly to the onset of instability in the Fermi state 1, each evolution of TT to 
higher order phases with growing a is manifested by a zero of cr(x), Eq. ( 11711 ). at some point < x, < X2„ different 
from the existing interfaces. If this occurs at the very origin, x, = 0, the phase number rises at TT by 1, corresponding 
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Figure 65: Instability point Xj and critical coupling a* as functions of screening. The regions of weak screening (WSR) and strong screening (SSR) are 
separated by the threshold value .v,;, . Note that a x,i,, a,i, is the triple point between the phases 1, 2, 3 in Fig. 1661 

to the opening of a void (passing from odd to even phase) or to emerging "island" (even — » odd). For x, > 0, either 
a thin spherical gap opens within a filled region or a thin filled spherical sheet emerges within a gap, so that the phase 
number rises by 2, living the parity unaltered. A part of the whole diagram shown in Fig. |66] demonstrates that with 
decreasing of xq (screening weakening) all even phases terminate at certain triple points. This, in particular, agrees with 
numerical studies of the considered model along the line jcq = 0.07 at growing a |137], where only the sequence of 
odd phases 1 — > 3 — > 5 — > . . . has been revealed (shown by the arrow in Fig. |65] |. The energy gain A (r) at TT as 
a function of small parameter t = a/a* - 1 is evidently proportional to t times the volume of a new emerging phase 
region (empty or filled). Introducing a void radius S and expanding the energy gain A{S) - z{n{x, S)} — zinpix)] in 6, one 
gets A = —I3\t6^ + ySa^^ + 0(6^), /3i,/32 > 0. As a result, the optimum void radius is 5 ~ a/t. Con sequ ently we have 



A(t) ~ T^^^ indicating a resemblance to the known "^/2-kind" phase transitions in the theory of metals 12371]. The peculiar 
feature of our situation is that the new segment of the Fermi surface opens at very small momentum values, which can 
dramatically change the system response to, e.g., electron-phonon interaction. On the other hand, this segment may have 
a pronounced effect on the thermodynamic al properties of ^He at low temperatures, especially in the case of P-pairing, 
producing excitations with extremely small momenta. 

For a TT with unchanged parity, the width of a gap (or a sheet) is found to be ~ t so that the energy gain is A (t) ~ 
and such TT can be related to the second kind. It follows from the above consideration that each triple point in the 
xo - a phase diagram is a point of confl uence of t wo ^/2-kind TT lines into one 2-kind line. The latter type of TT has 



already been discussed in the literature II135L 113711 . Here we only mention that its occurrence on a whole continuous 
surface in the momentum space is rather specific for systems with strong fermion-fermion interaction, while the known 
TT's in metals, under the effects of crystalline field, occur typically at separate points in the quasi-momentum space. It is 
interesting to note that in the limit xq — > 0, a — » 0, reached along a line a - kxo, we attain the exactly solvable model: 
U(p) — > (T.Tif'UQSip) with t/() = k/{2Mpf), which is known to display FC at all Uo > ll4lll . The analytic mechanism of 
this behavior is the disappearance of the poles of U (p), Eq. (I172l i. as po — > 0, restoring the analytical properties necessary 
forFC. Otherwise, the FC regime corresponds to the phase order 00, whe n the density of infinitely thin filled (separated 



by empty) regions approaches some continuous function < v(x) < 1 1113711 and the dispersion law turns flat according to 
Eq. (I170I I. A few remarks should be made at this point. 

First, the considered model formally treats xq and a as independent parameters, though in fact a certain relation 
between them can be imposed. Under such restriction, the system ground state should depend on a single parameter, 
say the particle density p^^, along a certain trajectory q'(xo) in the above suggested phase diagram. For instance, with the 
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Figure 66: Phase diagram in variables "screening-coupling". Each phase with certain topology is labeled by the total number of filled and empty regions 
(see Fig. I64t . Even phases (shaded) are separated from odd ones by "^/2-kind" topological transition (TT) lines, while 2-kind TT lines separate odd 
phases. Triple points, where two '/2-TT and one 2-TT meet, are shown by circles. 



simplest Thomas-Fermi relation for a free electron gas a(x(yi = this trajectory stays fully within the Fermi state 1 
over all the physically reasonable range of densities. Hence a faster growth of a(xo) is necessary for realization of TT in 
any fermionic system with the interaction, Eq. ( I172l i. 

Second, at increasing temperatures, the stepwise form of the quasiparticle distribution is melting. Therefore, as 
temperature moves away from its zero value, the concentric Fermi spheres are taken up by FC. In fact, these arguments do 
not work in the case of a few icebergs. Thus, it is quite possible to observe the two separate Fermi sphere regimes related 
to the FC and iceberg states. 

Ther e is a goo d reason to mention that neither in the FC phase nor in the other TT phases, the standard Kohn-Sham 
scheme 124011 is no longer valid. This is because in the systems with FC or TT phase transitions the occupation 
numbers of quasiparticles are indeed variational parameters. Thus, to get a reasonable description of the system, one has 
to consider the ground state energy as a functional of the occupation numbers E[(n{p)] rather than a functional of the 
density ^[pj IsllsellSTl]. 



16. Conclusions 

In this review, we have described the effect of FCQPT on the properties of various Fermi systems and presented 
substantial evidence in favor of the existence of such a transition. We have demonstrated that FCQPT supporting the 
extended quasiparticle paradigm forms strongly correlated Fermi systems with their unique NFL behavior. Vast body of 
experimental facts gathered in studies of various materials, such as high-r^ superconductors, heavy-fermion metals, and 
correlated 2D Fermi liquids, can be explained by a theory based on the concept of FCQPT. 

We have established that the physics of systems with heavy fermions is determined by the extended quasiparticle 
paradigm. In contrast to the stated Landau paradigm that the quasiparticle effective mass is a constant, within the extended 
quasiparticle paradigm the effective mass of new quasiparticles strongly depends on the temperature, magnetic field, 
pressure, and other parameters. The quasiparticles and the order parameter are well defined and can be used to describe 
the scaling behavior of the thermodynamic, relaxation and transport properties of high-Tc superconductors, HF metals, 
2D electronic and ^He systems and other correlated Fermi systems. The quasiparticle system determines the conservation 
of the Kadowaki- Woods relation and the restoration of the LFL behavior under the application of magnetic fields. 
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We have also shown both analytically and using arguments based entirely on the experimental grounds that the data 
collected on very different strongly correlated Fermi systems reveal their universal scaling behavior. This is because all 
above experimental quantities are indeed proportional to the normalized effective mass exhibiting the scaling behavior. 
Since the effective mass determines the thermodynamic, transport and relaxation properties, we conclude again that HF 
metals placed near their QCP demonstrate the same scaling behavior, independent of the details of HF metals such as their 
lattice structure, magnetic ground state, dimensionality etc. In other words, materials with strongly correlated fermions 
can unexpectedly be uniform despite their prominent diversity. 

We have also shown that in finite magnetic fields, in the NFL region formed by FCQPT the curve of any second order 
phase transitions passes into a curve of the first order ones at the tricritical point leading to the violation of the critical 
universaUty of the fluctuation theory. This change is generated by the temperature-independent term of the entropy formed 
behind FCQPT. The quantum and thermal critical fluctuations corresponding to second-order phase transitions disappear 
and have no effect on the behavior of the system at low temperatures, and the low temperature thermodynamics of heavy- 
fermion metals is determined by quasiparticles. 

We have found that the differential conductivity between a metal point contact and a HF metal or a high-Tc super- 
conductor can be highly asymmetric as a function of the applied voltage. This asymmetry is observed when a strongly 
correlated metal is in its normal or superconducting state. We have shown that the application of magnetic field restoring 
the LFL behavior suppresses the asymmetry. Correspondingly, we conclude that the particle -hole symmetry is macro- 
scopically broken in the absence of applied magnetic fields, while the application of magnetic fields restores both the 
LFL state and the particle-hole synnmetry. The above features determine the universal behavior of strongly correlated 
Fermi systems and are related to the anomalous low-temperature behavior of the entropy, which contains the temperature 
independent term. 

In the future, the realm of problems should be broadened and certain efforts should be made to describe the other 
macroscopic features of FCQPT, such as the propagation of zero-sound, sonic and shock waves. In addition to the already 
known materials whose properties not only provide information on the existence of FC but also almost cry aloud for 
such a condensate, there are other materials of enormous interest which could serve as possible objects for studying the 
phase transition in question. Among such objects are neutron stars, atomic clusters, ultra cold gases in traps, nuclei, and 
quark plasma. Another possible area of research is related to the structure of the nucleon, in which the entire "sea" of 
non- valence quarks may be in FCQPT. The combination of quarks and the gluons that hold them together is especially 
interesting because gluons, quite possibly, can be in the gluon-condensate phase, which could be qualitatively similar to 
the pion condensate proposed by A.B. Migdal long ago. We believe that FC can be observed in traps, where there is 
the possibiUty of controlling the emergence of a quantum phase transition accompanied by the formation of FCQPT by 
changing the particle number density. 

Overall, the ideas associated with a new phase transition in one area of research stimulates intensive studies of the 
possible manifestation of such a transition in other areas. This has happened in the case of metal superconductivity, 
whose ideas were successfully used in describing atomic nuclei and in a possible explanation of the origin of the mass of 
elementary particles. This, quite possibly, could be the case with FCQPT. 

Finally, our general discussion shows that FCQPT develops unexpectedly simple, yet completely good description 
of the NFL behavior of strongly correlated Fermi systems, while the extended quasiparticle paradigm constitute the 
properties inherent in strongly correlated fermion systems. Moreover, the extended paradigm can be considered as the 
universal reason for the NFL behavior observed in various HF metals, liquids, and other Fermi systems. 
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